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ABSTRACT 


The  problem  of  scheduling  a  fairly  large  number  of  jobs  through  a 
production  facility  containing  a  small  number  of  identical  processors 
is  explored.  The  application  of  an  analytic  procedure  is  studied,  but 
it  is  shown  to  produce  computational  problems  which  are  beyond  the  scope 
of  the  existing  computer  and  are,  therefore,  not  practical. 

Heuristic  scheduling  rules  are  then  developed  and  their  relative 
merits  tested  using  several  measures  of  effectiveness.  One  rule  is 
shown  to  be  best  according  to  each  measure  of  effectiveness.  The  result 
is,  knowing  the  quantity  he  wants  to  optimize,  one  can  choose  a  rule 
and  obtain  a  schedule  which  will  give  a  near  optimum  according  to  that 
measure . 

Algorithms  were  developed  for  testing  and  implementing  the  rules. 
The  time  shared  APL  system  at  the  University  of  Alberta,  which  permits 
remote  access  to  an  IBM  360/67  computer,  was  used  in  the  simulations. 
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CHAPTER  I 


THE  SCHEDULING  PROBLEM 

1 . 1  Definition  of  the  Scheduling  Problem 

The  scheduling  problem  is  to  allocate  a  set  of  jobs  to  a  shop 
or  production  facility  in  such  a  way  as  to  minimize  or  maximize  some 
measure  of  effectiveness.  For  our  purposes,  a  shop  will  consist  of  an 
integral  number  of  units  of  capacity  which  will  be  available  for  a  finite 
time.  Time  will  also  be  expressed  as  an  integral  number  of  units.  Thus, 
a  shop  will  be  specified  by  a  vector  C  where  the  N'th  element  C[N_[*  is 
the  capacity  available  in  time  period  N. 

A  job  is  defined  by  two  parameters,  capacity  and  duration . 

The  capacity  requirements  of  the  J'th  job,  X [j*l,  is  the  number  of  units 

of  shop  capacity  required  to  process  the  job  and  is  expressed  in  the 

same  units  as  the  capacity  of  the  shop.  The  dura t ion  of  the  J'th  job, 
p[j],  is  the  number  of  units  of  time  required  to  process  the  job  and  is 
expressed  in  the  same  units  of  time  used  to  define  the  shop.  Both 
capacity  and  duration  are  integers.  We  define  the  s_ize_  of  the  job  to 
be  the  product  of  its  capacity  and  duration  or  X  [j^  ’  P  [j  _j. 

The  set  of  jobs  to  be  processed,  or  job  roster,  is  defined  by 

the  two  vectors,  capacity  X  and  duration  P.  The  J’th  job  is  defined  by 

the  J'th  element  of  each  of  these  vectors.  No  jobs  are  added  to  the 
roster  while  the  shop  is  processing  the  jobs  it  originally  contained. 

A  schedule  is  defined  by  a  vector  NUM  where  NUM[jJis  the 
starting  time  for  job  J.  It  is  the  allocation  of  a  roster  of  jobs  through 
the  shop.  The  schedule  must  be  feasible ;  that  is,  in  any  time  period, 


*  The  APL  notation  used  throughout  is  described  in  Section  4.1.3. 
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the  total  capacities  of  the  jobs  in  process  must  not  exceed  the  capacity 
available  in  the  shop  in  that  period.  The  scheduling  process  is  not  con¬ 
strained  by  limitations  on  the  number  of  time  periods  the  shop  is  in 
operation . 

If  a  job  is  scheduled  to  occupy  a  unit  of  capacity  in  a  given 
time  period,  it  will  occupy  the  complete  unit  of  capacity  for  the  complete 
time  period.  The  model  may  be  visualized  as  a  shop  containing  machines 
which  can  do  only  one  job  at  a  time  and  can  only  be  set  up  for.  different 
jobs  at  certain  intervals  (e.g.,  each  morning)  making  it  impossible  to 
change  jobs  in  the  time  period  (during  the  day). 

A  measure  of  effectiveness  is  a  criteria  for  evaluating  a 
schedule.  It  is  a  mathematical  formulation  of  some  characteristic  of 
the  schedule  which  is  considered  to  be  important  and  becomes  the  objective 
function  in  the  analytic  approach  to  the  problem.  For  example,  if  it  is 
important  to  process  many  jobs  as  soon  as  possible,  mean  start  time  would 
be  a  suitable  measure  of  effectiveness. 

The  following  example  shows  the  problem  and  the  scheduling 
process  for  a  simple  case. 

1 . 2  Example  of  the  Scheduling  Process 

Example  Problem  1  below  illustrates  the  allocation  of  a  roster 
of  three  jobs  to  a  given  production  facility  or  shop.  The  shop  and 
rosters  to  be  scheduled  are  first  defined. 

The  shop  is  to  have  five  units  of  capacity  and  run  for  25 
time  periods.  It  is  represented  by  the  vector  C=5,5,...,5  or  25^5.  The 
jobs  are  defined  by  two  vectors: 


X  =  2,4,3 


P  =  1,3,6 


’  i 
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This  defines  three  jobs;  the  first  occupies  two  units  of  capacity  for 
one  unit  of  time  represented  by 

X[J>  2,  P[j>  1  J  =  1 

The  second  and  third  jobs  occupy  four  capacity  units  for  three  time  units 
and  three  capacity  units  for  six  time  units,  respectively,  and  are  repre- 


sented  as  follows: 

X[j]=  4, 

P[j>  3 

J  = 

X[J>  3, 

'-d 

c? 

1 _ 1 

II 

ON 

J  = 

Figure  1  shows  the  shop  and  the  allocation  of  the  jobs  to  it  on  a  "first- 
come-fir st -served"  basis. 

Figure  1 

First-Come-First-Served  Assignment  of  Example  Problem  1 


24  25 


The  schedule  is  defined  by  the  vector  NUM  which  has  the  following 


value: 


NUM  =  1,2,5 

This  means  job  1  starts  in  period  1,  job  2  starts  in  period  2,  and  job  3 
starts  in  period  5.  The  finish  times  are  related  to  the  start  times  and 
are.  represented  by  the  vector  F: 


F  -  1,4,10 


The  number  of  units  of  time  the  shop  is  occupied  is  defined  as 


Max  FCjJor  10  in  this  case.  The  unused  shop  capacity  after  the 
U  J^Jmax 

schedule  is  determined  is  C[XJ  l$N£Max.  f[\j!1 

The  number  of  jobs  in  process  is  represented  by  a  vector  JOB 
where  JOB  N  is  the  number  of  jobs  in  process  in  time  period  N.  For  this 
case 

JOB  DO  =  1,  UN  £10 

0,  N-S10 


The  sum  of  start  times  is  J  =  /  NUM  J  or  +  /NUM,  that  is, 

J=1 


14-2+5=8  for  this  schedule.  The  average  start  time  is  this  quanitity 
divided  by  Jmax  or  8/3. 

Another  measure  of  effectiveness  is  based  on  the  idea  that  the 
profit  from  a  job  is  proportional  to  its  size.  A  delayed  profit  is  then 
proportional  to  job  size  multiplied  by  delay  or  units  of  elapsed  time 
before  processing  begins.  This  is  also  called  a  weighted  disutility 
factor  and  is  described  subsequently  in  3.2.1.f. 

For  this  schedule,  the  disutility  factor  is  calculated  below. 


Job 

Delay 

Size 

Disutility  Factor 

1 

0 

22 

0 

2 

1 

12 

12 

3 

4 

18 

72 

Total  delayed  profit  or  disutility  factor  is  0  +  12  +  72  =  84. 

More  complex  rules  may  be  used  to  schedule  the  jobs  resulting 
in  different  schedules.  The  calculated  value  of  each  measure  of 


5 


effectiveness  will  then  be  different  than  that  calculated  above;  in  fact, 
the  numerical  value  of  each  measure  of  effectiveness  is  a  function  of  the 
schedule  and  therefore  depends  on  the  scheduling  rule.  A  heuristic  rule 
may  be  used:  take  the  largest  remaining  job  first.  Figure  2  shows  the 
schedule  resulting  from  applying  this  heuristic  rule  to  the  roster  and 
shop  in  the  example  problem. 


Figure  2 

Largest  Remaining  Job  Scheduled  First  -  Example  Problem  1 


This  schedule  is  defined  by  the  start  time  vector 
NUM  =  10,7,1 

and  the  related  finish  time  vector  is 
F  =  10,9,6 

The  three  measures  of  effectiveness  are  mean  start  time  which 
is  (+/NUM)  7  3  =  (10  +  7  +  1)  73-6,  shop  occupied  time  =  Max  F 
which  is  F  =  10  and  a  delayed  profit  of  (9  x  2)  +  (6  x  12)  +  (0  x  18)  =  90. 

The  next  heuristic  rule  applied  to  this  example  problem  is 
similar  to  the  above  rule  except  that  the  largest  capacity  job  will  be 
scheduled  first,  and  an  attempt  will  be  made  to  fill  empty  spaces  with 
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smaller  capacity  jobs.  This  resulting  schedule  is  represented  by 
Figure  3  and  the  vectors: 

NUM  =  4,1,4 
F  =  4,9,3 


Figure  3 

Largest  Capacity  First  and  Fill  Empty  Spaces  -  Example  1 


The  measures  of  effectiveness  are: 

Average  start  time  =  (+/NUM)  7  3  =  (4  +  1  +4)  7  3  =3 

Shop  occupied  =  max  F  =9 

Delayed  profit  =  (2  x  3)  +  (0  x  12)  +  (3  x  18)  =  60 

The  values  of  each  measure  of  effectiveness  for  each  scheduling 


rule  are  summarized  in  Table  1  below 
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TABLE  1 

Summary  of  Results  -  Example  Pro b 1 era  1 


Schedule 

1 

2 

3 


Measure  of  Effectiveness 
Shop  Occupied  Mean  Start  Time  Delayed  Profit 

10  8/3  84 
10  6  90 
9  3  60 


Schedule  2  is  a  poor  arrangement  according  to  these  measures 
of  effectiveness  as  it  yields  the  poorest  value  of  all  the  schedules  for 
each  measure  of  ef f ect ivenss .  However,  it  is  entirely  possible  that 
another  measure  of  effectiveness  is  of  importance;  and  Schedule  2  ranks 
high  according  to  it. 

The  example  illustrates  a  conflict  between  the  measures  of 
effectiveness.  Considering  Schedules  1  and  3,  Schedule  1  gives  a  better 
average  start  time  while  Schedule  3  gives  a  better  delayed  profit  loss 
and  reduces  the  time  the  roster  occupies  the  shop.  Either  Schedule  1 
or  3  will  therefore  be  considered  best  depending  on  whether  average  start 
time  or  delayed  profit  loss  is  considered  by  the  scheduler  to  be  more 
important . 


1 . 3  Economic  Significance 

Our  objective  is  to  develop  a  practical  scheduling  technique 
which  may  be  used  in  a  job  shop.  To  be  practical,  the  technique  must 
use  equipment  and  time  only  in  such  quantities  as  are  justified  by  the 
savings  resulting  from  using  it.  Also,  the  technique  must  be  fast 
enough  so  that  the  information  gathered  from  its  use  is  available  in 


time  to  be  usable. 
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It  must  also  be  possible  for  the  scheduler  to  choose  the 
quanitity  that  is  important  to  him  as  a  measure  of  effectiveness.  This 
requires  that  important  measures  of  effectiveness  must  be  expressed  in 
formula  form  if  an  analytic  technique  is  chosen;  or  if  heuristic  methods 
are  used,  a  particular  heuristic  rule  in  a  set  of  rules  must  be  shown 
to  be  associated  with  certain  behavior  measured  by  each  measure  of 
effectiveness . 

A  job  shop  typically  consists  of  more  than  one  department, 
and  each  job  travels  through  the  shop  visiting  each  department  in  a  given 
sequence.  The  sequence  of  movements  between  departments  and  time  spent 
in  each  department  is  usually  different  for  each  job. 

The  current  study  focuses  on  one  department .  There  is  a  queue 
or  roster  of  jobs  waiting  for  service  at  that  department,  and  no  new 
jobs  are  added  to  the  queue  until  the  entire  queue  or  roster  has  been 
processed  by  the  department.  The  department  can  handle  more  than  one 
job  at  a  time  as  long  as  the  combined  capacity  requirements  of  all  jobs 
running  do  not  exceed  the  capacity  of  the  department. 

While  this  is  a  simple  problem,  it  is  the  key  to  the  entire 
job  shop  scheduling  problem;  that  is,  the  individual  department. 

Behavior  within  the  department  must  be  understood  before  we  can 
analyze  flows  between  the  departments.  This  study  gives  some  results 
for  a  static  problem  where  no  new  jobs  enter  the  roster  during  processing; 
however,  the  model  developed  is  sufficiently  general  that  a  dynamic 
problem  involving  arrivals  to  the  queue  could  be  studied.  The  conclusions 
of  this  study  should  provide  some  insight  into  what  would  happen  in  a 
dynamic  case. 


This  study,  therefore,  is  a  start  on  the  solution  of  the 
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problem  of  allocating  jobs  with  fixed  capacity  and  duration  requirements 
through  a  shop  with  fixed  capacity  in  each  department.  The  time  period 
and  unit  of  capacity  which  the  job  must  occupy  entirely  or  not  occupy 
at  all  is  large  relative  to  the  size  of  the  job  meaning  a  good  "fit" 
between  the  jobs  in  the  shop  is  required  to  avoid  losing  large  amounts 
of  machine  time. 

This  model  is  reasonable  as  we  could  visualize  a  department 
containing  a  few  high  priced  identical  machines  which  require  considerable 
set-up  time  for  each  job.  Each  job  requires  one  or  more  machines  for 
processing.  A  problem  where  there  are  many  facilities  and  small  set-up 
time  has  more  of  a  flow  or  assembly  line  nature  and  is  not  really  a  job 
shop  in  the  sense  discussed  here. 

1 *4  Related  Problems 

A  model  often  analyzed  in  the  literature  consists  of  a  job 
shop  containing  a  number  of  different  machines  or  departments  and  jobs 
which  follow  a  predetermined  sequence  of  steps  through  the  shop.  [j3,4,5] 
The  jobs  are  defined  by  one  dimension,  duration,  and  occupy  the  entire 
department  when  they  are  in  it . 

This  is  in  contrast  to  the  model  discussed  above  where  the 
entity  studied  is  one  department.  In  many  models,  the  department  is 
indivisible;  in  this  case,  it  contains  nine  identical  units  of  capacity. 
Each  job  occupies  one  or  more  units  of  department  capacity  for  a  given 
number  of  units  of  time. 

Analysis  of  the  scheduling  problem  using  both  analytic  and 
heuristic  methods  has  been  presented  in  the  literature.  Due  to  analytic 
methods  becoming  impractical  because  of  the  volume  of  computational 
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requirements,  major  attempts  have  been  made  to  develop  heuristic 
methods  which  give  good  approximations  to  optima  for  certain  measures 
of  effectiveness. 

An  analytic  approach  to  a  job  shop  scheduling  problem  is 
presented  by  Brooks  and  White.  C3U  The  problem  consists  of  several 
machines;  each  job  is  routed  around  the  machines  in  a  given  sequence. 
In  the  solution,  a  matrix  A  is  developed  which  is  equivalent  to  the  A 
matrix  in  the  standard  zero  -  one  integer  programming  problem  form: 

min  cx 

subject  to  Ax  b  ;  x  =  0  or  1 

To  develop  this  matrix,  two  other  matrices  are  developed  and 

combined: 

1.  A  sequencing  array  consisting  of  the  routing  order  of 
each  job  around  the  machines. 

2.  A  facility  array  consisting  of  the  time  each  job  spends 
on  each  machine. 

The  problem  array,  or  A  matrix,  indicates  the  minimum  completion  time 
for  each  job  on  each  machine  and  is  derived  from  the  above  two  arrays. 

The  integer  programming  formulation  could  be  completed  by 
expressing  the  objective  function  or  measure  of  effectiveness  as  the  c 
vector  and  introducing  capacity  constraints  in  the  b  vector. 

However,  the  author  states  that  this  formulation  becomes  too 
massive  if  a  problem  of  practical  size  is  considered  and  could  not  be 
implemented  on  existing  computers.  Ten  small  problems  were  done  using 
this  method,  giving  an  optimum  value  for  three  different  objective 
functions,  and  the  results  are  then  compared  with  the  results  from 
application  of  three  heuristic  scheduling  methods. 


•t 
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The  three  heuristic  rules  are: 

1.  Random  selection 

2.  Shortest  operation  first 

3.  Longest  remaining  time  to  completion  first 

The  three  measures  of  effectiveness  are: 

A.  Lateness 

B.  Machine  idle  time 

C.  Total  time  to  complete  all  orders 

Table  2  shows  the  comparison  of  the  merits  of  using  each  of  the 
four  scheduling  methods  as  indicated  by  each  measure  of  effectiveness. 

A  low  number  always  indicates  a  better  result  than  a  high  one,  and  the 
numbers  represent  dimensioned  quantities  rather  than  dimensionless 
ranks.  ’  A  ranking  method  is  described  in  Section  3.4. 

TABLE  2 

Comparison  of  Brooks  and  White  Results 


Rule 

Measure 

of  Effectiveness 

A 

B 

C 

Integer  Programming 

4836  '* 

3333 

1180 

i 

Heuristic  1 

8377 

5586 

1473 

Heuristic  2 

7108 

5521 

1422 

Heuristic  3 

9531 

4911 

1323 

This  result  shows  that  for  at  least  some  measures  of  effective¬ 
ness,  heuristic  rules  can  be  developed  which  produce  results  not  far 
from  optimum. 

Using  a  rather  extensive  model,  Conway  has  studied  the  effect 
of  a  large  number  of  heuristic  rules  measured  by  several  measures  of 
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effect iveness .  [4]  The  situation  is  similar  to  that  studied  by  Brooks 
and  White  consisting  of  nine  independent  machines.  The  jobs  need  not 
be  processed  on  all  machines  and  may  return  to  a  given  machine  after  a 
subsequent  operation  on  another  machine.  The  number  of  operations  per 
job  is  a  random  variable  with  mean  9,  and  there  is  equal  likelihood  of 
the  job  going  to  any  other  operation  after  any  one  is  completed.  The 
time  to  perform  each  operation  on  each  job  is  exponentially  distributed 
with  mean  1 . 


The  scheduling 
processed  by  assigning  pr 
or  linear  combinations  of 


1. 

2. 

3. 

4. 

5. 

6 . 

The  measures 


Random  sele 
F irst -come- 
Shortest  or 
Greatest  or 
Largest  or 
Queue  which 

a.  Shortes 

b .  Queue  w 
of  effective 


rules  were  programmed  and  large  rosters  were 
iority  on  the  basis  of  the  following  criteria 
these  criteria, 
ct ion . 

first-served . 

longest  next  process  time, 
smallest  number  of  remaining  operations, 
smallest  total  work  required. 

job  enters  for  next  operation: 
t  queue. 

ith  least  work  required. 

ness  used  are  indicators  of  work  in  process 


inventory  and  were: 

1.  Work  remaining  and  not  completed. 

2.  Work  in  process  in  the  shop. 

3.  Work  content  in  queues. 

It  was  observed  that  there  are  important  differences  between  rules  that 
are  simple  and  easily  implementab le ,  and  different  rules  give  different 
performances  according  to  each  measure  of  effectiveness.  No  one  rule 


■ 
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is  best  according  to  all  measures  of  effectiveness. 


Fabrycky  and  Shamblyn  develop  and  test  a  particular  heuristic 
rule  with  a  model  similar  to  that  described  above . [5^  The  heuristic 


method  is  based  on  an  urgency  factor  which  is  the  probability  that  a 


job  will  be  completed  by  its  due  date  given  the  current  conditions. 


The  total  expected  time  for  the  iJth  order  to  complete  its 

routing  on  n  machines  is 

n 


where  tj  is  the  time  on  the  j  machine  and  consists  of 


representing  the  move  time,  queue  time,  set-up  time,  and  process  time, 
respectively,  associated  with  the  j'th  machine. 


T^  is  a  random  variable  whose  value  is  affected  by  the  multitude 


of  factors  affecting  the  values  of  the  individual  m,  q,  s,  and  p  and  the 
number  of  these  quantities  making  up  t^  for  any  particular  job.  By 
the  Central  Limit  Theorem,  we  can  say  T^  has  an  approximately  normal 
distribution  with  fixed  mean  and  standard  deviation. 

Since  T^  has  a  normal  distribution  at  each  point  in  time,  a 
normalized  Z  factor  can  be  calculated  for  each  job. 


n 


(Di  -  c)  -  ^cr j 


where  =  due  date  of  i'th  job 


C  =  current  date 


n 


•j 


\  =  mean  time  to  completion 
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n 


Z  Mi  =  variance  of  the  time  to  completion. 
j=k  ' 

Therefore,  Z^  is  a  measure  of  the  probability  that  the  i 1 th  order  will 
be  completed  by  its  due  date. 

The  algorithm  consists  of  computing  Z^  for  all  jobs  in  each 
time  period  and  scheduling  those  with  the  least  probability  of  completion 
or  greatest  urgency  factor.  The  urgency  factor  of  the  i'th  job  is  the 
calculated  above.  The  measure  of  effectiveness  is  the  percent  of 
orders  completed  late,  and  the  algorithm  reduced  this  by  22  percent 
compared  to  a  first-come-first -served  routine. 


1 . 5  Summary  of  Results  in  This  Thesis 

Fourteen  scheduling  rules  were  developed  and  their  effect 
studied  according  to  four  measures  of  effectiveness.  The  scheduling 
rules  were  applied  to  sixty  rosters  of  twenty-five  jobs  each.  Half  of 
the  rosters  were  generated  by  a  process  which  produces  numbers  uniformly 
distributed  between  two  bounds,  and  the  other  half  were  generated  by 
a  Poisson  process  with  fixed  mean.  Certain  undesirable  results  were 
then  eliminated  from  the  Poisson  rosters  as  described  in  Section  3.1. 


1.5.1  The  Scheduling  Rules 

The  scheduling  rules  assign  priority  to  job  charac¬ 
teristics,  that  is  according  to  largest  or  smallest  capacity,  duration, 
or  size.  In  order  to  form  a  unique  priority  order,  two  assign¬ 
ments,  called  a  primary  sort  and  tie  breaking  or  secondary  sort, 
must  be  made.  Some  of  the  rules  involve  only  one  sort  and  are 
called  primary  scheduling  rules.  These  assign  priority  according 
to  one  job  characteristic  only,  and  ties  are  broken  in  a  random 


■ 
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manner.  The  rules  which  form  a  unique  priority  order  are  called 
secondary  scheduling  rules. 

Except  for  the  first  scheduling  rule  in  Table  3,  a 
fitting  procedure  was  implemented.  If  at  any  time  the  occupied 
capacity  is  less  than  shop  capacity  and  occupied  capacity  plus 
capacity  of  the  highest  priority  unassigned  job  exceeds  shop 
capacity,  a  search  is  made  through  the  remaining  unassigned  jobs 
in  decending  priority  order  for  the  first  job  that  will  not  violate 
shop  capacity.  The  first  scheduling  rule  assigns  on  a  first-come- 
first-served  basis. 

Table  3  shows  the  names  of  the  scheduling  rules  and 
the  priorities  associated  with  each. 

TABLE  3 

Names  of  S cheduling  Rules 


Rule  Name 

Priorities 

Primary  Priority 

Tie-Breaking  Priority 

RANDOM 

first-come  first-served. 

RANDOMFIT 

first  job 

in  roster 

SS 

smallest  size 

LS 

largest  size 

SP 

smallest  duration 

LP 

largest  duration 

SPSC 

smallest  duration 

smallest  capacity 

SPLC 

smallest  duration 

largest  capacity 

LPSC 

largest  duration 

smallest  capacity 

LPLC 

largest  duration 

largest  capacity 
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Rule  Name 


Prior it ies 

Primary  Priority  Tie -Breaking  Priority 


SSSP 

SSLP 

LSSP 

LSLP 


smallest  size 


smallest  size 


largest  size 
largest  size 


smallest  duration 


largest  duration 


smallest  duration 


largest  duration 


These  criteria  do  not  necessarily  determine  the  order  in  which  the 
jobs  appear  in  the  final  schedule.  They  determine  the  order  in 
which  an  attempt  is  made  to  schedule  the  jobs  subject  to  capacity 


constraints  and  the  fitting  procedure  described  above. 


1.5.2  The  Measures  of  Effectiveness 

The  measures  of  effectiveness  used  are  defined  as 
follows.  The  definitions  are  expanded  in  Section  3.2.1. 

1 .5 . 2  .a  Measure  A  or  mean  start  time.  This  is  the  sum  of 
the  start  times  divided  by  the  number  of  jobs.  Where  NUM 
is  the  vector  of  start  times  and  NUM  is  the  number  of 
elements  in  this  vector,  mean  start  time  is  (+/NUM)  f^NUM. 

1 . 5 . 2 . b  Measure  B  or  units  of  time  the  shop  is  occupied. 

This  is  the  maximum  element  of  F,  the.  vector  of  finish  times, 
and  is  denoted  by  f/F. 

1 . 5 . 2 .  c  Measure  C  or  average  efficiency.  This  is  the  portion 
of  available  time  used  in  the  second  and  third  quarters  of 
the  period  of  time  the  shop  is  occupied. 

1 . 5 ♦ 2 . a  Measure  D  or  disutility  factor.  This  is  the  size 
of  each  job  multiplied  by  che  time  it  is  delayed  before 
processing  begins  summed  over  the  whole  roster. 


■ 


<■ 


■ 
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1.5.3  Results  of  the  Study 

The  results  of  this  study  are  summarized  in  Table  4 
where  the  ranking  of  each  of  the  four  primary  scheduling  rules  as 
measured  by  each  measure  of  effectiveness  are  shown.  The  rankings 
are  based  on  the  test  sample  of  sixty  rosters.  A  rank  of  1 
indicates  the  best  result  and  a  rank  of  4  indicates  the  poorest 
result  indicated  by  the  measure  of  effectiveness. 

TABLE  4 

Summarv  of  Results 


Scheduling  Rule 


Measure  of  Effectiveness 
A  ’  B  C  D 


SS  2 

LS  4 

SP  1 

LP  3 

The  conclusion  drawn  from  this  table  is  that  if  we 
wish  to  minimize  mean  start  time,  wTe  give  priority  to  jobs  with 
the  smallest  duration  requirement.  To  clear  all  jobs  out  of  the 
shop  in  minimum  time,  maximize  efficiency,  or  minimize  the  delayed 
profit  factor,  we  give  priority  to  the  largest  size  job. 

We  have  also  found  that  certain  tie-breaking  sorts  will 
improve  the  schedule  as  indicated  by  these  measures  of  effectiveness. 
When  priority  is  given  to  the  smallest  duration  job  which,  is  necessary 
to  minimize  mean  start  time,  ties  should  be  broken  by  giving 
priority  to  the  smallest  capacity  (or  largest  size)  job  to  further 
improve  the  mean  start  time.  When  priority  is  given  to  the  largest 


4  4  4 
111 
3  3  3 
2  2  2 
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size  job,  the  other  three  measures  are  optimized  further  by  giving 
priority  to  the  smallest  duration  (or  largest  capacity)  job  in  case 
of  ties  in  the  size  priority  assignment.  In  Section  3.5.2,  it  is 
shown  that  these  relationships  are  statistically  significant  (0.05 
level  of  significance  or  better). 

The  capacity  of  the  shop  was  constant  throughout  the 
scheduling  period.  The  schedule  was  constrained  by  the  requirement 
that  the  sum  of  the  capacities  of  the  jobs  in  process  must  at  no  time 
exceed  the  capacity  of  the  shop;  however,  the  shop  was  allowed  to  run 
for  a  sufficient  amount  of  time  so  that  the  schedule  was  not  con¬ 
strained  by  the  amount  of  time  the  shop  is  running. 

1.5.4  Summary 

One  rule  in  a  set  of  simple  rules  based  on  roster  charac¬ 
teristics  (capacities,  durations  and  sizes  of  jobs)  has  been  shown  to 
produce  a  best  result  according  to  each  of  the  four  measures  of 
effectiveness.  Conflicts  regarding  which  scheduling  rule  is  best 
exist  as  a  different  rule  ’has  been  shown  to  be  best  depending  on  the 
measure  of  effectiveness  which  is  held  to  be  important. 

The  rule  which  has  shown  to  be  the  best  in  the  set  has 
not  been  shown  to  be  optimal  or  near  optimal;  however,  a  more  optim¬ 
al  rule  would  be  more  complex  and  difficult  to  apply  in  practice. 

It  is  felt  that  the  set  of  rules  is  reasonably  comprehensive,  and  one 
rule  in  a  comprehensive  set  of  rules  is  likely  to  produce  a  near 
optimum  according  to  any  measure  of  effectiveness.  An  indication  of 
the  ability  of  the  rules  to  produce  optima  is  given  by  mid-range 
efficiency,  as  the  optimal  is  known  to  be  1;  and  for  mainy  rosters, 
the  efficiency  produced  by  the  best  rule  was  1. 


i 


CHAPTER  2 


ANALYSIS  OF  THE  SCHEDULING  PROBLEM 

2 . 1  Formulation  as  an  Integer  Programming  Problem 

We  will  now  present  an  analytic  method  which  could  be  used 
to  solve  the  scheduling  problem.  To  formulate  it  as  an  integer 
programming  problem,  we  must  express  it  in  the  form: 

min  cy 

subject  to 

Ay^b  y  =  0  or  1 

This  is  the  standard  zero-one  integer  programming  form  except 

the  variable  x  is  replaced  with  y  to  avoid  confusion  with  x  denoting  the 

capacities  of  the  jobs. 

The  zero-one  variable  y  is  defined: 

yj  n  =  (1  if  job  j  starts  at  time  n 

(0  if  job  j  does  not  start  at  time  n 

The  first  part  of  the  A  matrix  represents  the  constraint  that  a  job 

starts  once  and  only  once: 

N 

y,  n  =  1  (j  =  1.2,  ....  J) 

n=l 

where  there  are  J  jobs  and  N  periods  of  time  available. 

The  remainder  of  the  constrain  matrix  A  specifies  that  the 
total  capacity  requirements  of  all  jobs  running  in  each  time  period  must 
not  exceed  the  shop  capacity  in  that  time  period.  The  constraint  is: 


yj>n 


x .  Z  c 
J  “  n 


(n  =  1,2, 


N) 


The  xj  corresponds  to  X  M  defined  in  Section  1.1  as  the 
capacity  requirement  of  the  j*th  job,  and  cn  is  the  C  [^N^  representing 
shop  capacity  in  period  n. 
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In  the  following  matrix,  thg  first  constraint  is  represented 
by  the  first  J  rows  and  the  second  constraint  by  the  last  N  rows.  This 
is  the  detailed  representation  of  the  form  Ay ^  b; 


r  1  0  ..  0  10. .0  10..  0 

yl,l 

l 

01. .0  01  . .0  01. .0 

y2,i 

l 

•  •  •  ••  •  •  •  0 

yj,i 

• 

•  •  •  ••  •  ••  • 

yl ,  2 

• 

00..  1  00..  1  00  ..  1 

y2,2 

l 

x^X2  . .  xj  00. .0  00. .0 

• 

L 

ci 

0  0  ..  0 

• 

c2 

x^X2  . .  xj  .  . 

•  •  • 

•  •  • 

yj,2 

• 

• 

• 

• 

• 

•  •  • 

x^x2..  xj 

• 

• 

• 

t 

—  _ 

— I 

c 

o 

_ 1 

yl,n 

y2,n 

• 

• 

• 

_yJ,n 

To  show  that  the  first  J  rows  of  this  array  represent  the 
first  constraint,  we  evaluate  the  first  row  eliminating  the  terms  with 
coefficient  0, 

n,i  +  ?i,2  +  •••  +  yi>n  =  1 

or  in  summation  form 
N 

2-  yj,n  =  1  ■  J  “  1 

n=l 

The  same  result  is  obtained  from  evaluating  the  J'th  row 

yJ,l  +  yJ,2  +  *•*  +  yJ,N  =  1 
and  in  summation  form 


1 


j  =  J 
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Similarly,  each  of  the  first  J  rows  of  the  matrix  represents  this 
relationship  with  j  equal  to  the  row  index. 

The  last  of  the  constraint  matrix  or  the  (J  +  1)  to  (J  +  N) 
rows  represent  the  capacity  constraints. 

Evaluating  the  first  row,  we  find 

x-^y^  +  *2^2  1  +  •••  +  xjyj  \  +  0  +  ...  +  0^-C]_ 

or  in  summation  form 


J 

>  x  •  y .  ^  c  n  =  1 

j  J  J>n  “  n 

In  the  second  and  subsequent  periods  of  time,  allowances  must 
be  made  for  jobs  which  have  been  started  in  an  earlier  time  period  and 
have  extended  into  the  current  period.  This  is  done  by  defining  an  Xj 
as  a  column  vector  of  elements  Xj  and  order  pj  .  x j ,  therefore,  extends 
downward  in  the  A  matrix  (direction  of  increasing  n  or  time)  through  a 
number  of  rows  equal  to  p j ,  the  job  duration  or  the  time  it  must  run 
once  it  is  started. 

The  objective  function  cy  of  the  integer  programming  formulation 
is  used  to  introduce  a  measure  of  effectiveness  into  the  problem.  The 
vector  c  represents  the  objective  function  as  shown  in  two  examples 
below . 


If  we  want  to  minimize  mean  start  time,  we  minimize  the  sum 


N 


I  Z 

1=1  n=l 


n 


yj,n 


or  its  equivalent  cy  in  matrix  notation  where 

c  =  Q  ,  1 ,  ...  1  2,2,  ...  2  n ,  n ,  .  .  .  n^| 

yT  =(7i,i  y2,i-*-yj,i  yi,2  y2,2***yj,2  yi,N  y2N-*-yji(] 
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If  we  want  to  minimize  the  measure  of  delayed  profit,  we 
introduce  the  size  of  the  jobs  in  the  formulation. 

J  N 

H  yj.n  n  sj 

3=1  n=l 

where  Sj  =  xj*Pj  or  the  size  of  the  j  '  th  job.  In  the  formulation  cy, 
c  becomes 

c  =[sl5  S2,...,sj,  2s^,  2s2  >  •  .  • ,  2s  j,  .  ,  .Ns^ ,  Ns2 ,  .  .  .  jNsj^J 
An  example  of  the  formulation  of  a  simple  problem  will  now  be  given 
before  a  method  of  solving  the  problem  is  described. 


2.1.3  Example 

Example  Problem  1,  Section  1.2,  is  formulated  as  an 
integer  programming  problem  using  this  method.  First  we  observe 
that  the  roster  can  be  completed  in  ten  units  of  time  even  with 
a  very  poor  schedule.  Therefore,  N>10  and  we  can  choose  N=10  for 
simplicity . 

For  the  constraint  or  A  matrix,  we  construct  the  first 
J  rows  noting  J=3  as  there  are  three  jobs  in  the  roster.  Since 
N-10,  ten  unit  matrices  of  order  three  arranged  horizontally  are 
required . 

To  construct  the  last  N=10  rows,  observe  that  for: 


Therefore,  the  extended  vectors  representing  job  duration  as  well 
as  capacity  become: 

—  T 

XI  =£2] 
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x2t  =[a>4>4] 

x3t  =[3,3, 3,3, 3,3] 


The 

constraint  matrix 

is  as 

follows 

where 

a  blank  (-) 

indicates 

the 

cell  would  be 

a  zero  regardless 

of  the 

particular 

problem . 

1-- 

1-- 

1-- 

1-- 

1-- 

1-- 

1-- 

1-- 

1-- 

1-- 

-1- 

-1- 

-1- 

-1- 

-1- 

-1- 

-1- 

-1- 

-1- 

-1- 

--1 

243 

043 

--1 

--1 

--1 

--1 

--1 

--1 

--1 

--1 

--1 

243 

-  -  - 

— 

— 

— 

— 

— 

— 

— 

043 

043 

243 

— 

— 

--- 

— 

— 

— 

--  - 

043 

043 

043 

243 

___ 

— 

— 

— 

— 

— 

003 

003 

043 

043 

243 

— 

— 

— 

— 

— 

003 

003 

003 

043 

047 

243 

— 

— 

— 

— 

000 

003 

003 

003 

043 

043 

243 

— 

— 

— 

000 

000 

003 

003 

003 

043 

043 

243 

— 

— 

000 

000 

000 

003 

003 

003 

043 

043 

243 

— 

000 

000 

000 

000 

003 

003 

003 

043 

043 

243 

000 

000 

000 

000 

000 

003 

003 

003 

043 

043 

The  b  vector  in  the  constraint  Ay<b  is 
bT  =  [1  1  1  5  5  5  5  5  5  5  5  5  5] 


meaning  that  the  shop  has  five  units  of  capacity  in  each  of  ten 
periods  of  time. 

The  specific  y  vector  of  zero  one  variables  is 

yT  =  [yi,i  +  y2,i  +  y3,i  +  •••  +  yi,io  +  ^2,10  +  ^3,10] 

The  constraint  Ay<b  therefore  becomes: 
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1)  for  the  one  start  only  constraint 


10 


(j  =  1,2,3) 


n=l 


2)  for  the  capacity  constraints 


2y 

3y 


1,1  +  4y2,l  +  3y3,l  ^  5 


(n  =  1) 

+  4?  n  + 

<5  (n  =  10) 

The  objective  functions  of  mean  start  time  or  delayed  profit 


3,5  +  3y3 , 6  +  3y3 , 7  +  4y2,8  +  3y3,8  '  "3,9 
3y3, 9  +  2yl,10  +  4y2 , 10  +  3y3 , 10 


can  be  introduced  into  the  formulation  using  the  c  vector.  For 
mean  start  time: 

c  =  [l,  1;  1,  2,  2,  2,  ...,  10,  10,  10] 
and  for  delayed  profit  where  S  =  2,  12,  18 

c  =  [2,  12,  18,  4,  12,  36,  ...,  20,  120,  18 6] 


2 . 2  Balas  Algorithm 

2.2.1  Description  of  the  Algorithm 

The  Balas  Algorithm  provides  a  possible  solution  to 
this  problem.  EJ  Any  problem  expressed  in  the  form 

min  cy 

subject  to  Ay<b  and  y  =  0  or  1 

is  a  zero-one  integer  programming  problem  and  is  solvable  using  an 
additive  algorithm  due  to  Balas.  Our  problem  has  been  expressed  in 
this  form. 

Since  there  are  a  finite  number  of  ways  a  fixed  number 
of  jobs  can  be  scheduled  through  a  shop  given  the  capacity 
restriction,  there  are  a  finite  number  of  feasible  solutions  to 
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this  problem.  These  are 
Ay<b .  There  may  also  be 
schedules  which  meet  the 
funct ion . 


the  solutions  satisfying  the  constraint 
one  or  more  optimal  solutions  respresenting 
constraints  and  minimize  the  objective 


A  zero-one  integer  programming  formulation  is  needed 
because  y  is  constrained  to  be  0  or  1,  not  any  positive  value  as  in 
linear  programming.  In  this  problem,  the  shop  capacity  is  divided 
into  a  small  number  of  indivisible  units  whose  capacity  is  important 
relative  to  the  total  shop  capacity.  For  this  reason,  a  variable 
representing  an  "either  on  or  off"  situation  rather  than  a  continuous 
flow  is  needed  representing  the  fact  that  a  job  is  either  completely 
in  or  out  of  a  capacity  unit . 

The  Balas  Algorithm  starts  by  setting  all  n  zero-one 
y  variables  to  zero  and  by  systematically  assigning  certain  variables 
the  value  1,  a  feasible  optimal  solution  is  obtained  or  evidence 
is  obtained  that  no  feasible  solution  exists.  This  is  accomplished 
by  trying  only  a  small  set  of  the  possible  2n  combinations  of  y  =  0 


or  1  . 


Figure  4  shows  the  solution  tree  for  the  Balas  additive 


algorithm  where  each  horizontal  plane  of  intersections  represents 
y^  y 2  ...  as  we  proceed  down  the  diagram. 


Figure  4 

Solution  Tree  -  Balas  Algorithm 
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If  the  line  connecting  any  variable  to  the  previous 
variable  or  the  start  point  approaches  the  plane  representing  the 
variable  from  the  northeast,  the  variable  is  assigned  the  value  0. 
Approach  from  the  northwest  indicates  assigning  the  variable  the 
value  1.  For  example,  branch  A  approaches  the  y-^  plane  from  the 
northeast  and,  therefore,  indicates  y^  =  0.  The  end  point  'a' 
represents  the  solution  y^  =  0,  y^  =  1,  y^  =  0.  All  eight  ways  of 
assigning  y-^,  y2,  y^  =  0  or  1  are  represented  by  the  eight  terminal 
points  on  the  last  plane  of  the  diagram. 

If  each  end  point  were  tested  for  feasibility  and  the 
optimum  solution  or  solutions  picked  out,  we  would  eventually  solve 
the  problem.  However,  this  would  be  time  consuming  on  a  large 
problem;  and  the  Balas  Algorithm  identifies  branches  that  if  we 
proceeded  down  to  the  end  of,  we  would  not  find  any  feasible 
solutions  or  a  feasible  solution  better  than  the  best  one  that  has 
already  been  found.  For  example,  if  neither  solution  a  or  ^  is 
feasible,  this  is  recognized  at  the  point  of  intersection  of  branch 
A  and  D  and  branch  D  is  "cut  off." 

If  a  feasible  solution  f  has  been  found  and  solutions 
(T  ,  and  T'are  also  feasible  but  not  as  good  as  f,  this  is 
recognized  at  the  intersection  of  branch  A  and  B  so  the  entire  right 
side  of  the  diagram  would  be  cut  off  and  not  tried.  The  only 
solutions  which  would  be  tried  are  f  and  g. 

Balas  has  shown  that  there  is  no  chance  that  a  feasible 
solution  better  than  the  one  already  obtained  will  be  cut  off.  The 
method,  therefore,  guarantees  an  optimal  feasible  solution  to  the 
zero-one  integer  programming  problem. 


y 
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2.2.2  Applicability  of  the  Algorithm 

One  of  the  factors  affecting  the  practicability  of  the 
use  of  the  Balas  Algorithm  for  this  problem  is  the  size  of  the  arrays 
required.  These  are: 

Vector  c  -  objective  function  -  length  JN 
Vector  b  -  constraint  vector  -  length  J  +  N 
Matrix  A  -  constraint  matrix  -  JN 

J  +  N  rows 

Therefore,  the  total  number  of  elements  required  to  formulate  the 
problem  in  integer  programming  form  is  J  +  N  +  JN  +  JN(J+N)  in 
addition  to  JN  zero-one  variables. 

For  the  problem  to  have  any  practical  significance 
and  to  have  different  scheduling  rules  produce  different  schedules 
for  comparison,  considerably  larger  rosters  than  the  one  of  size  3 
in  the  example  are  needed.  For  any  given  measure  of  effectiveness, 
it  is  evident  how  to  schedule  a  small  roster;  but  we  want  a  practical 
scheduling  technique  for  a  roster  where  the  optimum  is  not  obvious. 

Rosters  of  twenty-five  jobs  with  average  capacities 
one-third  to  one-half  of  shop  capacity  were  chosen  as  representative 
of  the  situation  we  wish  to  simulate.  A  roster  of  this  size  will 
occupy  the  shop  for  about  sixty  units  of  time  when  the  jobs  are 
assigned  on  a  first -come-first -served  basis. 

Therefore,  in  the  notation  above,  J  =  25  and  N  =  60. 

The  size  of  the  input  matrices  for  a  worthwhile  problem  becomes: 

c  -  objective  function  -  JN  =  25  x  60  =  1,500  cells 
b  -  constraint  vector  -  J  +  N  =  25  +  60  =  85  cells 
A  -  constraint  matrix  -  (JN)(J+N)  =  1,500  x  85  =  127,500 


cells 
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The  elements  in  these  cells  form  an  orderly  pattern,  and  it  is 
possible  to  program  the  computer  to  generate  these  arrays  given 
the  roster  and  the  shop.  However,  a  severe  storage  problem  in 
the  computer  develops  as  explained  in  Section  4.1.  In  addition, 
the  algorithm  is  very  slow  when  the  input  A  matrix  exceeds  size 
30  x  30. 

Because  of  the  storage  and  speed  problems  encountered 
even  when  a  high  speed  IBM  360/67  computer  is  used,  future  attention 
is  to  be  focused  on  heuristic  methods  of  solution. 

2 • 3  Heuristic  Methods 

The  heuristic  methods  adopted  are  based  on  roster  character¬ 
istics;  that  is,  size,  capacity,  or  duration  of  the  jobs  in  the  roster. 

The  heuristic  rules  can  be  shown  to  produce  close  to  optimal  results 
according  to  one  measure  of  effectiveness,  and  there  is  strong  indication 
that  they  are  close  to  optimal  according  to  the  others.  They  are 
practical  rules  of  thumb  and  are  easier  to  apply  than  the  method  outlined 
above . 

The  rules  developed  schedule  the  jobs  in  a  priority  sequence 
determined  from  one  or  two  of  the  three  job  characteristics:  size, 
capacity,  or  duration.  The  primary  scheduling  rules  do  not  define  a 
unique  schedule  as  they  consider  only  one  characteristic  while  the  secondary 
rules  break  ties  by  considering  one  of  the  two  remaining  characteristics. 

The  schedule  is  unique  as  the  three  characteristics  are  linearly  dependent 
and  two  have  been  used. 

Scheduling  is  done  in  a  two-step  procedure:  the  first  step 
assigning  the  priority  and  the  second  step  incorporating  the  shop 
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capacity  constraint .  The  first  step  is  a  sorting  procedure  where  the 
roster  is  rearranged  in  priority  order  with  the  highest  priority  job 
first  and  the  lowest  priority  job  last . 

The  second  step  consists  of  a  routine  which  assigns  the  jobs 
in  the  roster  from  start  to  finish  taking  account  the  shop  capacity 
constraint.  Two  types  of  routine  are  used;  the  first  assigns  the  jobs 
in  order  of  appearance,  and  the  second  searches  ahead  when  necessary  for 
a  job  which  does  not  violate  the  capacity  constraint  and  then  returns 
to  the  point  where  it  was  originally  stopped.  A  more  detailed  description 
of  this  operation  is  given  in  Section  4. 2. 2.  C.  III. 

Reference  to  Section  1.5.2  shows  that  the  primary  scheduling 
rules  assign  priorities  using  size  or  duration  only.  A  third  possible 
priority  rule  based  on  capacity  was  not  used  since  it  was  evident  early 
in  the  study  that  the  second  of  the  two  ways  of  performing  step  2  produces 
far  superior  results  (see  Section  4.2.2).  This  is  verified  by  the 
difference  in  performance  of  rules  RANDOM  and  RANDOMFIT  where  the  first 
uses  method  1  on  an  unsorted  roster,  and  the  second  uses  method  2  on  the 
same  unsorted  roster.  (See  Section  3. 5. 2. I,  Table  14.) 

The  capacity  based  priority  rule  already  arranges  the  jobs  in 
capacity  order  defeating  the  purpose  of  the  second  method  of  assignment. 
Therefore,  the  capacity  based  priority  system  which  defeats  the  use  of 
the  fitting  procedure  is  eliminated. 

However,  a  capacity  sort  was  used  in  the  tie  breaking  sort. 

Once  all  items  of  a  set  have  one  of  three  linearly  dependent  character¬ 
istics  in  common,  it  does  not  matter  which  of  the  other  two  characteristics 
is  used  to  break  ties  as  the  schedule  will  be  identical.  The  capacity 
sort  was  used  because  it  is  equivalent  to  and  easier  to  program  than  a 
size  sort  with  jobs  of  equal  duration. 
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2.4  Examples  of  Scheduling  Rules  (Example.  Problem  2) 

The  four  primary  rules,  RANDOM  and  RANDOMFIT  are  applied  to  a 
roster  of  five  jobs  scheduled  through  a  shop  with  fixed  capacity  of 
four  units.  The  roster  is  defined  by  capacity  and  duration  vectors  as 
follows: 

x  =  2,  3,  1,  2,  4 
P  =  2,  4,  3,  1,  2 

Figure  5  shows  the  schedules  produced  and  the  order  in  which  an  attempt 
is  made  to  scheule  the  jobs. 
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The  numbers  in  the  diagrams  indicate  the  position  the  job  had 
in  the  original  roster. 

From  this  example,,  a  tendency  for  rules  based  on  a  smallest 
first  priority  to  pile  up  jobs  at  the  start  can.  be  seen.  A  combination 
of  the  fitting  routine  and  a  largest  first  discipline  tends  to  intermix 
large  and  small  jobs  resulting  in  a  more  compact  arrangement.  The 
measures  of  effectiveness  are  calculated  for  each  of  these  schedules  in 


Section  3.3o 
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CHAPTER  3 


EVALUATION  OF  SCHEDULING  RULES 


3  •  1  Data  Used  to  Evaluate  the  Rules 

The  data  used  consists  of  sixty  rosters  containing  twenty-five 
jobs  each  and  a  shop  with  nine  units  of  capacity  running  as  long  as 
necessary  to  process  the  roster.  For  a  given  roster,  the  same  process 
was  used  to  generate  both  the  capacity  and  duration  vector. 

The  rosters  were  generated  with  two  different  processes  and 
three  different  means,  the  result  being  the  set  of  sixty  rosters  is 


broken  into  six  groups  of  ten  rosters  each  as  follows: 

1.  Uniform  distribution  range  1  to  5 

2.  Uniform  distribution  range  I  to  7 

3.  Uniform  distribution  range  1  to  9 


Symbol  F3 
Symbol  F4 


Symbol  F5 


4.  Truncated  Poisson  distribution  mean  3  Symbol  P3 

5.  Truncated  Poisson  distribution  mean  4  Symbol  P4 

6.  Truncated  Poisson  distribution  mean  5  Symbol  P5 

The  uniformly  distributed  rosters  were  generated  using  a 

library  function  (in  the  computer)  which  generates  randomly  distributed 
numbers  within  a  specified  range.  To  generate  the  Poisson  rosters,  the 
random  distribution  from  this  library  function  was  used  in  a  formula 
which  transforms  a  random  distribution  into  a  Poisson  distribution. 

A  chi  square  test  is  performed  on  each  of  the  six  sets  of 
rosters  to  indicate  to  what  degree  the  generated  distribution  approaches 
the  intended  distribution.  The  results  of  the  test  are  calculated  in 


Tables  6  and  7  and  are  summarized  in  Table  5  below. 


, 
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TABLE  5 

Chi  Square  Test  Results  on  Rosters 


Roster 

Type 

(  X2  Capacities 

X2  Durations 

X2  Max.* 

Uniform 

1  to 

5 

.56 

5.32 

9.488 

Uniform 

1  to 

7 

9.85 

11.00 

12.592 

Uniform 

1  to 

9 

10.71 

12.29 

15.507 

Poisson 

Mean 

3 

7.22 

8.89 

11.070 

Poisson 

Mean 

4 

2.11 

3.52 

14.067 

Poisson 

Mean 

5 

4.99 

4.86 

15.507 

*  ( .05  level 

of 

s ignif icance) 

Table  5  shows  that  the  differences  between  the  observed  and 
theoretical  distributions  are  in  no  case  too  large  to  be  attributed  to 
chance,  and  it  is  a  reasonable  assumption  that  they  are  generated  by  the 
intended  distributions. 

The  Poisson  distribution  is  modified  by  truncations  at  both 
ends.  Values  of  zero  are  rejected  since  these  represent  a  null  job,  and 
values  greater  than  nine  are  rejected  since  these  represent  jobs  of 
greater  capacity  than  the  shop  can  handle.  Both  truncations  represent 
the  real  life  situation  as  a  job  must  have  a  greater-than-zero  size,  and 
a  shop  must  occasionally  reject  a  job  that  is  beyond  its  capabilities. 

Zero  elements  were  eliminated  by  adding  1  to  all  elements 
preserving  the  shape  of  the  distribution  but  shifting  it  horizontally  or 
increasing  the  mean  by  1.  Elements  of  greater  value  than  nine  were 
eliminated  by  either  regenerating  the  distribution  when  a  larger  element 
occurred  or  replacing  the  offending  element  by  a  random  element  from  one 


to  nine . 


■ 


■ 
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TABLE  6 

Chi  Square  Calculation  -  Uniform  Rosters 


Observation 


Capacities 


Durat ions 


Frequency 

2 

Expected  (F-E)  Frequency 

Expected 

(F-E) 

Range  1 

to  5 

E 

E 

1 

54 

50 

16 

63 

50 

169 

2 

49 

50 

1 

53 

50 

9 

3 

47 

50 

9 

46 

50 

16 

4 

51 

50 

1 

44 

50 

36 

5 

49 

50 

1 

44 

50 

36 

X2 

=  28/50  =  .56 

X2  = 

266/50 

X2  max . 

4  d.f . 

.05  level 

of  significance  =  9.488 

Range  1 

to  7 

1 

40 

35.6 

19.4 

34 

35.6 

2.6 

2 

46 

35.6 

108.2 

42 

35.6 

43.6 

3 

31 

35.6 

21.2 

34 

35.6 

2.6 

4 

32 

35.6 

13.0 

50 

35.6 

207  .4 

5 

37 

35.6 

2.0 

31 

35.6 

21.2 

6 

23 

35.6 

158.8 

25 

35.6 

112.4 

7 

41 

35.6 

29 .2 

34 

35.6 

2.6 

• 

X2  = 

351.8/35.6  =  9.85 

X2  =  392.4/356 

X  max . 

6  d.f 

.  .05  level 

of  significance  =  12.592 

Range  1 

to  9 

1 

28 

28 

0 

29 

28 

1 

2 

40 

28 

144 

20 

28 

64 

3 

21 

28 

49 

31 

28 

9 

4 

28 

28 

0 

35 

28 

9 

5 

31 

28 

9 

40 

28 

144 

6 

26 

28 

4 

20 

28 

64 

7 

31 

28 

9 

34 

28 

36 

8 

19 

28 

81 

24 

28 

16 

9 

26 

28 

4 

27 

28 

1 

X2 

=  300/28  =  10.71 

X2  - 

344/28  ^ 

X2  max . 

8  d.f 

.  .05  level 

of  significance  =  15.507 

=  5.32 


-  11.0 


=  12.29 


. 
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TABLE  7 

Chi  Square  Calculations  -  Poisson  Rosters 


Observation 

Capacities 

Durations 

Frequency 

Expected 

(F-E)2 

Frequency 

Expected 

(F-E)2 

Mean  3 

E 

E 

1 

36 

34 

.12 

43 

34 

2.40 

2 

56 

69 

2.45 

66 

69 

.13 

3 

74 

69 

.36 

68 

69 

.01 

4 

45 

46 

.02 

51 

46 

.54 

5 

30 

22 

3.04 

11 

22 

5.50 

6 

7 

9 

1.23 

9 

9 

.31 

7 

1 

3 

2 

3 

8 

1 

1 

0 

1 

X2 

=  7.22 

X2 

=  8.89 

X^  max .  5  d . 

f.  .05  level 

of  significance  = 

11.070 

Mean  4 

1 

12 

12 

12 

10 

12 

.33 

2 

42 

37 

.68 

39 

37 

.11 

3 

59 

57 

.07 

59 

57 

.07 

4 

49 

57 

1.12 

53 

57 

.16 

5 

42 

42 

37 

42 

.60 

6 

26 

25 

.04 

30 

25 

1.00 

7 

13 

12 

.08 

15 

12 

.75 

8 

7 

5 

.12 

3 

5 

.50 

9 

0 

2 

2 

2 

9 

1 

1 

1 

X2  =  2.11  X2  =  3.52 

X2  max.  7  d.f.  .05  level  of  significance  =  14.067 


Mean  5 


1 

6 

5 

.20 

3 

5 

.80 

2 

17 

18 

.06 

18 

18 

3 

29 

37 

1.73 

36 

37 

.03 

4 

55 

49 

.74 

47 

49 

.08 

5 

52 

49 

.18 

44 

49 

.51 

6 

39 

39 

48 

39 

2.07 

7 

24 

26 

.15 

30 

26 

.62 

8 

12 

15 

.60 

15 

15 

9 

16 

7 

1.33 

9 

7 

.75 

9 

X2  max .  8  d.f. 

.05  level 

5 

X2  =  4.99 
of  significance  = 

15.507 

5 

X2 

=  4.86 
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TABLE  9 

Comparison  of  Poisson  to  Uniform  Distribution 


Observation 

Capacities 

Durat ions 

Frequency 

Expected 

(F-E)2 

Frequency 

Expected  (F-E)2 

E 

E 

1 

5 

28 

484 

3 

28 

625 

2 

17 

28 

121 

18 

28 

100 

3 

29 

29 

1 

36 

28 

64 

4 

55 

28 

729 

47 

28 

361 

5 

52 

28 

576 

44 

28 

256 

6 

39 

28 

121 

48 

28 

400 

7 

24 

28 

16 

30 

28 

4 

8 

12 

28 

256 

15 

28 

169 

9 

16 

28 

144 

9 

28 

361 

X2  =  2448/28  = 

87.5 

x2= 

2340  =  83.5 

28 

Similarly  the 

uniformly 

generated 

rosters  of 

range  1 

to  9  are 

compared  to  a 

Poisson  di 

str ibut ion 

of  mean 

5: 

1 

28 

5 

29 

5 

2 

40 

18 

87.1 

20 

18 

29.2 

3 

21 

37 

6.9 

31 

37 

1.0 

4 

28 

49 

9.0 

25 

49 

11.7 

5 

31 

49 

6 . 6 

40 

49 

1.7 

6 

26 

39 

4.3 

20 

39 

9.3 

7 

31 

26 

1.0 

34 

26 

2.5 

8 

19 

15 

1.1 

24 

15 

5.4 

9 

26 

7 

16.3 

27 

7 

18.7 

9 

0 

5 

0 

5 

X2  =  132.3  X2  =  79.5 


The  generated  Poisson  distributions  are  also  tested  for  goodness 
to  fit  to  a  theoretical  uniform  distribution  with  the  same  mean  and  vice 
versa.  Note  that  a  uniform  distribution  with  range  1  to  9  has  a  mean  of 
5  and  is,  therefore,  compared  to  the  Poisson  mean  5  distribution.  The 
results  of  the  comparison  are  summarized  in  Table  8  below  and  the  cal¬ 


culations  shown  in  Table  9  above. 
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TABLE  8 


Comparison  to  Opposite  Distribution 


Comparison 


2  2  2 
X  Capacit ies  X  Pur at  ions  X  *  Max. 


Generated  Poisson  to  Uniform 


87.5 


83.5 


20.09 


Generated  Uniform  to  Poisson 


132.3 


79.5 


20.09 


*  .01  level  of  significance 

The  differences  between  actual  and  presumed  distributions  are 
too  great  to  be  attributable  to  chance;  and,  therefore,  we  conclude  that 
these  hypothesized  fits  are  not  good. 

It  is,  therefore,  reasonable  to  state  that  the  data  generated 
approximates  the  intended  distribution  for  each  group.  Each  roster  consists 
of  twenty-five  jobs  drawn  at  random  from  these  groups  simulating  the 
situation  where  the  job  characteristics  follow  some  known  distribution, 
and  the  shop  receives  a  random  sample  from  it . 

It  should  also  be  noted  that  the  truncation  and  shifting  does 
not  affect  the  Poisson  distribution  unduly.  The  truncations  actually 
occurred  rarely  as  for  a  Poisson  distribution  with  a  mean  of  4,  the 
theorectical  distribution  states  that  one  element  in  250  will  exceed  9. 
Therefore,  there  is  a  .1  chance  of  one  truncation  occurring  in  each 
roster.  For  a  mean  of  5,  there  is  a  .5  chance  of  one  truncation  in  any 
roster  assuming  only  the  durations  will  be  truncated. 


3 . 2  Cr iter ia  to  Measure  Effectiveness  of  Scheduling  Rules 
3.2.1  Definition  of  Measures  of  Effectiveness 

Nine  measures  of  effectiveness  were  developed  to 
evaluate  the  fourteen  scheduling  rules  described  in  Section  1.5.1. 
Some  of  these  measures  have  been  outlined  in  Section  1.5.2,  and  a 
more  complete  definition  of  them  follows. 
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3 . 2 . 1 . a  Mean  Start  Time 

This  is  the  measure  of  effectiveness  A  described 
in  Section  1.5.2  or  (+/NUM)  v^NUM. 


3 . 2 . 1 .  b 

Time  Shop  Occupied 

This  is  the  measure  of  effectiveness  B  of 

Section 

1.5.2  or  T/F. 

3 . 2 . 1 .  c 

Average  Number  of  Jobs  in  Process 

The  sum  of  the  number  of  jobs  in  process  each 

time  period  divided  by  the  number  of  time  periods  the  shop  is 
occupied  or  (+/JOB)  7  ([*/F). 

3 . 2 . 1 .  d  Overall  Efficiency 

The  number  of  units  of  work  required  by  the 
roster  divided  by  the  number  of  units  of  capacity  occupied  by 
it  or  (+/ (PxX) )  7  (CC  x  f/F) . 

3 . 2 . 1 .  e  Mid-Range  Efficiency 

This  is  the  measure  of  effectiveness  C  and  the 
same  as  3.2.1.d,  except  the  second  two  quarters  of  the  period 
of  time  the  shop  is  occupied  are  considered  instead  of  the  whole 


period . 

The  period  considered  is: 

1  - 

■  SLACK  7  CC  x  2  x  BEGIN 

where 

BEGIN  =  ( (+/F)  -  (4/  f/F))  7  4 

and 

SLACK  z  +/C  [[BEGIN  +  {  2  x  BEGIN^) 

3 .2 .1 .  f 

Disutility  Factor 

This  is  the  measure  of  effectiveness  D  and  is 

defined  as  +/PxXx  (NUM-1) . 

3 . 2 . 1 . g  Sum  of  Start  Times 


Or  +/NUM. 
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3.2.1.h  Standard  Deviation  of  Start  Times 


This  is  defined  as  ((+/(NUM-VECTMS)*2) 


NUM)*.5 


where  VECTMS  =  (p  NUM) 


~  o  NUM) . 


3 . 2 . 1 .  i  Square  Root  of  Variance  of  Start  Times  About  Zero 


This  is  similar  to  the  standard  deviation  except 


the  focal  point  is  zero  and  not  the  mean  of  the  distribution. 
It  is: 

((+/NUM*2)  f  ^  NUM)- .5 . 

The  quantities  in  the  above  definitions  are  as  defined 
in  Section  1.1  except  CC  is  C  CO.  Three  functions  are  defined  in 
terms  of  these  quantities  and  may  be  interpreted  by  reference  to  an 
APL  manual .  Cl]  They  are  defined  as  follows: 

BEGIN  -  the  lower  quartile  point  of  the  shop  occupied  period 
taken  to  the  nearest  integer. 

SLACK  -  the  number  of  unused  units  of  shop  capacity  between 
the  first  and  third  quart iles. 

VECTMS  -  a  vector  of  the  same  order  as  the  roster  vectors 


with  elements  equal  to  the  mean  start  time. 


In  this  notation,  the  operations  are  performed  according 


to  a  right -to-left  rule  except  where  otherwise  indicated  by  parentheses. 
The  basic  symbols  mean  the  following: 

+/  -  summation  over  the  entire  vector 
T/  -  maximum  element  in  vector 
^  -  number  of  elements  in  the  vector 
A*r  -  A  raised  to  the  power  r 

The  more  complex  notation  is  explained  in  Section  4.1.3  and  APL 


manuals . 
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3.2.2  Elimination  o £  Measures  of  Effectiveness 

Of  the  nine  measures  described  in  Section  3.2.1,  only 
the  four  labelled  A,  B,  C,  and  D  were  used  in  the  evaluation  of  the 
scheduling  rules.  Two  were  eliminated  because  they  were  not  relevant 
to  the  problem  at  hand,  and  three  were  eliminated  because  they  were 
linearly  related  to  others  which  were  retained. 

Measures  3.2.1.h  and  3.2.1.j  were  intended  to  show  how 
the  start  times  are  dispersed  and  point  to  schedules  where  one  job 
is  unduly  delayed  while  the  majority  are  processed  quickly.  In  the 
extension  of  this  study  to  a  dynamic  case  as  described  in  Section  1.3 
where  there  is  a  flow  of  jobs  into  the  roster  as  well  as  out,  those 
measures  would  be  relevant.  They  would  indicate  the  scheduling 
techniques  which  allow  a  job  to  be  kept  away  from  the  shop  by  the 
incoming  stream  of  higher  priority  jobs.  The  measures  are,  therefore, 
included  to  give  the  model  greater  flexibility  in  spite  of  their 
lack  of  usefulness  to  the  problem  at  hand. 

Measures  3.2.1a  and  3.2.1.g  are  trivially  related.  They 
are  (+/NUM)  ~  ^  NUM  and  +/NUM  where  NUM  is  the  order  of  the  roster 
vector,  a  constant.  The  mean  start  time  was  chosen  as  being  more 
illustrative  of  the  two . 

A  linear  relationship  also  exists  between  the  following 

measures: 

3 . 2 . 1  .b  =  B  =  f/F 

3.2.1. C  =  B'  =  (+/J0B)  i  (r/F) 

3.2 .1 . d  =  B"  =  (+/(PxX) )  -  (CC  x  P/F) 

To  compare  B  and  B"  note  that  in  B"  P,  X  and  CC  are  all  constant 
given  a  roster  and  shop  and  are,  therefore,  independent  of  the 
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scheduling  rule.  C /F  is  the  only  factor  dependent  on  the  scheduling 
rule  and,  therefore: 

B"  =  k  =  k  where  k  =  +/ (PxX) 

P/F  B  CC 

This  is  a  linear  inverse  relationship  between  measures  3.2.1.b  and 
3.2.1.d.  The  constant  is  known. 

To  show  that  B  and  B'  are  related,  we  must  show  that 
for  a  given  roster  and  shop  +/JOB  is  independent  of  the  scheduling 
rule . 

Each  job  runs  for  a  number  of  units  of  time  equal  to  its 
duration  requirement  and  will  occupy  this  time  sooner  or  later  in 
the  schedule.  No  matter  when  it  occupies  the  shop  or  how  many 
other  jobs  are  running,  each  job  will  eventually  increment 
JOB  £n,  N  +  1,...,M]  by  1  where  N  is  its  start  time  and  M  is  its 
finish  time.  Therefore,  +/JOB  is  the  sum  of  the  durations  of  the 
jobs  in  the  roster.  This  is  a  function  of  the  roster  and  not  the 
schedule  and,  therefore: 

B '  =  k  =  k  where  k  =  +/JOB 
T/F  B 

The  number  of  units  of  time  the  shop  is  occupied  was  chosen  as 
representative  of  these  measures. 

In  summary,  the  four  measures  of  effectivenss  chosen  as 
indicated  in  Section  1.5.2  are: 

A  -  3. 2.1. a  Mean  start  time 

B  -  3.2.1.b  Time  shop  occupied 

C  -  3.2.1.e  Mid-range  efficiency 


D  -  3 .2 . 1 .  f 


Weighted  disutility  factor 


' 
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3 • 3  Application  of  the  Measures  of  Effect iveness 

For  each  of  the  sixty  rosters,  the  nine  measures  of  effective¬ 
ness  were  calculated  and  displayed  by  computer  for  each  of  the  fourteen 
scheduling  rules.  As  an  example  of  this  procedure,  Example  Problem  2 
of  Section  2.4  is  considered.  The  four  measures  of  effectiveness  used 
are  calculated  for  each  of  the  six  schedules  in  Table  10. 

TABLE  10 

Measures  of  Effectiveness  -  Problem  2 


Common  Quantities 

+/ (PxX)  =4+  12+3+2+8=  29 
CC  =  C  [l]  =  4 


RANDOM  Rule 

NUM  =  1,  3,  3,  7,  8 
F  =  2,  6,  5,  7,  9 
BEGIN  =  2 
SLACK  =  1 

A  =  (1  +  3  +  3  +  7  +  8)  7  5  =  4.4 

B  =  Max.  [2,  6,  5,  7,  9]  =9 

C  =  1  -  (1/(4  x  2  x  2))  =  .94 

D  =  (4  x  0)  +  (12  x  2)  +  (3  x  2)  +  (2  x  6)  +  (8  x  7)  =  98 

RANDOMFIT  Rule. 

NUM  =  1,  3,  1,  7,  8 
F  =  2,  6,  3,  7,  9 
BEGIN  =*  2 


SLACK  =  3 


- 
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RANDOMFIT  Rule  (Cont'd.) 

A=(l+3+l+7+8)  -5=4.0 
B  =  Max.  [2,  6,  3,  7,  9]  =  9 
C  =  1  -  ( 3 / ( 4  x  2  x  2))  =  .81 

D  =  (4  x  0)  +  (12  x  2)  +  (3  x  0)  +  (2  x  6)  +  (8  x  7)  =  92 

SS  Rule 

NUM  =2,  6,  1,  1,  4 
F  =  3,  9,  3,  1,  5 
BEGIN  =  2 
SLACK  =  2 

A  =  (2  +  6  +  1  +  1  +  4)  r  5  =  2.8 
B  =  Max.  [3,  9,  3,  1,  5]  =  9 
C  =  1  -  (2/(4  x  2  x  2))  =  .88 

D  =  (4  x  1)  +  (12  x  5)  =  (3  x  0)  +  (2  x  0)  +  (8  x  3)  =  98 

LS  Rule 

NUM  =  7 ,  1,  1,  7,  5 
F  =  8,  4,  3,  7,  6 
BEGIN  =  2 
SLACK  =  1 

A  =  (7  +1  +  1+  7  +  5)  t5  =  4.2. 

B  =  Max.  £8,  4,  3,  7,  6]]  =  8 
C  =  1  -  (V(4  x  2  x  2))  =  .94 

D  =  (4  x  6)  +  (12  x  0)  +  (3  x  0)  +  (2  x  6)  +  (8  x  4)  =  68 

SP  Rule 

NUM  =  1,  5,  5,  1,  3 
F  =  2,  8,  7,  1,  4 
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SP  Rule  (Cont'd.) 

BEGIN  =  2 
SLACK  =  0 

A=  (1+5+5+  1+3)  t  5  =  3 . 0 
B  =  Max.  Q2,  8,  7,  1,  4D  =  8 
C  =  1  -  (°/(4  x  2  x  2))  =  1.00 

D  =  (4  x  0)  +  (12  x  4)  +  (3  x  4)  +  (2  x  0)  +  (8  x  2)  =  76 

LP  Rule 

NUM  =5,  1,  1,  5,  7 
F  =  6,  4,  3,  5,  8 
BEGIN  =  2 
SLACK  =  3 

A  =  (5  +  1  +  L  +  5  +  7)  =  3.8 
B  =  Max.  C.6  +  4+  3  +  5+  8j  =  8 

C  =  1  -  (3/(4  x  2  x  2))  =  .81 

D  =  (4  x  4)  +  (12  x  0)  +  (3  x  0)  +  (2  x  4)  +  (8  x  6)  =  72 

3 . 4  Ranking  of  the  Scheduling  Rules 

3.4.1  Description  of  the  Rankings 

For  each  roster,  the  measures  of  effectivenss  were 
calculated  for  each  of  the  fourteen  different  schedules.  Due  to 
variations  in  particular  rosters,  the  absolute  values  so  calculated 
would  vary  considerably  between  rosters.  As  the  object  is  to  derive 
a  scheduling  rule  and  test  its  relative  effectiveness  regardless  of 
the  roster,  a  system  of  ranking  the  rules  relative  to  each  other 
was  developed. 

To  illustrate  this  point,  consider  two  rosters,  A  and  B, 


■ 
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and  two  scheduling  rules,  1  and  2.  The  efficiencies  resulting  from 
using  the  two  rules  on  the  rosters  are  shown  in  Table  11. 

TABLE  11 


Scheduling  Rule _ Roster 


A 

B 

1 

.9 

.95 

2 

00 

00 

• 

.92 

We  may  compare  the  performance  of  Rule  2  on  roster  B 
to  rule  1  on  roster  A;  observe  that  .92  is  better  than  .9,  and 
conclude  scheduling  rule  1  is  better. 

However,  note  that  rule  2  performs  better  on  the  same 
roster,  and  the  following  method  of  ranking  the  rules  illustrates 
this  point.  Assign  a  value  0  to  the  best  rule  for  a  given  roster 
and  a  value  1  to  the  worst  rule.  Prorate  numbers  between  0  and  1 
for  the  other  rules  depending  on  their  relative  position  between 
the  extremes  as  measured  by  the  measure  of  effectiveness.  Table 
12  shows  the  application  of  this  method  to  the  absolute  value 
rankings  in  Table  11 . 


TABLE  12 


Scheduling  Rule 

1 

2 


Roster 

A 

0 

1 


B 

0 

1 


This  shows  the  consistent  better  performance  of  rule  1 


relative  to  rule  2  and  suppresses  the  difference  in  rosters.  It 


.  n 
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points  out  clearly  the  fact  that  we  will  get  better  results  on  any 
roster  by  using  rule  1  rather  than  rule  2. 

This  ranking  procedure  was  used  to  rank  the  calculated 
values  of  the  measures  of  effectiveness  obtained  by  application  of 
the  fourteen  scheduling  rules  for  each  roster.  The  best  performance 
was  identified  by  a  zero,  and  the  worst  performance  by  a  10,000.  The 
intermediate  values  of  the  measures  of  effectiveness  were  assigned 
a  value  between  0  and  10,000  which  was  a  linear  function  of  the 
difference  between  them  and  the  best  and  worst  values. 

3.4.2  Example  of  Ranking  Procedure 

The  procedure  described  above  is  illustrated  using  the 
calculated  values  of  the  measures  of  effectiveness  from  the  example 
of  Table  9.  Table  13  summarizes  these  results: 

TABLE  13 

Summary  of  Performance  -  Example  Problem  2 


Scheduling  Rule _ Measure  of  Effectiveness 


A 

B 

C 

D 

Random 

4.4* 

9* 

.94 

98* 

Randomf it 

4.0** 

9* 

.81* 

92 

Smallest  Size 

2 . 8** 

9* 

.88 

98* 

Largest  Size 

4.2 

g** 

.94 

68** 

Smallest  Duration 

3.0 

8** 

1.00** 

76 

Largest  Duration 

3.8 

8** 

V 

CO 

* 

72 

*  Worst 

**  Best 

The  table  is  considered  column  by  column  in  order  to 
compare  all  scheduling  rules  according  to  each  measure  of 


■ 
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effectiveness.  The  best  performance  in  each  column  is  indicated 
by  **  and  the  worst  by  *. 

Note  that  the  best  value  for  measures  A,  B,  and  D  is  a 
low  number  while  for  measure  C,  it  is  a  high  number.  Consistency 
is  obtained  by  assigning  the  best  value  0  no  matter  whether  it  is 
high  or  low. 

Table  13-A  shows  the  rankings  as  they  are  calculated 
using  this  method.  To  illustrate  the  method,  the  rank  for  the 
element  marked  (C<)  on  Table  12  is  calculated: 

10000  x  4.0  -  2.8  =  7500 
4.4  -  2.8 


TABLE  13-A 


Summary  Rankings 

-  Example 

Problem 

2 

Scheduling  Rule 

Measure  of  Effectiveness 

A 

B 

C 

D 

Random 

10000 

10000 

3333 

10000 

Randomf it 

7500°^ 

10000 

10000 

8000 

Smallest  Size 

0 

10000 

10000 

10000 

Largest  Size 

8750 

0 

3333 

0 

Smallest  Duration 

1250 

0 

0 

2667 

Largest  Duration 

s 

6250 

0 

10000 

1333 

3.4.3  Application  of  the  Ranking  to  the 

Scheduling  Rules 

Using  the  data  from  sixty  rosters,  a  tabulation  similar 
to  Table  13  was  made.  The  fourteen  scheduling  rules  were  evaluated 
according  to  each  of  the  four  measures  of  effectiveness  and  an 
average  over  the  sixty  rosters  computed.  The  standard  deviation  was 
also  calculated  and  used  as  explained  in  Section  3.5. 


m 
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TABLE  14 


MEAN  AND 

STANDARD  DEVIATION  OF 

RANKINGS 

1  Mean  Start 

Time  -  Measure  A 

Rule 

Number 

Mean 

St.d.  Dev. 

RANDOM 

60 

8731 . 25 

1583.83 

RANDOMFIT 

60 

3313.43 

863.42 

SS 

60 

527.32 

604.03 

LS 

60 

8825.82 

1382.28 

SP 

60 

311  .  5 

340 . 82 

LP 

60 

6986 . 75 

1409 . 24 

SPSC 

60 

231 . 25 

318.89 

SPLC 

60 

443.45 

439.89 

LPSC 

60 

6603  .  23 

1408 . 99 

LPLC 

60 

7651.68 

1379.98 

SSSP 

60 

434 . 63 

575.04 

SSLP 

60 

636 . 35 

578.44 

LSSP 

60 

8847.75 

1366.67 

LSLP 

60 

8679.9 

1489 . 56 

2.  Shop  Occupied  Time  - 

Measure  B 

Rule 

Number 

Mean 

Std.  Dev. 

RANDOM 

60 

9347.43 

1250 . 51 

RANDOMFIT 

60 

3554.4 

2029.84 

SS 

60 

6599.68 

2385 . 07 

LS 

60 

342 . 67 

727 . 45 

LS 

60 

4621 .48 

2550.19 

LP 

60 

2244 .93 

2062.68 

SPSC 

60 

4758.28 

2317 . 07 

SPLC 

60 

3962.85 

2475 . 43 

LPSC 

60 

3276 . 83 

2395 . 51 

LPLC 

60 

1411.72 

1788 . 45 

SSSP 

60 

6136 . 7 

2256.92 

SSLP 

60 

5684.9 

3306 

LSSP 

60 

242 . 32 

624.19 

LSLP 

60 

496 . 52 

857.62 

' 

1 
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The  average  relative  performance  of  the  fourteen  scheduling 
rules  as  measured  by  each  measure  of  effectiveness  presented  in 
Table  14.  From  these,  simple  rules  of  thumb  based  on  job  character¬ 
istics  for  scheduling  jobs  through  a  job  shop  will  be  obtained. 

The  standard  deviations  will  be  used  to  test  the  significance  of 
hypothesized  relationships. 

In  addition  to  Table  14,  Table  15  shows  the  position 
ranking  of  each  scheduling  rule  without  regard  for  the  varying 
differences  between  rules  as  shown  in  Table  14.  This  is  useful  in 
visualizing  initial  relationships  and  suggesting  hypotheses  for 
future  testing. 


TABLE  15 

Rankings  of  Scheduling  Rules 


Scheduling  Rule 

Measure 

of 

Effectiveness 

A 

B 

C 

D 

Random 

12 

14 

14 

14 

Randomf it 

7 

7 

8 

6 

SS 

5 

13 

12 

12 

LS 

13 

2 

2 

2 

SP 

2 

9 

9 

9 

LP 

9 

5 

5 

5 

SPSC 

1 

10 

10 

10 

SPEC 

4 

8 

7 

7 

LPSC 

8 

6 

6 

8 

LPLC 

10 

4 

4 

4 

SSSP 

3 

12 

11 

11 

SSLP 

6 

11 

13 

13 

LSSP 

14 

1 

1 

1 

LSLP 

3.4.4  Conclusions 

11 

From  the  Rankings 

3 

3 

3 

Examining  the  data  in  Tables  14  and  15,  several  patterns 
become  evident.  These  will  be  stated  now  and  tested  for  significance 


in  Section  3.5. 


•fr 
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3.4.4a  No  matter  what  measure  of  effectiveness  is  used  to 
evaluate  it,  the  f irst -come-first-served  rule  is  worst  or 
nearly  worst.  Note  that  it  is  always  worse  than  random 
fitting . 


3 .4 .4  .b 
markedly 
based  on 
durat ion) 
quantity . 


Two  different  groups  of  scheduling  rules  perform 
different  from  each  other.  One  of  the  groups  is 
priority  given  to  a  small  quantity  (size  or 
and  the  second  is  based  on  priority  to  a  large 


3.4.4. C  Random  fitting  assignment  ranks  between  the  two 
groups  according  to  all  measures  of  effectiveness. 

3.4.4. d  The  smallest  first  priority  group  of  scheduling 
rules  gives  better  mean  start  time  while  the  largest  first 
group  gives  better  shop  occupied  time,  efficiency,  and  weighted 
disutility  factor. 

The  idea  that  fitting  or  searching  ahead  in  the  roster 
for  a  job  which  will  utilize  the  unused  capacity  when  the  highest 
priority  job  exceeds  it  is  good  is  confirmed  by  conclusion  3.4.4.e. 

If  this  relationship  can  be  shown  to  be  statistically  significant, 
we  are  justified  in  omitting  the  capacity  based  priority  rules  as 
discussed  in  Section  2.3.  The  two  rules  RANDOM  and  RANDOMFIT  were 
introduced  to  illustrate  this  point  as  the  only  difference  between 
them  is  one  uses  the  fitting  procedure  while  the  other  does  not. 

We  have  chosen  rules  which  behave  in  different  ways 


and  measures  of  effectiveness  such  that  rules  rank  differently 


according  to  different  measures.  The  fact  that  random  fitting 
falls  between  the  two  groups  of  rules  according  to  each  measure  of 
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effectiveness  shows  that  we  have  chosen  scheduling  rules  with 
opposite  characteristics.  We,  therefore,  appear  to  have  a  fairly 
comprehensive  set  of  simple  scheduling  rules. 

Table  15  is  rearranged  in  order  to  illustrate  relation¬ 
ships  which  occur  as  a  function  of  the  second  or  tie  breaking  priority 
assignment .  Table  16  shows  the  same  data  as  Table  15  except  the 
two  secondary  scheduling  rules  of  the  same  primary  type  are  grouped 
around  the  primary  scheduling  rule. 

TABLE  16 

Rankings  of  Scheduling  Rules 


Scheduling  Rule 

Measure 

of 

Effectiveness 

A 

B 

C 

D 

SSSP 

3 

12 

11 

11 

SS 

5 

13 

12 

12 

SSLP 

6 

11 

13 

13 

SPSC 

1 

10 

10 

10 

SP 

O 

z. 

9 

9 

9 

SPLC 

4 

8 

7 

7 

LSSP 

14 

1 

1 

1 

LS 

13 

2 

2 

2 

LSLP 

11 

3 

3 

3 

LPSC 

8 

6 

6 

8 

LP 

o 

5 

5 

5 

LPLC 

10 

4 

4 

4 

Except  in  one  case, 

the  behavior 

of 

the  second  sort 

i; 

consistent.  This  is  shown  by  each  group  of  three  numbers  monotonically 
increasing  or  decreasing. 

The  conclusions  regarding  the  secondary  or  tie  breaking 
sort  are  presented  in  3.4.4.e,  3.4.4. f,  3.4.4.g,  and  3.4.4.h. 
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3.4.4. e  If  the  primary  scheduling  rule  is  smallest  size 
first,  priority  to  smallest  duration  improves  all  measures 
A,  B,  C,  and  D. 

3.4.4. f  If  the  primary  scheduling  rule  is  smallest  duration 
first,  priority  to  smallest  capacity  improves  measure  A  and 
priority  to  largest  capacity  improves  measures  B,  C,  and  D. 

3.4.4. g  If  the  primary  scheduling  rule  is  largest  size  first, 
priority  to  largest  duration  improves  measure  A  and  priority 
to  smallest  duration  improves  measures  B,  C,  and  D. 

3.4.4. h  If  the  primary  scheduling  rule  is  largest  duration 
first,  priority  to  smallest  capacity  improves  measure  A  and 
priority  to  largest  capacity  improves  measures  B,  C,  and  B. 

A  second  sort  priority  which  behaves  in  a  certain  manner 
according  to  measure  A  appears  to  behave  in  the  opposite  manner 
according  to  measures  B,  C,  and  D. 

Table  17  shows  only  the  primary  scheduling  rules  and 
their  rankings  as  taken  from  Table  15  while  Table  18  is  a  condensation 
of  it  using  the  ranking  numbers  1  to  4  only.  From  these  tables,  the 
performance  of  the  four  primary  groups  of  scheduling  rules  can  be 
easily  compared. 


TABLE  17 

Primary  Scheduling  Rule  Rankings 

Scheduling  Rule  ^  Measure  of  Effectiveness 

A_ _ B_ _ _ C _ D  _ 

5  13  12  12 

13  2  2  2 

2  9  9  9 

9  5  5  5 


SS 

LS 

SP 

LP 
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TABLE  18 

Primary  Scheduling  Rule  Rankings 


Scheduling  Rule 


SS 

LS 

SP 

LP 


Measure  of  Effectiveness 


A  B_ C 

2  4  4 

4  11 

13  3 

3  2  2 


D 

4 

1 

3 

2 


The  main  conclusions  of  the  study  are  drawn  from  these 


tables .  They  are  3.4.4. j  and  3.4.4.k. 


3.4.4. j  To  obtain  the  least  mean  start  time,  give  priority 
to  the  job  with  smallest  duration.  Ties  are  broken  with 
priority  to  smallest  capacity  according  to  3.4.4. f. 

3.4.4. k  To  obtain  the  least  time  to  complete  all  jobs,  high¬ 
est  efficiency  or  minimum  disutility  factor,  give  priority 

to  the  job  with  largest  size.  Ties  are  broken  with  priority 
to  smallest  duration  according  to  3.4.4.g. 

The  relative  value  of  other  priority  schemes  is  given  by  the  rankings 
in  Table  15 . 


There  are  apparent  conflicts  which  may  be  statistically 
significant  between  the  performance  of  a  priority  rule  in  the  primary 
and  secondary  sort.  As  an  example,  priority  to  smallest  duration 
results  in  better  mean  start  times  on  the  primary  sort  while  on  the 
tie  breaking  sort  of  jobs  in  size  order  the  reverse  is  true. 
Explanation  of  significant  discrepancies  is  made  in  Section  3.5.3 
following  the  statistical  tests. 
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3 . 5  Comparison  of  the  Scheduling;  Rules 

3.5.1  t-Test  for  Significance  of  Difference  of  Means 

A  t-test  for  the  difference  of  means  is  used  to  decide 
whether  the  observed  differences  in  the  scheduling  rules  is  sig- 
nificant  or  due  to  chance.  [6]  Basic  to  this  method  is  the  assumption 
that  we  have  drawn  our  sample  from  a  normally  distributed  population 
of  rankings. 

If  the  entire  population  of  rosters  were  formed  in  the 
same  way,  that  is,  had  a  set  of  jobs  which  would  interlock  into  the 
same  pattern,  a  unique  rank  as  indicated  by  each  measure  of  effective¬ 
ness  would  be  observed  (although,  different  for  each  measure  of 
effectiveness)  for  each  scheduling  rule  regardless  of  the  roster. 

The  ranks  are  different  for  each  roster  because  there  are  many 
chance  variables  affecting  the  makeup  of  each  roster  which  in  turn 
affects  the  interlocking  pattern. 

By  the  Central  Limit  Theorem,  the  cumulative  effect  of 
a  multitude  of  random  variables  approaches  normal  distribution.  In 
this  case,  the  cumulative  result  of  the  variables  affecting  the  inter¬ 
locking  pattern  is  the  calculated  value  of  the  measure  of  effectiveness 
It  is,  therefore,  reasonable  to  conclude  that  there  is  a  mean  rank 
for  each  scheduling  rule,  and  actual  ranks  for  each  individual 
roster  are  normally  distributed  about  it  with  a  fixed  standard 
deviation.  Our  sample  consists  of  sixty  such  actual  ranks. 

The  t-test  for  difference  between  means*  is  a  technique 
for  deciding,  given  a  sample  from  two  populations,  whether  there  is 


*  See  reference  CbU  p .  233. 
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a  difference  in  the  population  means.  In  this  problem,  the 
population  means  are  the  mean  ranks  assigned  to  the  scheduling  rules 
by  a  measure  of  effectiveness;  and  these  are  to  be  shown  to  be 
different  by  inference  from  the  sample  of  sixty  rosters.  When  this 
is  done,  we  have  established  that  one  rule  is  better  than  another 
when  we  are  trying  to  optimize  a  given  measure  of  effectiveness. 
Accordingly,  t-tests  described  in  the  above  reference  were  performed 
on  all  possible  combinations  of  two  items  from  each  column  of 
Table  14. 

3.5.2  Significance  Tests  on  Hypothesized  Conclusions 

The  relationships  that  were  hypothesized  from  Tables  14 
to  18  and  presented  in  Section  3.4.4  are  tested  for  statistical 
significance.  All  tests  are  performed  at  the  0.05  level  of  sig¬ 
nificance  unless  otherwise  stated,  and  the  number  of  degrees  of 
freedom  is  118.  (Sixty  items  in  each  sample  np  +  n2  -  2  =  60  +  60  - 
2  =  118.) 

3. 5. 2. a  To  test  that  the  random  first -come-first -served 
rule  ranks  last  according  to  measures  of  effectiveness 
B,  C,  and  D,  its  rank  is  compared  to  the  rule  which  ranks 
second  last  according  to  each  measure  (see  Section  3. 4. 4. a). 
Table  19  shows  that  there  is  a  significant  difference  at  the 
0.01  level  (t>  2.617)  between  this  rule  and  the  next  worst 
rule  according  to  all  three  measures  of  effectiveness. 
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TABLE  19 

Comparison  of  RANDOM  to  Next  Worst  Rule 


Measure 

RANDOM  Rank 

Next  Rule 

Next  Rank 

t -Value 

B 

9347 

SS 

6599 

7.83 

C 

8507 

SSLP 

7212 

3.02 

D 

9735 

SSLP 

5982 

10.78 

The  random  first- 

come-first -served 

rule  is  not 

the  worst 

rule  according  to 

measure  of  effectiveness  A. 

However , 

Table  20  shows  that  the  two  rules  which  are  worse  than  it 
are  not  significantly  worse. 


TABLE  20 

Comparison  of  RANDOM  to  Worst  Rules  -  Measure  A 


RANDOM  Rank 

Worse  Rule 

Worse  Rule  Rank 

t -Value 

8731 

LS 

8825 

0.35 

8731 

LSSP 

8847 

0.43 

The  fact  that  the  random  first -come-first -served  rule  is 
the  worst  according  to  nearly  all  measures  of  effectiveness 
has  therefore  been  established. 

It  has  also  been  shown  that  the  fitting  procedure  is 
necessary  regardless  of  the  measure  of  effectiveness  in 
Table  21  by  comparing  RANDOM  and  RANDOMFIT  and  showing  them 
to  be  different  well  beyond  the  0.001  level  of  significance 
(t  >  3.373)  . 
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TABLE  21 

Comparison  of  RANDOM  and  RANDOMFIT  Rules 


of  Effectiveness 

RANDOM  Rank 

RANDOMFIT  Rank 

t -Value 

A 

8731 

3313 

23.07 

B 

9347 

3554 

18.66 

C 

8507 

3963 

10.75 

D 

9735 

2309 

36.01 

3.5.2.b  The  primary  scheduling  rules  are  significantly 
divided  into  two  groups,  one  containing  LP  and  LS ,  the  other 
containing  SP  and  SS  ,  (See  Section  3.4.4.b.)  To  show  that 
the  rules  are  significantly  different  in  each  group,  the  rule 
is  chosen  which  ranks  closest  to  the  other  group  according  to 
each  measure  of  effectiveness.  The  ranks  of  the  two  rules 
thus  chosen  are  compared  in  Table  22  which  shows  there  is  a 
difference  between  them  at  the  0.001  level  of  significance 
(t >  3.373)  . 


TABLE  22 

Test _o f  Group  Separation  -  Primary  Rules  Only 

Shortest  Largest 

Measure  of  Effectiveness  First  Group  First  Group  t -Value 


Rule 

Rank 

Rule 

Rank 

A 

SS 

527 

LP 

6986 

32.36 

B 

SP 

4621 

LP 

2244 

5.56 

C 

SP 

4742 

LP 

2442 

5.12 

D 

SP 

3090 

LP 

1605 

4.70 

If  the  secondary 

rules 

based 

on  these  primary  rules 

are 

included  in  their 

respect ive 

groups,  an 

overlapping 

between 

■ 
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the  groups  occurs  according  to  measures  C  and  D.  For  measure 
A,  there  is  still  a  strong  separation  while  the  separation 
according  to  measure  B  may  still  be  due  to  chance.  This  is 
shown  in  Table  23. 


TABLE  23 


Test  of  Group  Separation  -  All  Rules  Included 


Measure  of  Effectiveness 


Shortest 
First  Group 

Rule  Rank 


Longest 
First  Group 

Rule  Rank 


t -Value 


A 

B 


SSSP  636 
SPLC  3962 


LPSC  6603 


LPSC  3276 


30.09 

1.52 


OVERLAP 


D 


OVERLAP 


The  primary  rules  are,  therefore,  shown  to  form  into  two 
groups  regardless  of  the  measure  of  effectiveness.  The 
groups  are  on  the  basis  of  opposite  performance  according 
to  any  measure  of  effectiveness.  By  utilizing  tie  breaking 
priority  rules,  a  degree  of  overlap  can  be  created  between 
the  groups  according  to  at  least  some  measures  of 
effectiveness . 

3.5.2 cC  To  show  that  the  RANDOMFIT  rule  lies  between  the 
two  groups  of  scheduling  rules  according  to  all  measures  of 
effectiveness,  its  rank  is  compared  to  that  of  the  primary 
scheduling  rule  closest  to  each  side  of  the  intergroup 
boundary.  (See  Section  3.4.4.C.)  These  are  the  same  rules 
that  were  compared  in  Table  22.  As  indicated  in  Table  24, 


•« 


the  RANDOM  rule  is  significantly  different  from  each  group 
at  the  0.05  level  except  for  a  marginal  case  indicated  by 

(*). 
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TABLE  24 

Location  of  the  RANDOMFIT  Rule  -  Primary  Only 

Measure  of  RANDOMFIT 

Effectiveness  Rank  Lower  Boundry  Upper  Boundry 


Rule 

Rank 

t-Diff . 

Rule 

Rank 

t -Value 

A 

3313 

SS 

527 

20.30 

LP 

6986 

17.07 

B 

3554 

LP 

2244 

3.47 

SP 

4621 

2.51 

C 

3963 

LP 

2442 

3.43 

SP 

4742 

1.65* 

D 

2309 

LP 

1605 

2.73 

SP 

3090 

2.49 

*  Significant  at  the  0.10  level 

If  the  secondary  rules  are  included  in  the  groups,  there  is 
no  significant  separation  for  RANDOMFIT  to  fall  into  according 
to  measures  B,  C,  and  D.  However,  according  to  measure  A,  it 
is  separated  from  the  lower  and  upper  boundaries  with  greater 
than  99.9  percent  confidence  (t>  3.373)  (t  =  19.78  and  15.29 
for  upper  and  lower,  respectively). 

RANDOMFIT  has,  therefore,  been  shown  to  produce  results 
midway  between  those  produced  by  the  two  major  groups  in  the 
set  of  primary  scheduling  rules  according  to  all  measures  of 
effect iveness . 

3.5.2.d  This  conclusion  is  elaborated  in  conclusions  3.4.4.e 
to  3. 4. 4. k  and,  therefore,  will  be  tested  in  Sections  3.5.2.e 


to  3 .5 .2 .k . 
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3.5.2.e  The  three  scheduling  rules  in  the  SS  (smallest  size) 
priority  group  were  examined  to  find  the  effect  of  the 
.  secondary  sort-  (See  Section  3.4.4.e.)  From  Table  25  which 
shows  the  best  and  worst  tie  breaking  priority  according  to 
each  measure  of  effectiveness,  we  see  that  SP  (shortest 
duration)  is  best  according  to  all  measures  of  effectiveness 
except  B. 

TABLE  25 

Tie  Breaking  Effect  -  SS  Priority 

Measure  of  Effectiveness  Best  Tie  Breaker  Worst  Tie  Breaker  t -Value 

Rule  Rank  Rule  Rank 


A 

SSSP 

434 

SSLP 

636 

1.90 

B 

SSL? 

5684 

SSSP 

6136 

0.87 

C 

SSSP 

6630 

SSLP 

7212 

1.20 

D 

SSSP 

5262 

SSLP 

5982 

1.55 

The  conclusion  that  SP  priority  tie  breaking  produces 
better  results  according  to  measures  of  effectiveness  A,  C, 
and  D  is  illustrated  by  a  weak  relationship  in  that  direction. 
The  opposite  effect  is  noted  according  to  measure  B  although 
it  is  very  weak. 

3.5.2.f  Table  26  compares  the  best  and  worst  rule  from  the 
set  of  SP  (smallest  duration)  based  rules  according  to  each 
measure  of  effectiveness.  (See  Section  3.4.4. f.) 


. 


■< 


60 


TABLE  26 


Tie 

Breaking 

Effect  -  SP 

Priority 

Measure  of  Effectiveness 

Best 

Tie  Breaker 

Worst  Tie  Breaker  t. 

-Value 

Rule 

Rank 

Rule 

Rank 

A 

SPSC 

231 

SPLC 

443 

3.00 

B 

SPLC 

3962 

SPSC 

4758 

1.80* 

C 

SPLC 

3748 

SPSC 

5530 

3.77 

D 

SPLC 

2339 

SPSC 

3809  ” 

4.17 

Given  that 

the 

primary  priority 

is  SP,  this  shows  that 

there 

is  a  significant  tendency  at  the  0.01  level  (t>  2.617)  for 
a  smallest  capacity  priority  to  improve  measure  A  and  a  largest 
capacity  priority  to  improve  measures  C  and  D.  There  is  also 
a  strong  indication  that  the  latter  is  true  for  measure  B 
as  well  (significant  at  0.1  level). 

3.5.2.g  Table  27  below  shows  that  the  evidence  supporting 
this  conclusion  is  relatively  weak  except  for  measure  of 
effectiveness  B.  (See  Section  3.4.4.g.) 

The  rosters  in  the  LS  (largest  size)  priority  groups 
are  compared  in  the  same  way  as  before. 

TABLE  27 

Tie  Breaking  Effect  -  LS  Priority 

Measure  of  Effectiveness  Best  Tie  Breaker  Worst  Tie  Breaker  t -Value 

Rule  Rank  Rule _ Rank _ 


A 

LSLP 

8679 

LSSP 

8847 

.63 

B 

LSSP 

242 

LSLP 

496 

1.84 

C. 

LSSP 

488 

LSLP 

702 

.95 

D 

LSSP 

258 

LSLP 

367 

.80 
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3.5.2.h  As  shown  in  Table  28  which  considers  the  rules  in 
the  LP  (longest  duration)  priority  group,  the  conclusion 
of  3.4.4.h  is  verified  very  strongly  past  the  0.001  level  of 
significance  (t  >  3.373). 

TABLE  28 

Tie  Breaking  Effect  -  LP  Priority 

Measure  of  Effectiveness  Best  Tie  Breaker  Worst  Tie  Breaker  t -Value 

Rule  Rank  Rule  Rank 


A 

LPSC 

6603 

LPLC 

7651 

4.08 

B 

LPLC 

1411 

LPSC 

3276 

4.79 

C 

LPLC 

1545 

LPSC 

3520 

5.27 

D 

LPLC 

1021 

LPLC 

2431 

4.60 

3 . 5 . 2 .  j  To  show  that  the  SPSC  rule  gives  the  best  result 
according  to  measure  of  effectiveness  A,  we  compare  SP  to 
each  of  the  other  three  primary  rules  to  confirm  that  it  is 
significantly  better  than  all  of  them.  (See  Section  3.4.4.j.) 
The  secondary  rule  is  verified  by  conclusion  3.5.2.f  which 
states  that  if  the  SP  rule  is  used,  measure  A  is  improved  by 
giving  priority  to  smallest  capacity  (SC).  Table  29  shows 
that  SP  is  the  best  primary  rule  at  the  0.05  level  of 
significance  (t>1.98). 

TABLE  29 

Comparison  of  SP  to  Primary  Rules  -  Measure  A 

SP  Rank  Other  Rule  t -Value 

Rule  Rank  _  _ 


SS 

527 

2.39 

LS 

8825 

45.93 

LP 

6986 

35.36 

311 
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3 . 5 . 2 . k  To  show  that  the  LSSP  rule  gives  the  best  result 
according  to  measures  of  effectiveness  B,  C ,  and  D ,  we 
compare  LS  to  each  of  the  other  three  primary  rules  to 
confirm  that  it  is  better  than  all  of  them.  (See  Section  3.4.4.k.) 
Conclusion  3.5.2.g  states  that  if  the  LS  rule  is  used  on  the 
primary  priority,  there  is  a  strong  indication  that  priority 
should  be  given  to  smallest,  duration  to  give  the  best  result 
according  to  measures  of  effectiveness  B,  C,  and  D. 

Table  30  compares  LS  to  each  of  the  other  primary  scheduling 
rules  for  each  measure  of  effectiveness  and  shows  that  this 
rule  is  superior  to  the  other  primary  rules  at  the  0.001 
level  of  significance  (t>  3.373). 

TABLE  30 


Measure  of  Effectiveness 

LS  Rank 

Other 

Rule 

t -Value 

Rule 

Rank 

B 

342 

SS 

6599 

19.27 

SP 

4621 

12.39 

LP 

2244 

6 . 68 

C 

700 

SS 

6833 

15.06 

SP 

4742 

10.61 

LP 

2442 

5.07 

D 

354 

SS 

5660 

15.18 

SP 

3090 

9.76 

LP 

1605 

5.79 

< 


:  ' 
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Therefore,  a  one  best  primary  rule  has  been  isolated 
and  found  to  be  significantly  different  from  all  the  others 

according  to  each  measure  of  effectiveness.  A  secondary 
rule  (SPSC)  which  produces  a  unique  schedule  has  been  shown 
to  produce  a  significantly  best  result  according  to  measure 
A.  There  is  strong  indication  that  the  unique  secondary 
rule  LSSP  produces  the  best  result  according  to  measures  B, 

C,  and  D.  While  the  results  of  Table  27  which  analyzes  the 
second  sort  effect  when  the  primary  priority  is  largest  size 
(LS)  are  not  significant,  they  are  consistent  with  each  other 
and  with  other  results  which  have  been  verified  more  strongly. 

Table  18  of  Section  3.4.4  shows  the  relative  rankings 
of  each  primary  scheduling  rule  according  to  each  measure  of 
effectiveness.  To  show  that  there  is  a  significant  difference 
between  all  primary  rules,  Table  31  shows  the  t-value  of  the 
difference  between  each  rule  and  the  next  worst  rule. 

TABLE  31 

Significance  Test  on  Table  18 

Scheduling  Rule  Measure  of  Effectiveness 


A  B  C  D 


ss 

32.36 

WORST  RULE 

WORST  RULE 

WORST  RULE 

LS 

WORST  RULE 

6.68 

5.07 

5.79 

SP 

2.39* 

4.35 

4.19 

6.14 

LP 

7.16 

5.56 

5.12 

4.70 

Therefore,  the  separation  of  the  primary  scheduling 
rules  is  very  strong  as  it  is  significant  at  the  0.001 


. 
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level  (t  ^3.373)  except  for  one  case  (*)  which  is  significant 
at  the  0.01  level  of  significance. 

3.5.3  Analysis  of  General  Results 

It  seems  reasonable  that  giving  priority  to  the  smallest 
jobs  should  minimize  mean  start  time.  However,  there  is  a  significant 
difference  in  the  result  depending  on  whether  ’'smallest."  means 
smallest  size  or  smallest  duration.  On  the  tie  breaking  sort,  it  is 
noted  that  either  smallest  or  largest  duration  first  gives  the  best 
result  according  to  measure  A  depending  on  whether  the  primary  rule 
is  smallest  or  largest  size  first,  respectively .  The  reason  for  this 
shift  is  not  immediately  evident.  When  the  primary  rule  is  based  cn 
duration  rather  than  size,  a  tie  breaking  priority  to  the  smallest 
size  or  capacity  job  gives  the  best  result  according  to  measure  A 
regardless  of  whether  the  primary  rule  is  largest  or  smallest 
duration  first. 

Measures  B,  C,  and  D  all  give  the  same  general  rankings 
to  the  scheduling  rules.  It  is  reasonable  that  an  action  which 
improves  the  time  the  shop  is  occupied  will  also  improve  efficiency, 
but  it  is  not  as  evident  why  it  should  improve  the  weighted  disutility 
or  profit-loss  factor. 

For  measure  D,  the  weighted  disutility  factor,  it  is 
advantageous  to  complete  as  many  large  jobs  as  soon  as  possible. 

To  do  this,  one  of  the  rules  in  the  "largest  first"  group  is  used 
and  is  shown  to  be  superior  to  any  in  the  "smallest  first  group. 

In  the  tie  breaking  sort  when  the  primary  rule  is  largest  duration, 
priority  to  largest  capacity  further  improves  measure  D.  However, 
when  the  primary  rule  is  largest  size,  a  seemingly  opposite  effect 


‘ 
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occurs  as  priority  to  smallest  duration  improves  measure  D.  Note 
that  largest  duration  first  is  better  than  smallest  duration  first 
in  the  primary  set  of  rules  . 

For  measures  B  and  C,  shop  occupied  period  and  mid-range 
efficiency,  a  largest  duration  or  size  priority  is  also  best.  It 
is  not  immediately  evident  why  this  should  be  so  as  the  same  result 
should  occur  whether  the  schedule  or  its  mirror  image  is  considered. 

In  addition,  the  same  conflict  between  primary  and  secondary  scheduling 
rules  exists  as  does  for  measure  D. 


In  summary,  three  points  must  be  explained  or 


rationalized: 


1.  Why  do  the  "largest  first"  rules  optimize  measures  B  and 
C,  shop  occupied  time  and  mid-range  efficiency. 

2.  Why  is  measure  A  optimized  by  a  sort  according  to  duration 
while  the  other  three  measures  B,  C,  and  D  are  optimized 
by  a  sort  according  to  size. 

3.  Is  there  a  real  conflict  between  the  behavior  of  the  primary 
and  secondary  rules  and  what  is  the  reason  for  it. 

To  facilitate  the  analysis,  a  full  display  of  the 
schedule,  number  of  jobs  in  process,  and  amount  of  capacity  utilized 
in  each  time  period  as  well  as  the  summary  results  was  made  for 
five  rosters,  one  from  each  group.  An  example  of  this  display  is 
Table  32  made  for  a  roster  of  ten  jobs. 
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TABLE  32 

SAMPLE  COMPUTER  OUTPUT 
APPLICATION  OF  LSSP  RULE  TO  ROSTER  OF  TEN  JOBS 


CAPACITIES ,  PROCESS  TIMES  AND  NUMBER  OF  JOBS  IN  ROSTER 
8462739257 
5825924951 
10 

REARRANGED  ROSTER  IN  ORDER  SCHEDULING  IS  ATTEMPTED 


CAPACITIES  7 

8  9 

4  5 

2  6  2 

7  3 

PROCESS  TIMES  9 

5  4 

8  5 

9  2  5 

1  2 

SCHEDULE  IS  AS 

FOLLOWS 

PROCESS  TIME 

CAPACITY 

START 

FINISH 

9 

1 

1 

9 

9 

2 

1 

9 

5 

8 

10 

14 

4 

.9 

15 

18 

8 

4 

19 

26 

5 

5 

19 

23 

5 

2 

24 

28 

2 

3 

24 

25 

2 

6 

27 

28 

1 

7 

2-9 

29 

JOBS  IN  PROCESS  22222222211 
11111112  2.  22233 

2  2  2  1 

UNUSED  CAPACITY  00000000011 
11100000000000 
3  112 


SUMMARY  CHARACTERISTICS 

16.9  29  1.724137931  0.9540229885  0.9603174603 

2766  169  9.565040512  19.4190628 

SUMMARY  CHARACTERISTICS  ARE: 

Mean  Start  Time 
Shop  Occupied  Time 

Average  Number  of  Jobs  in  Process 

Total  Efficiency 

Mid-Range  Efficiency 

Disutility  Factor 

Sum  of  Start  Times 

Standard  Deviation  of  Start  Times 

Square  Root  of  Variance  About  O  Start  Time 
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To  explain  why  a  largest  duration  or  size  first  type  of 
priority  rule  is  best  according  to  measures  of  effectiveness  B  and  C, 
we  examine  the 
cycle . 


number 


of  jobs  running  in  each  time  period  in  the 
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FIGURE  5 
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FIGURE  6 
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FIGURE  7 
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Figures  5  to  8,  inclusive,  show  the  number  of  jobs  in 
process  in  each  time  period  for  the  four  primary  scheduling  rules 
using  a  typical  roster  from  each  of  the  groups  of  roster  types. 

The  difference  between  the  plots  for  smallest  duration 
or  size  first  rules  and  largest  duration  or  size  first  rules  is  that 
the  former  have  a  pronounced  peak  at  the  start  while  the  latter  are 
relatively  uniform  throughout  the  time  period.  This  means  that  a 
"shortest  first"  rule  tends  to  bunch  up  a  large  number  of  jobs  at 
the  start  while  the  other  type  of  rule  runs  the  same  number  of  jobs 
at  all  times. 

The  jobs  which  a  shortest  first  routine  schedule  at  the 
start  are  all  the  small  jobs  leaving  the  large  jobs  to  be  run  at  the 
end  of  the  cycle.  The  largest  first  routine  will  give  priority  to 
a  large  job  when  space  is  available  rather  than  fill  it  with  several 
small  jobs.  Often  more  than  one  large  job  cannot  be  run  at  the 
same  time  because  of  the  capacity  restriction;  and  when  a  smallest 
first  scheduling  rule  is  used,  the  large  jobs  are  trailed  out  one 
at  a  time  after  the  small  jobs  are  scheduled.  The  largest  first 
routine  will  schedule  a  small  job  in  the  unused  capacity  while  a 
large  job  is  running;  and,  therefore,  the  amount  of  shop  time  which 
would  have  been  used  scheduling  small  jobs  only  is  saved. 

The  mid-range  efficiency,  measure  C,  is  improved  by 
priority  to  large  jobs  as  these  rules  permit  the  unused  capacity 
which  would  occur  if  large  jobs  only  are  scheduled  to  be  filled  with 
small  jobs.  Therefore,  the  reason  a  largest  first  routine  improves 
both  measures  B  and  C  is  because  it  is  able  to  fit  the  jobs  together 
better  by  utilizing  small  jobs  to  fill  the  gaps  between  large  ones. 
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Measure  of  effectiveness  A  has  been  shown  to  be  optimized 
by  giving  priority  to  smallest  duration  rather  than  smallest  size. 

If  the  primary  rule  was  priority  to  smallest  size,  a  longer  than 
necessary  duration  is  introduced  at  the  start  of  the  schedule  as 
some  small  size  jobs  may  have  long  durations.  The  result  is  an 
advancement  of  the  start  and  finish  times  of  all  jobs  subsequently 
scheduled.  This  means  that  the  mean  start  time  will  be  longer  than 
if  the  smallest  duration  jobs  were  scheduled  first. 

Measures  of  effectiveness  B,  C,  and  D  require  that  the 
largest  size  rather  than  largest  duration  job  be  given  priority  on 
the  primary  sort.  As  we  have  seen,  a  good  fit  is  essential  to  good 
performance  according  to  time  shop  is  occupied  and  mid-range  efficiency 
(measure  B  and  C) .  It  must  be  then  that  a  size  rather  than  duration 
based  routine  produces  a  better  fit. 

For  a  good  fit,  there  must  be  a  great  variety  of  jobs 
available  at  each  decision  point  or  point  at  which  a  new  job  may  be 
scheduled.  Since  both  capacities  and  durations  are  randomly  generated 
and  size  is  the  product  of  the  two,  there  will  be  a  greater  variety 
of  capacities  in  an  equal  size  group  than  in  an  equal  duration  group. 
This  is  b  ecause  there  are  two  random  variables  affecting  the  size 
group  while  only  one  affects  the  duration  group. 

As  an  example  of  why  a  good  mix  is  important  to  good 
fitting,  consider  what  may  happen  when  a  three-unit  capacity  bloc 
becomes  available  in  the  first  part  of  the  cycle  and  there  are  three 
jobs,  one  requiring  two  units  and  two  requiring  one  unit  of  capacity. 

If  the  two-unit  job  occurs  first,  it  is  scheduled;  and  the  one-unit 
jobs  are  easily  scheduled  later.  If  both  one-unit  jobs  occur  ahead 
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of  the  two-unit  job,  they  are  both  scheduled  leaving  a  more  difficult 
job  to  schedule  later.  If  they  occur  mixed,  it  is  possible  to 
schedule  the  two-unit  job  immediately  as  well  as  one  one-unit  job. 

The  largest  size  first  scheduling  rule,  therefore, 
results  in  a  more  varied  roster  making  a  good  fitting  pattern  more 
likely  than  if  one  parameter  were  monotone  increasing  or  decreasing 
throughout . 


According  to  measures  of  effectiveness  B,  C,  and  D, 
largest  duration  first  is  a  better  rule  than  shortest  duration  first 
on  the  primary  sort  while  to  tie  break  it  is  better  to  choose  smallest 
duration  first  after  a  primary  sort  according  to  size.  This  is 
equivalent  to  giving  priority  to  the  largest  capacity  job  on  the 
secondary  sort.  By  doing  this,  if  a  large  pocket  of  unused  capacity 
occurs  at  the  start  of  the  schedule,  it  is  filled  immediately  by  a 
large  capacity  job.  If  it  was  filled  by  a  small  capacity  job,  this 
would  leave  the  other  larger  job  until  later  where  only  small 
capacity  vacancies  may  occur.  The  fit  would  be  poorer  as  these 
vacancies  would  be  left  unoccupied  because  the  small  jobs  were 
scheduled,  and  the  larger  job  would  run  by  itself  at  the  end. 


3.6 


Unusual  Rosters 

3.6.1  C lassif icat ion  of  Unusual  Rosters 

An  unusual  roster  is  one  which  consistently  produces 
unusual  results  according  to  any  or  all  of  the  measures  of  effective¬ 
ness.  The  ranks  of  each  scheduling  rule  according  to  each  measure 
of  effectiveness  for  the  sixty  rosters  are  considered  and  the  mean 
and  standard  deviation  calculated.  From  this,  a  95  percent 
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confidence  interval  can  be  calculated;  and  an  unusual  result  is  a 
rank  that  falls  outside  this  confidence  interval. 

Table  33  shows  the  frequency  of  occurrence  of  unusual 
ranks  for  each  roster  according  to  each  measure  of  effectiveness, 
and  Table  34  is  a  composite  for  all  measures  of  effectiveness.  Both 
horizontal  and  vertical  scales  in  the  tables  are  used  to  identify 
the  rosters  where  the  horizontal  scale  identifies  the  roster  type 
as  follows: 


10 

-  Uniform 

roster 

range 

20 

-  Uniform 

roster 

range 

30 

-  Uniform 

roster 

range 

40 

-  Poisson 

roster 

mean 

50 

-  Poisson 

roster 

mean 

60 

-  Poisson 

roster 

mean 

The  numbers  on  the  vertical  scale  represent  the  individual  roster 
in  each  of  the  groups  above.  For  example,  the  code  number  for  the 
roster  represented  by  column  40  and  row  5  is  45  and  is  a  Poisson 
roster  with  mean  3.  As  it  is  in  the  sixth  row,  it  is  the  sixth  in 
the  matrices  containing  Poisson  rosters  with  mean  3.  By  referring 
to  Section  4.3.3,  the  location  of  the  capacities  for  this  roster  is 
R0SCP3  [6;]  and  the  durations  are  located  in  the  matrix  ROSDP3[6j. 
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TABLE  33 

Frequency  of  Unusual  Rosters 


Average  Start _ Time 


10 

20 

30 

40 

50 

60 

2 

2 

5 

5 

2 

1 

1 

2 

2 

1 

2 

1 

2 

4 

1 

2 

1 

6 

Shop  Occupied 

Time 

1 

1 

1 

5 

1 

1 

1 

1 

2 

1 

2 

1 

2 

1 

2 

2 

2 

2 

1 

4 

2 

2 

1 

1 

Average  Efficiency 

1 

3 

1 

5 

2 

2 

1 

2 

1 

5 

2 

1 

1 

1 

1 

1 

3 

1 

1 

2 

1 

3 

1 

Disutility  Factor 

1 

1 

3 

1 

1 

1 

1 

2 

2 

3 

2  . 

1 

2 

2 

2 

1 
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7 
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TABLE  34 

Composite  of  Table  33 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

Total 


10 _  20  30 


4  4  2 

5  1  13 

5  3  4 

5  6 

3  4 

1  2 

7  2 

2  8 

1  6 

23  30  35 


40  50  60 


1  1  6 

5  1 

3  3  4 

12  8 

2 

3  1 


4 

2  15 

3  3  6 

22  13  40 


The  maximum  element  in  any  cell  in  Table  33  is  14  which 
would  occur  if  every  scheduling  rule  applied  to  the  roster  gave  an 
unusual  result  according  to  the  measure  of  effectiveness.  If  the 
technique  has  isolated  5  percent  of  the  cases  as  deviant  (corresponding 
to  95  percent  confidence  interval),  there  should  be  5  percent  of 
14  x  4  x  60  or  168  observations  in  Table  34.  There  are  actually 
163  observations. 

Table  35  shows  the  number  of  rosters  with  each  fequency 
of  unusual  results.  Note  that  the  average  roster  would  produce  2.7 
unusual  results  according  to  any  of  the  four  measures  of  effectiveness. 
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TABLE  35 

Number  of  Rosters  and  Frequency  of  Unusual  Results 


Frequency 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Number 

19 

8 

6 

8 

6 

4 

4 

1 

2 


1 

1 


60  Rosters 

There  are  five  rosters  whose  performance  could  be  con¬ 
sidered  unusual;  of  these,  three  are  moderately  unusual,  and  two 
are  extremely  unusual.  The  matrix  location  codes  (see  Section  4.3.3, 
Table  47)  for  the  moderately  and  extremely  unusual  rosters  are  shown 
in  Table  36 . 


TABLE  36 


Location  of  Unusual  Rosters 


Code 

Moderate  27 

38 

63 

Extreme  31 

68 


Matrix  Location 
R0SF4  [8;] 
R0SF5  [9;'| 
ROSE 5  [4;] 
R0SF5  [2;] 
R0SP5  [9;] 


Description 
Uniform  Range  1-7 
Uniform  Range  1-9 
Poisson  Mean  5 
Uniform  Range  1-9 
Poisson  Mean  5 
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The  last  two  rosters  were  analyzed  to  try  and  find  a  simple  reason 
for  their  unusual  performance. 

Another  roster  which  is  of  interest  is  code  60  or 
R0SP5  DO  •  This  roster  has  the  greatest  number  of  unusual 
occurrences  of  any  roster  according  to  measure  A  and  none  according 
to  any  other  measure  of  effectiveness. 

3.6.2  Analysis  of  Unusual  Rosters 

Referring  to  Table  36,  it  may  appear  that  a  roster  with 
large  average  size  relative  to  the  size  of  the  shop  is  prone  to 
unusual  results.  However,  Table  27  shows  that  seven  out  of  twenty 
or  35  percent  of  the  largest  rosters  (code  30  and  60)  have  no  unusual 


o 


ccurrences  at  all,  while  nineteen  in  sixty  or  32.5  percent  of  the 


whole  sample  have  no  unusual  occurrences.  It,  therefore,  seems  as 
likely  that  a  large  roster  wTill  have  no  unusual  occurrences  as  an 
average  sized  one. 

The  column  totals  of  Table  34  show  the  frequency  of 
unusual  occurrences  for  each  roster  class.  By  using  a  chi  square 
test  for  goodness  of  fit,  we  compare  these  to  a  uniform  distribution 
of  unusual  occurrences  to  see  if  there  is  any  significant  relation¬ 
ship  between  roster  size  or  type  and  the  frequency  of  unusual 
occurrences.  A  chi  square  value  is  calculated  as  X  =  14.3;  and, 
therefore,  the  hypothesis  that  a  uniform  distribution  applies  to 
the  data  is  rejected  at  the  0.05  level  of  significance  (X5  max  =  11.07) 
However,  when  the  two  extremely  unusual  rosters  are  deleted  from 
consideration,  X~  =  8.2;  and  the  fit  is  acceptable  at  the  5  percent 
level  of  significance. 

Therefore,  the  proneness  to  unusual  results  does  not  seem 
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to  be  associated  with  size  of  the  jobs  in  the  roster.  If  a  simple 
causal  factor  is  to  be  found,  it  must  be  related  to  the  distribution 
of  jobs  in  the  roster;  e.g.,  dispersion  of  jobs,  the  way  capacities 
and  durations  are  related,  etc. 

The  unusual  uniformly  distributed  roster  is  considered 
first.  It  has  been  drawn  out  of  a  population  which  has  been  shown 
in  Table  6,  Section  3.1,  to  follow  a  uniform  distribution  with 
range  1  to  9 .  The  unusual  roster  is  compared  in  the  same,  manner  to 

O 

a  theoretical  uniform  distribution,  and  X  values  of  9.90  and  8.40 
are  calculated  for  capacity  and  duration,  respectively.  At  the 
5  percent  level  of  significances  (X-  max  =  15.507),  we  cannot  reject 
the  hypothesis  that  this  roster  forms  a  uniform  distribution.  In 
fact,  it  is  somewhat  closer  to  the  theoretical  distribution  than 
the  population  from  which  it  is  drawn.  Therefore,  no  unusual 
distribution  of  the  roster  has  been  found. 

The  capacities  and  durations  may  associate  in  an  unusual 
manner;  i.e.,  all  large  capacities  with  large  durations,  etc. 

Table  37  shows  the  capacities,  durations,  and  sizes  of  jobs  in  this 
roster . 

TABLE  37 

Capacities,  Durations  and  Sizes  -  Roster  R0SF5  D J 

Capacity  628293312678 

Duration  637225886732 

Size  36  6  56  4  18  15  24  8  12  42  21  16 

7776163372755 
5365  4  344  55431 

35  21  42  30  4  18  12  12  35  10  28  15  5 


Capacity 

Duration 

Size 
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The  number  of  times  a  large  or  small  capacity  has 
associated  with  a  large  or  small  duration  has  been  tabulated  in 
Table  38  where  "large"  and  "small"  denote  falling  above  and  below 
the  mean,  respectively. 

TABLE  38 

Capacity  -  Duration  Association  -  Roster  R0SF5 


Duration 

Capacity 

Large 

Small 

Large 

6 

4 

Small 

8 

7 

This  shows  a  slight  tendency  to  have  a  greater  number 
of  large  capacity  and  small  duration  jobs  although  the  effect  is  not 
pronounced.  There  is  also  a  tendency  for  the  jobs  to  be  smaller 
than  normal,  as  eight  exceed  the  normal  size  of  25  while  seventeen 
jobs  are  less  than  it. 

The  standard  deviations  of  the  capacities  and  durations 
are  a  measure  of  their  dispersion.  The  capacities  of  the  roster  in 
question  have  a  standard  deviation  of  2.41  compared  to  an  average 
for  the  set  of  ten  of  2.53.  It  is  fourth  from  the  lowest.  The 
durations  have  the  lowest  standard  deviation,  1.88,  compared  to  an 
average  of  2.47.  This  is  somewhat  low,  but  it  may  be  noted  that  the 

next  lowest  duration  standard  deviation  roster  has  only  six  unusual 

rankings . 

Table  39  shows  the  scheduling  rules  for  which  this 
roster  had  unusual  rankings.  There  may  be  a  complex  mix  of  jobs  in 

^Le  roster  arranged  in  such  a  way  as  to  make  one  scheduling  rule 
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extremely  useful  or  useless  distorting  the  results  according  to  any 
one  measure  of  effectiveness.  In  Table  39,  the  first  number  indicates 
the  rank  of  the  roster  in  question;  and  the  second  is  the  average 
rank  for  all  rosters. 

TABLE  39 

Comparison  of  Unusual  Results  of  RQSF5  C20  to  Normal 


Scheduling  Rule 

Measure  of 

Effectiveness 

A 

B 

C 

D 

RANDOM 

RANDOMFIT 

SS 

LS 

t 

2143/342 

3604/700 

SP 

LP 

8571/2244 

9352/2442 

7300/1605 

SPSC 

SPLC 

LPSC 

8571/3276 

9352/3250 

7551/2431 

LPLC 

8571/1411 

9352/1545 

7139/1021 

SSSP 

SSLP 

LSSP 

LSLP 

2143/496 

3604/702 

The  fact  that  only  one  class  of  scheduling  rules  is 
subject  to  unusual  occurrences  suggests  that  the  unusual  occurrences 
are  not  due  to  random  placement  of  jobs  in  the  roster  which  would 
make  a  random  rule  very  good  or  bad  distorting  the  rest  of  the  results. 
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A  characteristic  of  this  roster  interacts  V7ith  the  LP  scheduling 
priorities  producing  much  worse  than  usual  results  according  to 
measures  of  effectiveness  B,  C,  and  D. 

These  are  the  measures  which  require  that  the  jobs  be 
fit  well  together  to  optimize.  Since  there  seems  to  be  no  unusual 
roster  characteristics  and  there  are  sufficient  small  jobs  to  make 
the  fitting  procedure  work,  we  must  conclude  that  this  unusual  result 
cannot  be  attributed  to  prominent,  simple  roster  characteristics. 

It  is  due  to  a  complex  interaction  of  jobs  which  affects  their 
ability  to  nest  into  the  shop  without  leaving  large  blocs  of  unused 
capacity . 

The  other  extreme  roster  is  a  Poisson  distributed  roster 
with  mean  5  and  is  drawn  from  a  population  which  has  been  shown 
(Table  7,  Section  3.1)  to  have  these  characteristics.  A  c.hi  square 
comparison  of  this  roster  frequencies  to  theoretical  frequencies  for 
this  distribution  is  given  in  Table  40. 


TABLE  40 

Frequency  Comparison  -  Roster  R0SP5  C9 p 


Observation 


Capacity 

Frequency  Expected  (F-E) ^ 

E 


Duration 

Frequency  Expected  (F-E) 

E 


10-5 

2  0  1.8 

3  0  3.7 

4  5  4.9 

3  8  4.9 

6  5  3.9 

7  5  2.6 

8  0  1.3 

9  2  .7 

9  0  .5 


x^  limit  for  8  degrees  at  freedom 


.50 

1 

.5 

.50 

1.80 

0 

1.8 

1.80 

3.70 

3 

3.7 

.13 

0.00 

1 

4.9 

3.10 

1.96 

12 

4.9 

10.25 

0.31 

5 

3.9 

.31 

2.22 

2 

2.6 

.14 

1.50 

1 

1.5 

.17 

.53 

0 

.7 

1.20 

0 

.5 

=  12.52 

X2  = 

17.70 

and  .05 

level  of 

significance  = 

15.507 
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This  roster  does  not  follow  the  Poisson  distribution 
as  well  as  the  average  as  is  high  and,  in  fact,  too  high  for  the 
durations.  Examination  of  the  roster  shows  it  is  more  tightly  packed 
around  the  mean  than  is  the  Poisson  distribution.  The  standard 
deviation  of  the  capacities  is  the  second  lowest  in  the  group;  and 
for  the  durations,  it  is  the  lowest.  The  chi  square  test  confirms 
that  the  tendency  of  the  durations  to  be  5  or  6  is  significantly 
different  than  that  which  would  be  expected  from  a  Poisson  distribution. 

Table  41  is  similar  to  Table  39  and  illustrates  the 
scheduling  rules  which  produce  the  unusual  ranks.  Again,  the  first 
number  indicates  the  rank  of  R0SP5  [jl ;]  and  the  second  number  is 
the  average  rank  for  all  rosters. 

TABLE  41 

Comparison  of  Unusual  Results  from  R0SP5  C9  ;1  to  Normal 
Scheduling  Rule  Measure  _o f  Effectiveness 


A  BCD 


RANDOM 

7392/9735 

RANDOMFIT 

SS 

2266/527 

LS 

5962/354 

SP 

LP 

8239/2242 

6201/1605 

SPSC 

SPLC 

0/2339 

LPSC 

•  9453/6603 

10000/3276 

10000/3250 

10000/2431 

LPLC 

SSSP 

2277/434 

5455/1411 

SSLP 

LSSP 

LSLP 

2344/636 

5835/258 

3790/367 
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From  this  table,  a  general  tendency  to  produce  poorer 
results  is  noted  for  any  scheduling  rule  type  measured  by  any  measure 
of  effectiveness.  There  is  a  strong  tendency  for  a  random  rule  to 
be  better  than  the  usual  best  sorting  and  fitting  rule;  and  for 
measure  of  effectiveness  D,  four  fitting  rules,  SS ,  LPSC,  SSSP,  and 
SSLP,  are  worse  than  random  first -come-first-served . 

For  a  typical  roster,  the  random  first-come-first-served 
rule  is  the  worst;  and  the  other  rules  which  incorporate  the  fitting 
procedure  show  a  marked  degree  of  improvement  over  it.  In  the  roster, 
the  jobs  must  be  arranged  in  such  a  way  as  to  permit  random  first- 
come-first-served  to  produce  nearly  as  good  a  fit  as  a  fitting  routine 
in  other  words,  the  jobs  are  already  in  nearly  the  order  in  which  a 
fitting  routine  would  assign  them.  Therefore,  the  other  rules  do 
not  make  much  of  an  improvement  over  first-come-first-served;  and 
their  ranks  show  considerably  worse  than  average. 

Also,  the  very  tight  compaction  in  this  roster  means  the 
jobs  are  tending  toward  being  identical.  In  the  extreme  case  where 
the  jobs  are  identical,  it  would  not  matter  which  rule  was  used  as 
the  schedules  would  all  be  the  same. 

The  roster  R0SP5  D0;U  is  of  interest  because  all 
unusual  occurrences  were  indicated  by  measure  A  and  none  were 
indicated  by  the  other  measures.  While  SP  remains  the  best 
scheduling  rule  for  this  roster  according  to  measure  A  as  is  the 
case  in  the  general  result,  the  results  with  the  LS  group  are 
improved;  and  with  the  SS  group,  they  are  worse. 

Again  in  this  roster,  there  is  no  unusual  pattein  to 
the  manner  in  which  capacities  and  durations  are  associated,  the 
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distribution  of  the  capacities,  and  durations  or  the  size  of  the 
jobs  in  the  roster.  There  appears  to  be  no  simple  explanation  for 
the  unusual  behavior  of  this  roster  as  indicated  by  measure  A. 


The  above  approach  has  been  to  isolate  an  unusual  roster 
and  try  to  find  some  characteristic  which  causes  the  unusual  behavior. 
Except  for  deciding  that  it  does  not  matter  what  rule  is  applied  to 
a  roster  of  identical  jobs,  the  effort  is  unsuccessful  and  the 
deviations,  therefore,  can  only  be  attributed  to  chance  arrangements 
of  the  jobs  in  the  roster.  A  practical  method  of  predicting  what 
roster  will  behave  unusually  and  how  it  will  behave  before  it  is 
scheduled  has  not  been  found. 

The  apposite  approach  is  now  taken,  that  is,  isolate  a 
possible  factor  which  may  cause  unusual  results  and  see  if  it  in 
fact  does.  The  most  promising  causal  factor  appears  to  be  standard 
deviation  of  the  items  in  the  roster.  The  Poisson  and  uniform 
rosters  with  highest  and  lowest  standard  deviations  are  examined  to 
see  if  they  result  in  an  inordinate  number  of  unusual  occurrences. 
Table  42  shows  the  extreme  rosters  and  the  number  of  unusual  ranks 
occurring  in  each. 
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TABLE  42 

Roster  Identification  Standard  Deviation  Unusual  Ranks 


Uniform  Rosters 

Capacity:  High 

-  F5  [1;] 

(30) 

2.85 

2 

Low 

"  F3  £8;] 

(17) 

1.28 

0 

Duration:  High 

-  F5  [7;] 

(36) 

2.72 

2 

Low 

-  F3  £8;3 

(17) 

1.26 

0 

Poisson  Rosters 

Capacity:  High 

-  p5  £3;] 

(62) 

2.21 

4 

Low 

-  P3  £4;] 

(43) 

1.11 

1 

Duration:  High 

-  P5  £2;] 

(61) 

2.09 

1 

Low 

-  P3  £6;-.] 

(45) 

1.04 

3 

From  this,  it  is  seen  that  the  rosters  with  the  most 
unusual  standard  deviations,  either  high  or  low,  have  about  the 
expected  number  (2.7)  of  unusual  occurrences.  Therefore,  there  appears 
to  be  no  relationship  between  standard  deviation  and  frequency  of 
unusual  occurrences.  It  has  also  been  shown  at  the  start  of  this 
section  that  there  is  no  relationship  between  job  size  and  frequency 
of  unusual  occurrences . 

3 . 7  Reliability  of  the  General  Rules 

If  a  rule  for  identifying  rosters  which  will  behave  unusually 
cannot  be  found,  we  are  interested  in  the  probability  that  the  rule 
which  has  statistically  been  shown  to  be  best  will  actually  be  best 
for  any  given  roster  and  a  specified  measure  of  effectiveness.  lo 
find  this  probability,  we  find  out  how  many  times  a  roster  produces 
an  unusual  result  with  the  set  of  scheduling  rules  associated  with 


« 
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the  primary  rule  which  optimizes  according  to  each  measure  of 
e f f ect ivens s .  In  other  words,  if  the  optimizing  rule  generally  is 
LS,  unusual  occurrences  with  LS,  LSSP,  and  LSLP  are  considered;  and 
unusual  results  with  the  other  twelve  rules  are  ignored. 

Table  43  shows  the  number  of  unusual  occurrences  for  each 
roster  with  only  the  set  of  three  optimizing  scheduling  rules 
considered  for  each  measure  of  effectiveness.  It  in  effect  shows 
the  number  of  times  the  best  scheduling  rule  does  not  work. 

TABLE  43 

Number  of  Unusual  Occurrences  With 
Optimizing  S cheduling  Rules _ 


Roster 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

*  Unusual  according 
**  Unusual  according 
Other  -  Unusual 


10  20  30 


2  4 

3**  i* 

5  4 


2* 

1* 

measure  A  only, 
all  measures, 
to  measures  B 


40  50 


1  1 

2 

1* 

1 

1* 

2 

1*  2 

C,  and  D  only. 


to 
to 

according 


60 


2* 

6 


It  is  interesting  to  note  that  with  the  exception  of  one  roster, 


rosters  which  are  not  well  behaved  according  to  measure  A  are  always 
well  behaved  according  to  measures  B,  C,  and  D;  and  rosters  which 
are  not  well  behaved  according  to  B,  C,  or  D  are  always  so  according 
to  measure  A. 

This  shows  that  the  scheduling  rules  are  fairly  reliable.  For 
measure  A,  only  eight  of  sixty  rosters  show  any  extreme  results;  and 
these  do  not  show  extreme  results  according  to  the  other  measures. 
For  criteria  B,  C,  and  D  combined,  only  eleven  of  sixty  rosters  show 
any  extreme  results;  and  only  four  have  more  than  three  extremes 
which  corresponds  to  one  for  each  measure  of  effectiveness.  A 
totally  ill-behaved  roster  would  have  a  total  of  twelve  in  its  cell 
in  Table  42,  nine  unstarred  and  three  starred.  Only  two  unstarred 
rosters  exceed  even  half  of  nine,  and  two  starred  numbers  exceed 
half  of  three. 

Therefore,  the  scheduling  rules  are  fairly  reliable  as  they 
show  a  tendency  not  to  work  at  all  on  only  four  out  sixty  or  7 
percent  of  the  rosters.  The  best  scheduling  rule  has  also  been 
shown  to  be  statistically  significantly  better  than  the  others. 


CHAPTER  4 


PROGRAMMING  OF  SCHEDULING  MODELS 


4.1 


4.1.1  The  APL  System 

APL/360  is  a  version  of  Iverson's  language  executable 
on  the  IBM  360/67  computer  at  the  University  of  Alberta.  LG  The 
user  has  direct  access  to  the  computer  via  a  terminal  consisting  of 
a  modified  typewriter,  a  small  box  of  circuitry,  and,  in  some  cases, 
a  dataphone.  The  computer  can,  therefore,  be  used  at  any  place 
served  by  telephone  lines. 

To  gain  access  to  the  system,  the  user  turns  the  terminal 
on  or  if  it  is  serviced  by  dataphone,  dials  the  number  of  the 
telephone  line  allotted  to  the  computer  in  the. same  way  a  standard 
telephone  call  is  made.  He  then  enters  his  identification  number 
using  the  keyboard  and  is  recognized  by  the  system  and  given 
access  to  his  private  data  and  programs  and  the  public  libraries. 

All  programs  and  data  are  stored  in  "workspaces,"  and 
each  user  is  assigned  one  or  more  of  these.  He  can  define  a  program 
and  store  it  in  his  private  workspace  or  can  copy  and  store  programs 
from  the  public  libraries.  The  programs  in  the  public  libraries 
are  also  stored  in  similar  workspaces.  Each  workspace  is  quite 
small,  however,  as  it  contains  about  30,000  "bytes"  while  two  or 
four  bytes  are  required  to  store  a  number. 

Each  user  has  a  library  of  workspaces  which  is  identified 


by  his  identification  number  (the  same  one  used  to  gain  access  to 
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the  system).  Each  workspace  within  a  user's  private  library  and 
programs  and  data  in  the  workspace  are  named  by  the  user.  The 
public  libraries  are  similarly  arranged  with  named  workspaces  and 
named  programs  within  these  workspaces. 

Once  a  user  is  signed  on  by  typing  his  number  into  the 
sj^stem,  he  frequently  wants  to  do  one  of  the  things  described 
below . 

1.  Use  a  previously  defined  workspace.  If  the  workspace  has 
previously  been  defined  and  labelled  WSNAME,  he  types  )LOAD 
WSNAME  into  the  system.  The  system  responds  by  indicating  the 
date  WSNAME  was  last  used  and  stored  and  loads  its  contents 
into  what  is  known  as  an  "active  workspace."  There  is  one  of 
these  for  each  terminal.  The  contents  of  WSNAME  are  left 
unaltered  as  a  result  of  following  action  on  the  terminal  which 
only  alters  the  contents  of  the  active  workspace. 

2.  Save  the  active  workspace.  The  contents  of  the  active 
workspace  replace  the  contents  of  WSNAME  as  stored  permanently 
in  the  computer  when  the  command  )SAVE  WSNAME  is  entered  into 
the  system.  The  user's  work  in  a  current  session  essentially 
consists  of  alterations  to  the  active  workspace.  If  the  "save" 
command  above  is  not  entered  before  the  session  is  ended,  the 
contents  of  the  active  workspace  are  erased  when  the  terminal 
is  turned  off,  and  the  work  done  in  the  session  is  lost. 

3.  Obtain  another  workspace.  Apply  to  the  computer  centre  for 
another  workspace.  At  the  next  session  type  )SAVE  ANOTHER  which 


identifies  the  new  workspace  by  the  label  ANOTHER.  The  system 
responds  by  typing  ANOTHER  SAVED  "DATE"  or  if  this  is  attempted 
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without  requesting  another  workspace,  by  SAVE  RATION  EXCEEDED. 

The  user  then  types  )LOAD  ANOTHER,  and  the  contents  of  ANOTHER 
are  loaded  into  the  active  workspace  for  use,  modifications, 
and  storage. 

4.  Copy  a  program  or  data  from  another  workspace.  If  ANOTHER 
is  currently  loaded  into  the  active  workspace  and  a  program 
PROGNAME  is  defined  and  stored  in  the  other  workspace  WSNAME, 
the  command  )COPY  WSNAME  PROGNAME  copies  the  program  into  the 
active  workspace. 

If  a  program  to  which  access  is  required  is  stored  in 
another  user's  private  workspace  and  his  number  is  NUMBER  and 
his  workspace  name  is  WORKSPACE,  the  command  )COPY  NUMBER 
WORKSPACE  PROGNAME  transfers  the  program  from  the  other  user's 
private  library  to  the  active  workspace. 

If  a  public  library  is  called  LIBRARY  and  a  workspace 
in  this  public  library  is  called  PACKAGE,  the  command  )COPY 
LIBRARY  PACKAGE  PROGNAME  is  used  to  copy  a  program  from  the 
public  library  to  the  active  workspace.  Data  is  handled  in  the 
same  manner  as  programs . 

The  procedure  necessary  to  start  up  the  system  and  prepare 
it  for  operation  is,  therefore,  very  simple.  In  spite  of  the  small 
size  of  the  individual  workspace,  it  is  possible  to  arrange  each 
workspace  so  that  only  programs  and  data  relevant  to  the  problem 
at  hand  are  actually  stored  in  it. 

Since  a  number  of  active  workspaces  are  being  processed 
at  any  one  time,  this  is  actually  a  time  sharing  system  even  though 
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each  user  has  an  illusion  of  continuous  response.  At  present,  there 
is  space  for  31  terminals  to  be  processed  by  a  no  priority  sequencing 
method . 

4.1.2  Operating  Techniques 

Several  features  are  incorporated  in  the  system  which 
make  it  possible  to  program,  debug,  and  alter  programs  and  data 
rapidly.  Direct  access  to  the  computer  makes  it  possible  to  try  a 
program  as  soon  as  it  is  written  and  make  alterations  to  it  while 
it  is  being  executed. 

Once  the  system  has  been  set  up  in  a  way  described  in 
4.1.1,  the  user  may  wish  to  do  the  following  operations: 

1.  Define  a  program.  Type  any^PROGNAME  which  defines  the 
program  to  be  entered  as  PROGNAME .  A  [l]  is  typed  by  the 
system,  and  the  user  enters  the  first  line  of  his  program 
followed  by  a  carriage  return.  A  CO  is  then  typed  and  so  on 
until  the  user  enters  the  last  line  followed  byV 

2.  Define  data.  If  the  item  to  be  entered  is  to  be  labelled 
DATA,  type  DATA  <—  followed  by  the  data  to  be  so  labelled. 

3.  Execute  a  program.  Type  the  program  name  PROGNAME  into 
the  system.  Provided  all  the  input  required  by  the  program  is 
defined,  the  program  is  executed  and  information  is  generated, 
stored,  and  displayed  as  required  in  the  program.  If  an  error 
is  encountered,  the  execution  is  suspended,  the  offending 
statement  is  displayed  with  a  diagnostic  message,  and  control 
passed  to  the  user  so  a  correction  may  be  made. 

4.  Display  data  .  Type  the  date  name  DATA  into  the  system. 

5.  Display  a  program.  Type  the  command^PROGNAME  [Q3  V  into 
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the  system.  The  program  is  displayed  as  it  was  entered  but  is 
not  executed. 

6.  Correct  a  program.  If  the  user  wishes  to  correct  statement 
N  of  program  PROGNAME,  he  enters  the  command  ^PROGNAME  CnU 
followed  by  a  carriage  return.  The  system  responds  by  typing 

DO  and  the  user  then  enters  the  corrected  statement  followed 
by  a  y  . 

7.  Resume  execution  of  a  program.  If  execution  of  a  program 
has  been  suspended,  a  signal  PROGNAME  M  is  displayed  at  the 
time  of  suspension;  and  N  is  the  statement  number  at  which 
suspension  occurred.  To  resume  execution,  type-£>N  followed 
by  a  carriage  return. 

8.  Delete  a  program  or  data.  Deletion  is  done  using  the 

comma  nd^PROGNAMEV  or  V  DATA  \!  . 

9.  Stop  execution  of  a  program.  A  program  can  be  suspended 

by  pushing  the  ATTN  button  on  the  keyboard.  This  causes  the 
statement  PROGNAME  to  be  typed  and  control  passes  to  the 

keyboard.  The  contents  of  the  workspace  can  be  stored,  and 
the  user  can  start  the  next  session  at  the  point  he  finished 
the  last  one  by  typing.  -~>  N .  If  the  program  is  stopped  and 
it  is  not  desired  to  run  it  to  completion,  typing  — ^>,0  until 
a  )A  results  in  no  output  will  clear  the  system. 

The  above  indicates  that  the  mechanics  of  the  operation 
of  the  system  permits  very  rapid  programming  due  to  instant  ei i or 
feedback  and  the  ease  with  which  corrections  are  made.  A  more 
detailed  explanation  for  someone  wanting  to  use  the  system  is 
given  in  the  reference  manual.  D'J 
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4.1.3  The  APL  Language 

The  APL  language  is  so  constructed  that  an  algorithm  can 
be  expressed  in  a  form  close  to  its  mathematical  statement.  A  few 
symbols  are  predefined  as  operators  while  any  other  combination  of 
letters  can  be  used  to  identify  variables  or  programs.  In  the 
execution  of  statements,  a  right  to  left  precedence  rule  is  employed 
except  when  overruled  by  parentheses.  That  is,  the  right-hand 
argument  for  any  operator  is  the  entire  expression  to  the.  right, 
and  the  left-hand  argument  is  the  immediate  symbol  to  the  left. 

Vectors  and  matrices  are  handled  in  the  same  way  as 
numbers  by  most  of  the  operators.  Test  loops  for  the  end  of  the 
array  are  unnecessary  greatly  simplifying  the  programs  compared  to 
what  would  be  required  in  some  other  programming  languages.  For 
example,  the  sum  of  two  matrices  is  defined  as  another  matrix 
consisting  of  the  sums  of  corresponding  elements.  The  APL  command 
to  add  quantities  A  and  B  and  store  the  result  in  C  is  C  A  3-  B. 

If  A  and  B  are  matrices,  this  automatically  specifies  C  according 
to  this  definition  of  matrix  addition;  and  if  A  and  B  are  different 
sizes,  an  error  is  indicated. 

The  operators  in  APL  are  either  monadic  or  dyadic 
meaning  one  or  two  arguments,  respectively,  are  required.  The 
same  symbol  may  mean  two  different  things  when  used  as  a  monadic 
or  dyadic  operator;  for  example,  if  A  is  a  number  and  B  is  a  vector, 
the  operator  %  in  the  monadic  expression  2  A  generates  a  vector  and 
in  the  dyadic  expression  B^A  generates  a  number. 

Basic  to  any  programming  language  is  the  ability  to 
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express  mathematical  relationships  precisely  and  the  ability  to 
conditional  branch.  Each  of  these  will  be  discussed  in  turn. 

The  basic  dyadic  symbols  for  the  five  elementary 
operations  are  -,  +,  x,  -,  and  *  where  the  latter  means  exponentiation. 
Parentheses  (  )  are  used  to  enclose  expressions  in  order  to  supercede 
the  right  to  left  rule.  The  dyadic  symbol^-  sets  the  left  argument 
equal  to  the  right  argument. 

The  conditional  branching  is  done  using  the  standard 
symbols  ^  ^  ^  and>  -enclosed  in  the  form  — >  (A  =  B)/STATMENT. 

This  means  that  if  A  equals  B,  the  next  statement  to  be  executed 
is  that  labelled  STATEMENT.  Otherwise,  the  next  statement  in  the 
program  is  executed  whether  it  is  labelled  or  unlabelled. 

The  ability  to  handle,  vectors  and  matrices  is  due  to 
several  other  operators  not  found  in  some  other  languages.  The 
number  of  elements  in  a  vector  A  isj)A;  and  if  A  has  been  previously 

defined, ^  A  is  automatically  set  equal  to  the  order  of  vector  A. 

•  • 

A? B  is  a  vector  of  order  A  and  all  elements  equal  to  B  while  2  A  is 

a  vector  of  order  A  and  elements  1,  2,  ...,  A.  Certain  elements 
are  picked  out  of  arrays  very  easily.  For  example,  f/A  and  L/A 
are  the  maximum  and  minimum  elements,  respectively,  of  A,  A  certain 
element  in  a  vector  is  located  by  the  expression  AtB  which  is  the 
index  of  the  first  element  of  A  equal  to  B.  A  £  <  r/A)  is  then  the 
index  of  the  maximum  element  of  a  vector . 

The  element  or  elements  of  a  vector  are  called  or 
identified  by  A[[  3 where  the  expression  in  the  brackets  is  the 
index.  It  may  be  one  number  or  an  expression  which  produces 
integers.  For  example,  A(V,3]  is  the  fifth  and  third  element  of 


. 


’ 


A,  respectively.  Also,  A  [l  + £  iQ  represents  the  3rd,  4th,  5th, 
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B,  B+  1,  B+  2  elements  of  vector  A. 

The  variables  in  APL  may  be  specified  as  either  local 
or  global  at  the  time  the  program  generating  them,  is  written.  A 
global  variable  label  exists  for  the  entire  workspace  it  is  in,  and 
execution  of  any  program  defining  a  variable  using  that  label  will 
redefine  it.  In  contrast,  if  a  variable  is  local  within  a  program 
and  the  same  label  is  used  for  a  global  variable  in  the  same  work¬ 
space,  executing  the  program  containing  the  local,  variable  does 
not  use  or  respecify  the  global  variable  in  the  workspace. 

An  example  of  local  variables  is  illustrated  by  the 
program  DIFFMEANS  (Appendix  A. 3. a.).  Variables  after  the  semicolon 
are  local,  that  is,  they  are  defined  and  used  within  the  function 
or  program  and  not  used  outsider  of  it. 

Any  program  can  be  written  as  a  monadic  or  dyadic 
operator  with  input  and  output  arguments  specified.  The  program 
GENFLT  (Appendix  A. 1  .a.)  illustrates  this  point.  The  input  arguments 
are  MAX  and  NUM,  and  the  command  20  GENFLT  50  is  equivalent  to 
specifying  MAX  =  20  and  NUM  =  50  and  then  executing  the  program. 

The  output  argument ‘ is  V;  and  if  the  result  is  to  be  stored  in 
RESULT,  RESULT  <c~  20  GENFLT  50  is  equivalent  to  executing  the 

program  and  then  setting  RESULT  =  V. 

It  is,  therefore,  very  easy  to  specify  the  input  and 
output  arguments  as  this  can  be  done  in  the  same  statement  thac 
orders  the  program  to  be  executed.  If  a  program  requires  more  than 
two  input  arguments,  the  two  labels  which  are  available  are 
specified  as  scalars,  vectors,  or  matrices  as  required.  An  example 
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of  this  is  DIFFMEANS  where  six  inputs  are  required  and  are  introduced 
via  the  two  labels  A  and  B.  A  and  B  are  both  vectors  of  order 
three  containing  three  input  arguments  each. 

4.1.4  Usefulness  of  APL 

The  APL  system  was  used  in  the  simulation  model  for  three 
main  reasons:  ease  of  access,  direct  link  with  the  computer,  and 
its  adaptability  to  the  model. 

1.  Ease  of  access.  Because  the  system  is  time  shared,  it  is 
possible  for  many  people  to  use  the  system  at  the  same  time. 
There  was  little  trouble  in  getting  adequate  computer  time. 

2.  Direct  link  to  computer.  This  permitted  the  type  of 
instant  error  feedback  and  instant  correction  feature  described 
in  Section  4.1.2.  This  was  useful  as  the  model  required  a 
large  number  of  programs  which  are  dependent  on  each  other. 

If  a  day  or  two  delay  occurred  between  an  error  and  its 
correction,  it  would  have  been  impossible  to  program  the 
model  in  the  time  allotted. 

3.  Adaptability  to  the  model.  A  main  feature  of  APL  is  its 
ability  to  handle  arrays  with  a  minimum  of  ''housekeeping1' 
programming  to  specify  the  arrays  and  flag  the  end  of  them. 
Section  4.2.2.c.ii  illustrates  this  point  in  the  explanation 
of  the  program  SUMARY. 

As  there  was  a  large  amount  of  data  to  be  handled  (the 
initial  rosters  contain  3,000  numbers)  and  similar  operations 
are  to  be  done  on  large  groups  of  data,  it  is  useful  to  store 
and  use  the  data  in  matrix  form.  The  matrix  is  then  processed 


.. 
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in  an  orderly  fashion,  and  this  language  is  easily  adapted 
to  processing  matrices. 

Several  disadvantages  in  the  system,  as  it  exists,  were 
also  found.  Since  the  system  is  only  in  operation  for  a  short  time 
each  day,  it  is  desirable  to  be  able  to  actually  use  the  computer 
as  much  as  possible  while  connected  to  the  system;  however,  data 
must  also  be  entered  into  the  system,  and  this  must  be  done  while 
it  is  in  operation.  There  is  no  provision  for  slow  entry  of  data 
on  a  relatively  inexpensive  machine  (such  as  a  keypunch)  which  gets 
it  into  a  form  which  can  be  rapidly  entered  into  the  main  computer; 
however,  the  model  was  designed  in  such  a  way  as  to  minimize  the 
data  handling  requirements . 

The  second  disadvantage  lies  in  the  small  workspace 
size  and  this  eliminated  the  Balas  algorithm  from  further  considera¬ 
tion.  While  the  algorithm  is  not  practical  due  to  inherent  difficulties 
(inefficient  when  matrix  size  exceeds  30  x  30),  it  could  have  been 
used  to  find  optima  which  would  be  useful  in  checking  the  results 
from  the  heuristic  rules.  It  may  be  possible  to  store  10,000  to 
15,000  numbers  in  one  workspace  while  the  Balas  algorithm  requires 
a  main  matrix  of  over  100,000  numbers  for  the  size  of  problem  con¬ 
sidered  here.  (See  Section  2.2.2.) 

In  spite  of  these  difficulties,  the  system  was  considered 
to  be  useful  for  this  model,  and  the  heuristic  rules  were  simulated 
in  APL  and  the  results  obtained  from  this  system. 

4 . 2  Programs  Used  in  the  Simulation 


The  programs  used  in  this  model  fall  inio  three  categories. 
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1.  Programs  which  generate  and  process  the  input  data  (rosters). 

2.  Programs  which  schedule  the  rosters  and  compute  summary 
statistics  from  the  schedules. 

3.  Programs  which  process  the  summaries  generated  above. 

Each  set  of  programs  will  be  discussed  in  turn  in  Sections 
4.2.1,  4.2.2,  and  4 .-2.3. 

4.2.1  Data  Generators  and  Processors 

•  The  data  generators  are  the  routines  GENPOISSON  and 
GENFLT .  The  processors  are  CHECK,  STRINGOUT,  SUMMARY,  and  TABFREQ . 
These  programs  are  displayed  in  Appendix  Al. 

GENPOISSON  and  GENFLT  generate  a  vector  of  order  L  and 
mean  M  or  a  vector  of  order  L  and  range  1  to  R  when  the  commands 
L  GENPOISSON  M  or  L  GENFLT  R  are  entered  into  the  system  with  the 
workspace  REGINA  loaded  into  the  active  workspace.  These  programs 
are  independent  of  any  other  privately  defined  or  public  library 
functions;  therefore,  they  are  executable  without  any  other  functions 
defined  or  loaded  into  the  workspace. 

The  summary  characteristics  of  the  rosters  are  computed 
by  typing  SUMMARY  followed  by  F3,  4,  or  5;  or  P3,  4,  or  5,  depending 
on  whether  we  wish  to  compute  the  characteristics  of  the  uniform  or 
the  Poisson  distributed  rosters  with  mean  3,  4,  or  5 .  For  example, 
SUMMARYP5  computes  the  characteristics  of  the  ten  Poisson  rosters 
with  mean  5.  The  output  is  described  in  Section  4.3.3. 

To  use  this  routine,  two  conditions  must  be  met: 

1.  The  rosters  must  be  in  appropriate  form 

2.  Several  other  programs  must  be  defined 
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For  point  1,  the  roster  capacities  must  be  in  a  matrix 
labelled  ROSC  and  durations  in  a  matrix  labelled  ROSD  followed  by 
the  same  suffix  as  SUMMARY,  This  matrix  must  have  as  many  columns 
as  elements  in  the  roster  (25)  and  as  many  rows  as  rosters  (10). 
This  is  elaborated  in  Section  4.3.3. 

For  point  2,  conditions  defined  by  the  hierarchy  of 
programs  stated  in  Table  44  must  be  met. 

TABLE  44 

Hierarchy  of  Programs  in  Workspace  SASK 


1 

!  CHECK 

i  RND 

1 

I  STRINGOUT 

SUMMARY 

. 

i  TABFREQ 

CHECK 

v' 

RND 

STRINGOUT 

SUMMARY 

✓ 

✓ 

TABFREQ 

A  check  (i/)  indicates  that  the  program  on  the  vertical 
axis  must  be  defined  in  order  to  execute  the  program  on  the  hor¬ 
izontal  axis.  This  convention  will  be  kept  in  future  hierarchy 
tab les  . 

The  programs  STRINGOUT  and  TABFREQ  are  independent  of 
any  other  programs,  and  no  others  depend  on  them  as  they  have  no 
(w')  in  either  vertical  or  horizontal  columns.  RND  is  independent 
of  any  others  but  others  are  dependent  on  it  (checks  in  horizontal 
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column  only) .  SUMMARY  is  dependent  on  other  programs  but  none  are 
dependent  on  it  while  CHECK  is  both  dependent  on  other  programs  and 
others  are  dependent  on  it . 

The  rosters  are  checked  using  the  program  SUMMARY  which 
feeds  each  roster  in  turn  to  the  program  CHECK  by  allocating  each 
row  of  the  roster  matrices  to  the  input  arguments  of  CHECK.  It 
accumulates  the  summary  data  generated  by  CHECK  in  a  matrix  SUMR 
followed  by  the  same  suffix  as  follows  SUMMARY. 

CHECK  computes  the  ranges,  averages,  and  standard 
deviations  for  the  roster  assigned  to  it  and  using  RND  rounds  these 
off  to  a  suitable  number  of  places  of  decimal.  It  uses  the  program 
TABFREQ  to  compute  how  many  of  each  element  there  are  in  the  roster 
and  prints  the  roster.  It  then  passes  control  back  to  SUMMARY 
which  stores  the  results  in  SUMR  and  assigns  it  another  roster 
which  is  the  next  row  of  the  matrices  containing  the  rosters.  The 
n'th  row  of  SUMR  contains  the  characteristics  of  the  roster  repres¬ 
ented  by  the  n’th  row  of  the  roster  matrix.  When  all  the.  rosters  are 
processed,  SUMMARY  prints  the  summary  characteristics  and  passes 
control  back  to  the  keyboard. 

4.2.2  Main  Scheduling  Routine 
4 . 2 . 2 . a  Main  Routine 

The  scheduling  of  the  roster  and  computation 
of  the  summary  characteristics  for  each  roster  is  coordinated 
by  a  program  MAIN  which  accepts  definition  of  the  data 
(roster  and  shop),  incorporates  the  scheduling  routines, 
compiles  data  from  them,  and  prints  the  summary  statistics 


(STATS) . 
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4 . 2 . 2 . .  b  Primary  Scheduling  Programs 

The  programs  which  incorporate  the  fourteen 
scheduling  rules  are  called  RANDOMS,  RANDOMFITS,  SSRANDOMFIT, 
LSRANDOMFIT ,  SPFITS,  LPFITS,  SPSCFITS,  SPLCFITS,  LPSCFITS, 
LPLCFITS ,  SSSPFIT,  SSLPFIT,  LSSPFIT,  and  LSLPFIT .  These  are 
the  program  names  of  the  scheduling  rules  described  in 
Section  1.5.1  and  are  in  the  same  order  as  they  are  defined 
in  Table  3.  lable  45  shows  the  dependence  of  these  routines 
on  the  sorter,  scheduler,  and  summary  routines  which  will  be 
described  in  Section  4.2.2.C. 

TABLE  45 

Dependence  of  Primary  Scheduling  Programs 


1 

!  ASGNXT 

1  ASGNAHD 

DECSORT 

1 

!  ASSORT 

j 

I  DECSORT  3 

r— - - 

\ 

1  ASSORT  3 

& 

CO 

RANDOMS 

✓ 

RANDOMFITS 

S 

\s 

SS RANDOMFITS 

IS 

V 

LSRANDOMFIT 

s 

is 

V' 

SPFITS 

\S 

\S 

LPFITS 

IS 

SPSCFITS 

IS 

s' 

SPLCFITS 

*S 

s 

V' 

LPSCFITS 

s 

\S 

\s 

LPLCFITS 

i' 

is 

SSSPFIT 

V 

iS 

\s 

IS 

SSLPFIT 

*s 

IS 

S 

LSSPFIT 

kS 

IS 

is 

LSLPFIT 

IS 

jarA 
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The  seven  routines  upon  which  the  fourteen 
scheduling  programs  are  dependent  fall  into  three  categories 
and  are  described  in  Section  4.2.2.C. 

1.  Scheduler  -  ASGNXT,  ASGNAHD 

2.  Sorter  -  DECSORT,  ASSORT,  DECS0RT3 ,  ASS0RT3 

3 .  Summary  -  SUMARY 

The  sequence  of  operations  in  a  typical  scheduling 
routine  is  to  accept  the  input  shop  and  roster  defined  by 
MAIN,  rearrange  the  roster  using  0,  1,  or  2  sorters  routines, 
schedule  the  jobs  using  one  of  the  schedulers  and  compute 
summary  statistics.  As  an  example,  the  method  is  described 
in  detail  for  one  scheduling  routine,  LSSPFIT,  which  assigns 
jobs  with  priority  to  largest  size  and  breaks  ties  according 
to  smallest  duration  first. 

FIGURE  10 
LSSPFIT 

Largest  Size  -  Smallest  Duration 
 Scheduling  Program 


. — 1 

f— 1 

l _ 1 

CAP  «c— 

INPCAP 

i — i 

1 — ! 

DECSORT  3 

[2] 

DUR  <- 

INPDUR 

! - 1 

00 

t _ 1 

X  <—  BB 

C3] 

NR  JOB  <- 

—  INPNRJOB 

f — 1 

1 _ 1 

P  << —  AA 

! — 1 

! — J 

ASSORT 

D°'J 

C  <-  SHOP 

rsD 

DUR 

AA 

Qn] 

ASGNAHD 

re] 

CAP  <^— 

BB 

[12] 

SUMARY 
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Statements  1,  2,  and  3  specify  the  roster  from 
the  INP  variables  defined  by  the  program  MAIN  as  global 
variables.  In  this  way,  each  of  the  fourteen  scheduling 
routines  including  this  one  is  guaranteed  an  unaltered  set 
of  input  data. 

Statement  4  sorts  the  roster  in  ascending 
duration  order  as  required  by  LSSP.  Note  that  this  is  the 
tie-breaking  sort  and  is  done  first.  Statements  5  and  6 
place  the  roster  capacity  and  duration  vectors  which  result 
from  M  in  suitable  order  for  statement  7. 

Statement  7  sorts  the  roster  in  descending 
size  order  which  is  the  primary  sort  required  by  LSSP . 

Because  the  sorting  method  chooses  the  first  number  in  the 
sequence  in  case  of  ties,  the  second  sort  has  already  been 
done  according  to  the  specified  tie-breaking  rule  (in  statment 
4) .  Statements  8  and  9  put  the  roster,  which  is  now  in  order 
according  to  the  required  priority,  into  the  input  arguments 
of  the  scheduler  routine;  and  £l0^]  specifies  the  shop  the 
same  as  that  specified  by  MAIN. 

Statement  11  is  the  scheduler  routine  which 
incorporates  the  fitting  procedure  and  produces  the  schedule. 
Statement  12  computes  summary  statistics  from  the  schedule 
produced. 

4 . 2 . 2 . c  Dependent  Sorter,  Scheduler,  and  Summary  Routines 

As  indicated  by  Table  45  in  4.2.2.b,  these 
routines  are  necessary  for  the  execution  of  the  main  scheduling 


routines  . 


Each  of  these  will  be  described  in  turn. 
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4 . 2 . 2 . c . i  Sorter  Routines 

These  consist  of  four  programs  - 
DECSORT,  ASSORT,  DECS0RT3 ,  and  ASS0RT3 ,  none  of 
which  are  dependent  on  other  predefined  functions. 

The  input  arguments  are  vectors  DUR  and  CAP,  and 
the  output  is  vectors  AA  and  BB.  These  contain  the 
same  elements  as  inputs  DUR  and  CAP  but  are  rearranged 
in  some  specific  order.  If  DUR  Ci]  becomes  AA  m , 
CAP  Li]  will  become  BB  Lj]  for  all  I  and  J. 

The  routines  ASSORT  and  DECSORT 
rearrange  the  input  vector  DUR  in  ascending  or 
descending  order,  respectively,  while  ASS0RT3  and 
DECS0RT3  generate  a  new  vector  SIZE  such  that 
SIZE  [i]  =  DUR  [i]  x  CAP  CO  for  all  I  and  then 

sort  the  vector  SIZE  in  ascending  or  descending 
order,  respectively. 

Again,  if  DUR  [  i]  x  CAP  [i]]  generate 
SIZE[l]and  SIZE  DG  becomes  SIZE  CjU  after  sorting, 
then  DUR  f_Jj  and  CAP  C  iJ  become  AA  CJ3  and  BB  Cd  1 
for  all  I  and  J.  The  sequence  is  shown  in  Figure  10. 


FIGURE  10 
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4 . 2 . 2  .  c . i i  Summary  Routine 

This  is  one  program  SUMARY  which 
computes  the  nine  measures  of  effectiveness  according 
to  the  formulae  outlined  in  Section  3.2.1,  puts  them 
in  a  vector  called  SUM,  and  prints  this  vector.  (In 
addition,  MAIN  adds  the  vector  SUM  calculated,  for 
each  schedule  to  a  matrix  which  is  further  worked 
into  the  ranked  summary  statistics.)  This  program 
is  independent  of  any  others. 

Figure  11  compares  the  APL  statements 
in  Section  3.2.1  to  more  conventional  notation 
showing  the  limits  of  summation  over  the  arrays. 

Figure  11 

Comparison  of  APL  Summary  Calculations  to 
More  Conventional  Notation 

J. 

SUM  Cl]  =  mean  start  time  =  (-t/NUM^  f  0  NUM  =  ^  NUM  i  =  MS 

l  IfJ _ 

J 

where  NUM  i  is  start  time  of  i'th  job  J  is  number  of  jobs  in 
the  roster. 

SUM  £2j  =  time  shop  is  occupied  =  P/F  =  Max.  F ^ 

1£  i^J 

where  F^  is  the  finish  time  of  the  i'th  job. 

SUM  ['.3]  =  average  jobs  in  process  =  (+/J0B)  7  (P/F)  - 

N_ 

L  J0Bi 

i=l _ 

Max .  F 


where  J0Bi  is  the  number  of  jobs  running  in  time  i 
N  is  the  number  of  time  periods  any  job  is  running 
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SUM  [4]]  =  average  efficiency  =  +/(P  x  X)  ~  (CC  x  f/F) 

J 

PiXi 

_  i=l _ _ 

CC • (Max  Fi) 

14i4J 

where  Pj_  and  X-^  are  the  durations  and  capacities  of  the  i'th  job 
CC  is  the  initial  shop  capacity  in  each  period 
SUM  [jf]  =  mid-range  efficiency  =  1  -  SLACK  -  CC  x  2  x  BEGIN 


=  1  = _ SLACK _ 

2  * (CC)  •  (BEGIN) 

where  BEGIN  =  ((F/F)  -  (4  I  F))  7  4 

=  next  lowest  max.  F^  that  is  divisible  by  4 
divided  by  4 

SLACK  =  +  /C  ChEGIN  +  l2  x  BEGIN!] 

2 -BEGIN 


C£  ,  Ci  =  unused  shop  capacity  in 

period  i 

i  =  BEGIN  +  1 

The  notation  A  I  B  represents  the  remainder  from  the  operation 
B  ;■  A;  e  .g. ,  4  19  =  1 

SUM  [6  1  =  weighted  delay  factor  =  +/PxXx(NUM  -  1)  = y  Pi^i^i 

i=l 

where  d^  is  the  delay  of  job  i 


SUM  £7j  ~  total  start  time  =  +/NUM  -  J  NUM^ 

i=l 

SUM  [g]  =  start  time  standard  deviation  = 


where  MS  is  mean  start  time 


defined  in  SUM  [1] 


. 


. 
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The  expression  (^NUM)  MS  develops  a  vector  of  length 
NUM  =  J  and  each  element  MS.  The  subtraction  NUM  -  MS 
would  not  be  possible  directly  as  these  vectors  must  have 
equal  lengths  for  an  element  by  element  subtraction  to  be  made. 


SUM  C9D  =  start  time  deviation  about  0  =  ((+/NUM*2) 


V  J 

4 . 2 . 2 . c . i i i  Scheduling  Routines 

These  are  the  key  routines  in  the 
model  and  are  dependent  on  three  routines  which 
provide  more  detail  about  the  generated  schedule  which 
will  be  described  in  Section  4.2.2.d. 


ASGNXT  and  ASGNAHD  require  as  inputs 


a  vector  representing  shop  capacity  (C)  and  vectors 
representing  the  durations  (P)  and  capacities  (X)  of 
the  jobs  in  the  roster.  The  routines  will  assign 
each  job  i  (Pi,  Xi)  to  the  shop  so  that  the  sum  of 
the  capacity  requirements  of  all  jobs  running  does 
not  exceed  the  shop  capacity  at  any  time.  The  output 
is  a  schedule;  i.e.,  a  vector  of  start  times  (NUM) 
and  finish  times  (F)  for  the  roster,  a  vector  (C) 
indicating  unused  shop  capacity  in  each  time  period 
and  a  vector  JOB  indicating  the  number  of  jobs  in 
process  in  each  time  period. 

The  routine  ASGNXT  schedules  the  jobs 
in  the  same  order  as  they  appear  in  the  input  roster. 


Ill 

At  the  first  time  period  t  =  1,  the  first  A  jobs 
such  that: 

A  A  +  1 

X  *i  4CX  and  XI  >  Cx 
i=l  i=l 

are  scheduled.  Then  the  clock  is  advanced  to  t  =  2 
and  the  A  +  l'th  job  is  tried.  If  there  is  less 
unoccupied  capacity  at  t  =  2  than  the  requirement 
of  the  A  +  l'th  job,  the  clock  is  again  advanced; 
and  the  A  +  l'th  job  is  tried  again.  This  continues 
until  it  is  scheduled.  This  procedure  is  repeated 
with  each  succeeding  job  in  the  roster  in  order 
until  they  are  all  scheduled.  Figure  12  illustrates 
the  procedure  diagramat ically . 
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FIGURE  12 

Flowchart:  for  ASGNXT 
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The  variable  NN  is  introduced  to  mark 
the  time  period  in  which  each  assingment  occurs. 

When  the  unused  capacity  is  decremented  beyond  the 
period  of  time  in  which  the  job  is  initially 
assigned,  N  is  incremented;  but  we  wish  to  go  back 
to  the  time  period  the  job  is  assigned  to  see  if 
the  next  job  will  fit.  Otherwise,  only  one  job 
would  be  assinged  in  any  time  period.  Figure  13 
illustrates  this  method  of  keeping  track  of  the  base 


t  ime 
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FIGURE  13 

Decrementing  Procedure  for  N 
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The  output  format  from  ASGNAI1D  is 
the  same  as  ASGNXT,  but  the  scheduling  method  is 
different.  This  routine  is  the  basis  of  the  fitting 
procedure  mentioned  several  times  previously. 

An  attempt  is  still  made  to  schedule 
the  jobs  in  the  order  they  appear  in  the  input 
roster;  however,  if  a  job  will  not  fit,  a  search 
ahead  in  the  roster  is  made  until  the  shop  is  filled 
or  the  roster  exhausted.  Then  the  clock  is  advanced 
to  the  next  time  period,  and  the  first  job  which 
was  encountered  last  time  period  but  was  not  scheduled 
is  tried  again. 

If  the  situation 

A_  Ajy  1 

^  X.j_^-Cp  and  Xi>  Ci 

i=l  i-1 

occurs,  we  try 

A 

YL  xi  +  xA+2  ?C 
i=l 

If  ?  is  Z.  ,  we  schedule  job  A  +  2  and  consider  XA+^ 
in  the  same  manner  as  above.  If  ?  is  equal,  we 
schedule  job  A  +  2,  increment  the  clock,  and  try  with 
job  A  +  1.  If  ?  is  we  ignore  job  A  +  2  and 

try  job  A  +  3.  Jobs  A  +  2,  A  +  3,  - -  J  are  tried 

in  this  manner;  and  when  job  J  is  reached,  the  clock 
is  incremented,  and  the  first  unassigned  job  in 

the  roster  tried  again. 

To  avoid  having  jobs  assigned  more 


■ 
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than  once,  each  job  is  initially  marked  with  a  0; 
and  when  it  is  scheduled,  the  mark  is  changed  to  a 
1.  Before  the  procedure  above  is  tried,  the 
appropriate  element  of  the  marking  vector  is  tested 
for  0  or  1.  If  it  is  0,  an  attempt  is  made  to 
schedule  the  job;  and  if  it  is  1,  the  job  has  already 
been  scheduled  and  the  next  job  in  the  roster  is 
considered.  Figure  14  illustrates  the  ASGNAHD 
algorithm. 


FIGURE  14 


Flowchart  for  ASGNAHD 
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4 . 2 . 2 .  d  Roster  Detail  Routines 

These  routines  can  be  added  to  ASGNXT  and 
ASGNAHl)  to  provide  more  detail  about  the  schedule  than  its 
summary  characteristics  produced  by  the  program  SUMARY 
described  in  4.2.2.c.ii. 

Statement  Cll  of  either  ASGNAHD  or  ASGNXT 
can  be  either  a  — ->2  or  call  for  program  HEADING.  Statement 
[jL6U  of  ASGNXT  can  be  a-— >11  or  call  for  program  DISPLAY 
and  statement  \ll]  of  ASGNAHD  can  be — >18  or  DISPLAY. 

If  the  —>alternat ive  is  chosen,  only  the 
summary  characteristics  of  the  schedule  are  printed;  if 
HEADING  and  DISPLAY  are  used,  the  schedule  is  printed  in  the 
order  the  jobs  are  scheduled  including  the  start  and  finish 
time,  capacity,  and  duration  of  all  the  jobs  in  the  roster. 
The  routine  ASGNDETL  is  attached  at  the  end  of  either  ASGNAHD 
or  ASGNXT  and  causes  the  unused  capacity  and  number  of  jobs 
running  in  each  time  period  to  be  printed. 

4 . 2 . 2 .  e  Summary  Ranking  Routine 

The  program  MAIN  generates  a  matrix  MAT  con¬ 
sisting  of  all  the  summary  data  produced  by  SUMARY  after  each 
of  the  fourteen  scheduling  routines  is  applied  to  the  roster. 
The  ranking  routine  STATS  applies  the  ranking  method  of 
Section  3.4.1  to  this  9  x  14  matrix.  STATS  scans  each  of 
the  nine  columns  in  turn  and  picks  out  the  maximum  (MAX) 
and  minimum  (MIN)  element  in  each  and  finds  their  difference 
(DIFF).  For  each  element  in  the  column,  ELEMENT  -  MIN  is 


DIFF 
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calculated ,  rounded  to  four  places  using  RND,  and  multiplied 
by  10,000.  The  element  is  then  replaced  by  this  quantity 
resulting  in  a  9  x  14  matrix  of  rankings  from  0  to  10,000. 

The  hierarchy  of  programs  used  in  the  model  to  compute 
schedules  and  print  summary  data  for  each  roster  is  given  in 
Table  46.  A  ( \/)  in  the  horizontal  row  beside  any  program 
indicates  it  is  dependent  on  the  program  in  the  top  heading. 
The  hierarchy  shows  program  dependence  only;  it  may  be  that 
a  particular  routine  requires  data  from  another  to  give 
meaningful  results.  This  is  not  indicated  as  all  programs 
then  would  be  interdependent. 


TABLE  46 

Hierarchy  of  Programs  in  Workspace  REGINA 


MAIN 

RANDOMS 

CO 

£-i 

w 

Fn 

g 

n 

Z 

g 

jSSRANDOMFIT 

LSRANDOMFIT 

SPFITS 

LPFITS 

jSPSCFITS 

CO 

H 

H 

Fu 

CJ 

Z 

pc, 

CO 

I  LPSCFITS 

LPLCFITS 

SSSPFIT 

ISSLPFIT 

LSSPFIT 

LSLPFIT 

ASGNXT 

ASGNAHD 

Idecsort 

{assort 

IDECS0RT3 

ASS0RT3 

SUMARY 

STATS 

HEADING 

K— ' 

Pu 

CO 

1—1 

i - • 

H 

W 

Q 

z 

c 

CO 

< 

RND 

MAIN 

✓ 

is 

4/ 

U 

\S 

Is 

is 

Vs 

IS 

Is 

IS 

l' 

t- 

IS 

is 

is 

O 

o 

p 

RANDOMS 

IS 

is 

RANDOMFITS 

v' 

ts 

SS RAN DOME IT 

t/ 

V 

is 

LSRAND0MFIT 

V 

V' 

is 

SPFITS 

l' 

V 

Is 

LPFITS 

l* 

•s 

is 

SPSCFI.TS 

is 

V 

I- 

SPLCFITS 

is 

s 

iS 

w* 

LPSCFITS 

V' 

is 

V 

LPLCFITS 

»' 

V" 

c- 

SSSPFIT 

is 

l' 

\S 

\s 

SSLPFIT 

i' 

\s 

U" 

t- 

LSSPFIT 

IS 

s 

\S 

V- 

LSLPFIT 

IS 

S 

ASGNXT 

o 

o 

o 

ASGNAHD 

© 

o 

DEC SORT 

1 

TABLE  46  (Cont'd.) 
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MAIN 

RANDOMS 

|  RANDOMFITS 

H 

w 

>d 

o 

Q 

a 

s 

cn 

cn 

LSRANDOMFIT 

!  SPFITS 

|  LPFITS 

!  SPSCFITS 

|  SPLCFITS 

1  LPSCFITS 

1  LPLCFITS 

SSSPFIT 

SSLPFIT 

!  LSSPFIT 

[  LSLPFIT 

!  ASGNXT 

ASGNAIID 

DECSORT 

i  ASSORT 

DECS0RT3 

|  ASS0RT3 

j  SUMARY 

STATS 

!  HEADING 

DISPLAY 

i  ASGNDETL 

RND 

ASSORT 

DEC SORT 3 

AS SORT 3 

SUMARY 

STATS 

\/ 

HEADING 

DISPLAY 

ASGNDETL 

RND 

(v^)  compulsory  dependence 
0  optional  dependence  (see  Section  4.2.2.d) 

4.2.3  Sumary  Processing  Programs 

The  key  program  in  this  group  of  programs  is  labelled 
SYN  prefixed  by  A,  B,  C,  or  D,  corresponding  to  the  four  measures 
of  effectiveness.  This  is  the  program  which  generates  the  chart 
of  means  and  standard  deviations  which  is  displayed  in  Table  14  of 
Section  3.4.3  and  is  the  main  result  of  this  experiment. 

The  program  SYN  accpets  the  rankings  computed  by  STATS 
for  each  scheduling  rule  for  each  measure  of  effectiveness  in  matrix 
form,  essentially  a  60  x  14  matrix  for  each  measure  of  effectiveness 
This  corresponds  to  sixty  rosters  and  fourteen  scheduling  rules  run 
together.  The  output  is  a  matrix  MSD(  )T0TAL  where  (  )  contains 
the  same  letter  A,  B,  C,  or  D  as  prefixed  SYN  corresponding  to  the 
measure  of  effectiveness.  For  example,  BSYN  computes  the  mean  and 
standard  deviation  of  the  shop  occupied  time  (measure  B)  for  the 
sixty  rosters  for  each  of  the  fourteen  scheduling  rules  and  places 
the  number  of  items  (60),  mean  and  standard  deviation,  in  a  3  x  14 


matrix  MSDBTOTAL. 
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SYN  is  dependent  on  a  program  MSD  which  performs  the 
actual  calculations  required.  The  statement  MSD  VECTOR  yields 
a  vector  consisting  of  the  number  of  elements  in  VECTOR,  their  mean, 
and  their  standard  deviation  calculated  according  to 


N 


]>_  (VECTOR^  -  VECTOR) 2 
i=l 


N 


where  N  is  the  number  of  elements  in  VECTOR 


VECTOR  is  the  mean  of  the  elements  in  VECTOR 


The  chart  of  t  values  for  the  differences  between  means 
in  MSD(  )TOTAL  is  generated  by  a  program  labelled  C  Oil?  A  RE .  This 
compares  the  first  row  of  MSD(  )TOTAL  to  the  other  rows  in  turn, 
calculating  t  values  for  the  differences  between  means  then  compares 
the  second  row  to  all  subsequent  rows  and  so  on  until  all  possible 
pairs  of  means  are  compared. 

This  program  is  dependent  on  a  routine  DIFFMEANS  and  is 
really  an  orderly  method  of  feeding  the  input  arguments  to  DIFFMEANS. 
This  program  incorporates  a  routine  for  calculating  a  t  value  from 
Nl’  si>  xl>  N2,  S2,  and  X2  where  N  is  the  number  of  items  in  the 
sample;  S  is  the  standard  deviation,  and  the  X's  are  the  means  for 
which  we  are  trying  to  establish  significance  of  difference. 

DD 

[83  7-5.5 

[93  7-5.8 


ine  lormuiae  aic. 


Z  =  Ni  +  No  -  2 


NjSj2  +  N2S2‘ 


z 


sD  - 


S  /  Nl  +  N2 

NfN2 


Xj  -  x2 
"  Sn 


T 


D OD  7-5.1 
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where  the  numbers  in  £  3  indicate  the  statement  number  in  DIFFMEANS 
(see  Appendix  A. 3. a)  and  following  is  the  formula  number  in  the 
reference .  D1 

The  program  RANK  is  independent  of  other  programs  and 
there  is  no  routine. dependent  on  it.  Given  an  initial  vector,  this 
will  arrange  the  elements  in  ascending  order  and  print  only  the 
indices  which  the  elements  had  in  the  initial  vector  in  this  order. 

The  input  is  a  row  of  a  matrix,  and  each  row  is  automatically  processed 
in  turn.  The  procedure  is  illustrated  in  Figure  15.  VECTOR  is  the 
vector  to  be  ranked. 
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FIGURE  15 

Flowchart  for  RANK 
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4 . 3  Data  Handling 

4.3.1  Workspace  Allocation 

The  final  model  is  contained  in  three  workspaces  con¬ 
taining  programs  and  data.  These  are  defined  as  follows  where  the 
block  letters  indicate  the  workspace  name. 

1.  REGINA  -  contains  the  scheduling  routines  described  in 
Section  4.2.2,  the  roster  generating  routines  and  rosters  in 
process . 

2.  SASK  -  contains  the  routines  for  checking  rosters,  the 
rosters  and  their  summary  characteristics. 

3.  YORKTON  -  contains  the  t-test  and  ranking  routines.  The 
resulting  rankings  for  the  four  measures  of  effectiveness  are 
also  stored  here  as  well  as  their  means  and  standard  deviations 
for  the  group  of  sixty  rosters. 

4.3.2  Data  Flow 

The  data  is  generated  in  the  workspace  REGINA  from  the 
roster  generator  programs.  From  this,  there  is  a  two-way  branch; 
the  rosters  are  stored  in  another  workspace  (SASK)  due  to  lack  of 
storage  capacity  in  the  initial  workspace,  and  the  summary  results 
are  dumped  on  paper  while  each  roster  is  processed.  The  summary 
rankings  for  the  four  measures  of  effectiveness  chosen  are  then 
re-entered  into  the  system  in  another  workspace  (YORKTON).  This 
could  possibly  have  been  done  automatically  while  each  roster  was 
processed;  however,  the  need  was  not  anticipated  at  that  time. 

The  data  flow  in  each  workspace  was  well  automated.  In 
SASK,  the  rosters  were  passed  automatically  through  their  checking 
routine  and  the  results  printed  and  stored  for  future  access.  In 
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YORKTON ,  the  mean  and  standard  deviation  were  computed  from  the 
summary  data  and  printed  and  stored  in  matrices.  The  result  of  the 
t-test  was  dumped  on  paper  only;,  however,  to  retrieve  this,  the 
user  has  only  to  execute  the  program  (  )SYN  where  (  )  contains 
either  A,  B,  C,  or  D,  depending  on  the  measure  of  effectiveness 
desired . 

4.3.3  Data  Storage 

Advantage  was  taken  of  the  ability  of  APL  to  handle 
matrices  by  storing  and  handling  data  in  matrix  form.  An  exception 
to  this  is  the  initial  processing  of  each  roster  through  the  scheduling 
routines.  In  this  case,  each  roster  was  processed  individually 
since  the  operation  frequently  took  ten  to  twenty  minutes  per  roster. 

As  the  allowable  session  on  a  terminal  is  twenty  minutes,  it  is 
pointless  to  have  a  system  which  automatically  proceeds  to  the  next 
roster.  However,  a  minor  modification  to  the  input  part  of  MAIN 
would  enable  this  procedure  to  be  automated. 

In  general,  the  model  is  arranged  so  that  any  procedure 
will  take  from  five  to  ten  minutes  and  then  give  control  back  to 
the  operator.  If  the  time  was  any  shorter,  a  high  proportion  of 
the  terminal  session  would  be  spent  on  "housekeeping"  and  if  it 
were  longer,  there  would  be  a  lack  of  error  control  and  the  frequent 
nuisance  of  suspending  programs  when  the  session  is  over  resulting 
in  a  chopped  up  output. 

When  the  rosters  were  transferred  to  SASK,  they  were 
put  into  matrices  which  were  identified  by  the  letters  ROS .  The 
capacities  and  durations  were  put  in  separate  matrices,  and  each 
set  of  two  matrices  contained  the  ten  rosters  which  were  generated 
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by  the  same  process.  Table  47  shows  the  names  of  the  vectors  in 
which  the  rosters  are  stored. 


TABLE  47 


Roster  Storage 

Vector  Names 

Roster  Class 

Capacities 

Durat ions 

Uniform  distributed 

1  to  5 

R0SCF3 

R0SDF3 

Uniform  distributed 

1  to  7 

R0SCF4 

R0SDF4 

Uniform  distribution 

1  to  9 

R0SCF5 

R0SDF5 

Poisson  distribution* 

Mean  3 

R0SCP3 

R0SDP3 

Poisson  distribution* 

Mean  4 

R0SCP4 

R0SDP4 

Poisson  distribution* 

Mean  5 

R0SCP5 

R0SDP5 

The  coding  in  Table  47  is 

interpreted  as  follows: 

ROS 

X 

A 

Y  Z 

*  4' 

/ 

Indicates  a  roster' 

/ 

\  \ 

\  Mean  of  d 

istribut ion 

1 

C  -  capacity,  D  -  duration 

\ 

F  -  uniform,  P 

-  poisson 

*  See  Section  3.1 

All  roster 

matrices 

are,  therefore,  of  size 

10  rows 

(representing  ten  rosters)  and  25  columns  (representing  25  jobs  per 
roster) . 


The  command  SUMMARY  Y(YZ)  automatically  computes  the 

summary  characteristics  of  the  group  of  ten  rosters  ROSC  (YZ)  and 

D 

Pi  aces  the  summary  in  a  matrix  SUMR(YZ) .  The  capacity  average, 
standard  deviation,  maximum,  minimum;  duration  average,  standard 
deviation,  maximum  minimum;  and  roster  work  required  are  stored  in 
the  rows  of  these  matrices  and,  therefore,  the  matrix  size  is  ten 
rows  (representing  ten  rosters)  by  nine  columns  (representing  the 
nine  items  above) . 
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The  summary  data  for  each  roster  generated  by  the  main 
processing  routine  was  entered  in  YORKTON  in  matrix  form.  Each 
matrix  represents  ten  rosters  in  the  same  groupings  as  before  and 
one  measure  of  effectiveness.  Each  matrix  contains  fourteen  rows 
representing  the  fourteen  scheduling  rules  and  ten  columns  representing 
ten  rosters -in  the  group.  The  data  is  stored  in  matrices  labelled 
as  in  Table  48. 


TABLE  48 


Roster 

Summary  Dat 

a  Matrix 

Storage 

Labels 

D 

Type 

Me 

asure  of 

Effect iveness 

A 

B 

C 

Uniform 

1  to 

5 

AF3 

BF3 

CF3 

DF3 

Uniform 

1  to 

7 

AF4 

BF4 

CF4 

DF4 

Uniform 

1  to 

9 

AF5 

BF5 

CF5 

DF5 

Poisson* 

Mean 

3 

AP3 

BP  3 

CP  3 

DP  3 

Poisson* 

Mean 

4 

AP4 

BP4 

CP4 

DP4 

Poisson* 

Mean 

5 

AP5 

BP5 

CP5 

DP  5 

*  See  Section  3.1 

The  ranking  matrix  for  any  roster  is  contained  in  these 
summary  matrices.  For  example,  if  we  wanted  the  ranking  matrix  for 
the  seventh  Poisson  distributed  roster  of  mean  4,  we  would  take  the 
seventh  column  of  AP4,  BP4,  CP4,  and  DP4  •• 

The  average  rank  of  each  scheduling  rule  according  to 
each  measure  of  effectiveness  is  stored  in  the  second  column  of  the 
matrices  MSD(W) TOTAL  where  W  =  A,  B,  C,  or  D,  depending  on  the  measure 
of  effectiveness.  The  first  column  of  this  matrix  indicates  the 
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number  of  items  in  the  sample;  and  the  third  is  their  standard 
deviation.  These  are  the  matrices  displayed  in  Table  14. 


A  summary  of  the  data  in  storage  is  then: 

Location 

Data  Matrix  Name 

SASK 

Rosters  ROS  (X  YZ) 

SASK 

Roster  Summary  SUM  R(YZ) 

YORKTON 

Rankings  (WYZ) 

YORKTON 

Ranking  Summary  MSD(W)TOTAL 

where  the  variables  denote  the  following  as  mentioned  above 


W  - 

measure  of  effectiveness 

X  - 

capacity  or  duration 

Y  - 

poisson  or  uniform 

Z  -  mean  of  the  distribution 


CHAPTER  5 


FURTHER  RESEARCH 

In  this  paper,  only  a  very  special  case  of  a  wide  variety  of 
scheduling  problems  was  studied;  and  practical  rules  of  thumb  were 
developed  for  dealing  with  this  case. 

However,  the  programs  are  sufficiently  general  that  with  some  mod¬ 
ifications,  other  studies  could  be  undertaken.  The  data  storage  system 
is  such  that  comparisons  of  performance  of  the  same  rule  operating  on 
rosters  generated  in  different  ways  could  be  explored.  This  was  initially 
attempted,  but  the  sample  sizes  were  too  small  to  get  significant  results. 
It  may  be  that  some  rules  will  show  greater  or  lesser  effectiveness 
depending  on  whether  the  rosters  are  Poisson  or  uniformly  distributed  or 
are  generated  by  some  other  process. 

The  field  of  dynamic  scheduling  problems  (see  Section  1.3)  was  not 
explored.  The  current  model  could  be  made  to  simulate  a  continual  flow 
of  jobs  into  the  shop  by  stopping  the  process  at  suitable  time  periods, 
injecting  a  number  of  jobs  generated  in  a  specified  manner,  resorting  the 
roster  and  starting  over.  The  major  trouble  Conway  found  with  the 
shortest  operation  discipline  would  now  occur;  i.e.,  large  jobs  would 
have  very  long  waiting  times.  It  would  not  have  the  very  strong 
superiority  in  optimizing  mean  start  time  which  was  observed  here. 

The  capacity  of  the  shop  could  change  over  time,  and  it  could  also 
run  for  only  a  fixed  number  of  time  periods.  An  increasing  shop  could 
be  handled  without  modification  while  for  a  decreasing  shop,  the  main 
scheduling  model  would  have  to  be  modified  to  stop  a  job  of  greatei 
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capacity  requirements  than  the  new  shop  capacity  being  scheduled  too 
close  to  the  change.  The  model  is  not  set  up  to  deal  with  a  shop  of 
fixed  restrictive  duration.  An  integer  programming  formulation  may  be 
required  for  problems  involving  shops  of  decreasing  capacity  or  restricted 
duration . 

A  situation  could  occur  where  it  is  possible  to  "crash"  a  job  at  a 
certain  expense.  The  expense  could  be  such  that  if  say  a  capacity  =  5, 
duration  =  5,  job  is  crashed,  it  may  require  capacity  =  7,  duration  =  4; 
the  expense  would  be  28  -  25  =  3  capacity  time  units.  As  with  the 
restrictive  shop,  the  model  is  not  equipped  to  handle  this.  This  would 
probably  be  an  interesting  stud)7  as  "crashing"  is  frequent  in  practice; 
and  some  simple  rules  given  the  penalties  for  crashing  may  be  developed 
based  on  simple  conditions  existing  at  the  time  a  decision  is  to  be  made. 
Profitability  of  expanding  or  shutting  down  part  of  the  shop  could  also 
be  explored. 

The  work  done  in  this  paper  has  shown  that  simple  rules  can  be 
applied  to  a  job  shop  scheduling  problem  where  the  following  conditions 
apply: 

a.  The  shop  is  of  fixed  capacity  and  nonrestr ict ive  in  duration. 

b.  The  entire  roster  is  fixed  and  known  at  the  beginning  of  the 
scheduling  operation. 

c.  The  jobs  are  of  unalterable  capacity  and  duration. 

These  conditions  are  probably  too  strict  for  the  model  to  have  any 
direct  practical  importance.  However,  they  should  indicate  the  direction 
for  further  research  to  follow  by  giving  future  researchers  a  "feel"  for 
what  may  be  the  best  rule  when  the  restrictions  are  lifted. 


BIBLIOGRAPHY 


Adams,  W.  S.,  An  Introduction  to  APL/360,  Department  of  Computing 
Science,  University  of  Alberta,  Publication  No.  7, 
September,  1967. 

Balas,  Egon,  "An  Additive  Algorithm  for  Solving  Linear  Programs 

With  Zero-One  Variables,"  Operations  Research,  Vol.  13, 
pp.  516-549. 

Brooks,  G.  H.,  and  White,  C.  R.,  "An  Algorithm  for  Finding  Optimal 
or  Near-Optimal  Solutions  to  the  Production  Scheduling 
Problem,"  Journal  of  Industrial  Engineering,  Vol.  16, 
pp.  34-40. 

Conway,  R.  W.,  "Priority  Dispatching  and  Work  in  Process  Inventory 
in  a  Job  Shop,"  Journal  of  Industrial  Engineering. 

Vol.  16,  pp.  123-130. 

Fabrycky,  W.  J. ,  and  Shamblyn,  J.  E.,  "A  Probability  Based 

Sequency  Algorithm,"  Journal  of  Industrial  Engineering. 
Vol.  17,  pp.  308-312. 

Freund,  J.  E.,  and  Williams,  F.  J.,  Modern  Business  Statistics, 
Prentice-Hall,  Inc.,  Englewood  Cliffs,  N.  J.,  1958. 

Richmond,  S.  B„,  Statistical  Analysis,  Second  Edition,  The  Ronald 


Press  Company,  New  York,  1964* 


. 


A-  1 


APPENDIX  A 


A~  1  Data  Generators  and  Processors 
A-la.  Generators 

VGENFLTIU]  V 

V  V+MAX  GENFLT  NUM\K\I 

Cl]  7^0p0 

[2]  0 

[3]  GENFLT 1 :K+?MAX 

[4]  V+VtK 

[5]  I+I+l 

[6]  •>(  I<NUM )  /  GENFLT1 

V 


S/GENPOISSON  [Q]V 

V  V+LAMBDA  GENPOISSON  N;I 

[1]  F^OpO 

[2]  J«-  0 

[3]  GENP0ISS0N1 :+(N<I+I+l ) /GENPOISSON^ 

[4]  £M  ?1 0000  0  )vl  00001 

[5]  A+B+l 

[6]  K+ 0 

[7]  F«-* *( -LAMBDA) 

[8]  GENPOIS SON 3  :+(  ( F*A ) >U )/ GENPOISSON 2 

[9]  iW+P+BxLdtfBZM  1 

[10]  -+GENP0ISS0N3 

[11]  GENP0ISS0N2 : V+VtK 

[12]  -*S£WP0ISS0W1 

[13]  GENPOISSON^  :  ->-0 

V 


Documentation 


*  A  vector  V  is  generated.  The  elements  of  V  have  either  a 
uniform  or  Poisson  distribution 

*  GENFLT  specifies  a  uniform  distribution  with  range  of  one 
to  MAX.  The  order  of  the  vector  is  NUM 

*  GENPOISSON  specifies  a  Poisson  distribution  with  mean  LAMBDA. 
The  order  of  the  vector  is  N 


1  •  ■ 
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A  -  l.b  Processors 

V SUMMARYF 3[Q] V 

V  SUMMARYF 3 

[1]  N+l 

[2]  SUMRF  3«-(10  9)p0.. 

[3]  SUMMARYF31 :N 

[4]  SUMRF3LN;  1+R0SCF31N-,  ~\CHECK  R0SDF3[N  ;  ] 

[5]  iM7+l 

[6]  +(  /!/  =  1 1  )  /  0 
[7]  -+SUMMARYF3 1 

V 


V5™i?yf4[mv 

V  SUMMARYF* * 

[  1  ]  N+-1 

[2]  SUMRF*  +  (10  9 ) p  0 

[3]  SUMMARYF*! :R 

[4]  S£/Mi?F4[F;  >i?0SCF4  [  ;  ^CHECK  ROSDF*lN  ;] 

[5]  +  l 

[6]  -►(#  =  ll)/0 

[7]  + SUMMARYF 41 

V 


V  SUMMARYF 5 [ □] V 

V  SUMMARYF 5 

[1]  tf«-l 

[2]  S’[/Mi?F5^(10  9  )p  0 

[3]  SUMMARYFbliN 

[4]  S£/Mi?F5[tf ;  ]«-i?0SCF5  [  W ;  1CHECK  ROSDF  5[F;] 

[5]  iM7+l 

[6]  +  ( /7=1 1  )  /0 

[7]  +SUMMARYF 51 

V 


Documentation 


*  Computes  summary  statistics  for  the  uniform  rosters 

*  Requires  input  matrices  ROS  CF  and  ROS  DF  representing  capac¬ 
ities  and  durations  of  the  jobs  in  the  rosters 

*  Generates  an  output  matrix  SUM  RF  consisting  of  the  summary 
characteristics  of  the  input  rosters 


‘ 


VSUMMARYP3 [□] V 

V  SUMMARYP 3 

[1]  N+l 

[2]  SUMRP 3«-(10  9 ) p  0 

[3]  SUMMARYP31 :SUMRP3lN ; l+ROSCPZlN; 1CHECK  ROSDP3[N ; ] 

[4]  /M/+1 

[5]  ->(  71^  =  1 1  ) /0 

[6]  +SUMMARYP3 1 

V 


V  SUMMARYP^ [ □] V 

V  SUMMARYP 4 

[1] 

[2]  SUMRPU,  +  (  10  9  )  p  0 

[3] 

[4]  S£/m4PJP41  -.SUMRPHlN ;  ]«-PtfSCP4  [  P  ;  IPPPCK  P<9SPP4[W:] 

[5]  m/+i 

[6]  -KP=11  )/0 

[7]  -+ SUMMARYP  41 

V 


1  SUMMARYP 5[[]]  V 

V  SUMMARYPb 

Cl] 

[2]  SUMRP5-<r(  1  0  9  )  p  0 

[3]  SUMMARYP51:N 

[4]  5PMPP5CP;  >P£SCP5[W;  ~\CHECK  ROSDPStN ;] 

[5] 

[6]  +(ff=ll)/0 

[7]  +SUMMARYP 51 

V 

Documentation 


*  Computes  summary  statistics  for  the  Poisson  rosters 

*  Requires  input  matrices  ROS  CP  and  ROS  DP  representing 
capacities  and  durations  of  the  jobs  in  the  rosters 

*  Generates  an  output  matrix  SUMRP  consisting  of  the  summar 
characteristics  of  the  input  rosters 


- 
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VCHECKiniV 

V  A+X  CHECK  P 

[1]  4«-9p0 

[2]  4[  1  ]«-(  +/X)  ipX 

[3]  VECTAVCAP+(pX)p(Alll ) 

[4]  4  [  2  ]  -<-4  RND(((+/(X-VECTAVCAP)* *2)ipX)*0.5) 

[5]  4[3>[/J 

[6]  4[4]«-r/Z 

[7]  4[5>(  +/P)tPP 

[8]  7PC,!Z,4ra?C>^(pP)p  (4[5]  ) 

[9]  4[6>4  P/l/P(((+/(P-7PCT47PPC>)*2)*pP)*0.5) 

[10]  4[7]«-L/P 

[11]  4[8>r/P 

[12]  4  [  9  ]*<-  +  /  (  PxJ  ) 

[13]  (  •  »  ) 

[14]  (' CAPACITIES  ARE  »;*;) 

[15]  'FREQUENCIES' 

[16]  Y4PFPF3  J 

[17]  (»») 

[18]  ('DURATIONS  ARE  »;P;) 

[19]  'FREQUENCIES' 

[20]  TAB FREQ  P 

[21]  ( » • ) 

[22]  4 

[23]  (’») 

[24]  (  »  *  ) 

[25]  (  »  »  ) 

V 

Documentation 


*  Incorporated  in  the  SUMMARY  programs  to  perform  the  actual 
computation  of  the  summary  characteristics 


*  Input  is  two  vectors,  X  and  P.  X  is  the  vector  of  capacities  of  the 
jobs  in  the  roster,  and  P  is  the  duration  vector 

*  Output  is  a  vector  A  consisting  of  the  nine  summary  results 

*  Displays  the  vector  A  of  summary  characteristics 

*  Incorporates  the  routine  TABFREQ  to  display  a  frequency 
distribution  of  the  capacities  and  durations  of  the  jobs  in  each 
roster 


.  .  > '  ~ 
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VTABFREQIUI V 

V  TABFREQ  R;S;I 
Cl]  1+ 1 

[2]  £«-100p0 

[3]  T4BFf?£’Ql:S'[I>+/i?  =  l 
C4]  J.5CJ] 

[5]  +((+/S)=pi?)/0 

[6]  I+-I+1 

[7]  +TABFREQ1 

V 

Documentation 


*  Input  is  vector  R 

*  Computes  vector  S,  where  S  m  is  the  number  of  elements  of  R 

equal  to  I 

*  Displays  I  and  S  Q  I]] 


V STRINGOUT [□] V 

V  STRINGOUT  ROS 

[1]  VECT+ OpO 

[2]  N+-1 

[3]  STRINGOUT 1 : VECT+VECT ,ROSlN ; ] 

[4]  N+N  +  l 

[5]  -* *(/i  = ll)/0 

[6]  +STRINGOUT 1 

V 

Documentation 


*  Input  is  a  10  by  N  matrix  ROS 

*  Output  is  a  vector  VECT  consisting  of  the  rows  of  ROS 


Vi? /!/£>[[]]  V 

V  R+N  RND  X 

[1]  i?«-(  10*-tf)xL  0.  5+ZxlO*/7 

V 


Documentation 


*  Library  defined  rounding  function 

*  R  is  a  number  equal  to  X  rounded  to  N  places  of  decimal 


A-  2  Main  Scheduling  Programs 


A- 2a.  Main  Routine 

VMAINl □] V 

V  MAIN 

[1]  MAT+{ 14  9 ) p  0 

[2]  ' ENTER  DURATIONS  CAPACITIES  AND  NUMBER  OF  JOBS' 

[3]  INPDUR*- □ 

[4]  INPCAP+U 

[5]  INPNRJOB+U 

[6]  * SPECIFY  SHOP ’ 

[7]  SHOP*-  □ 

[8]  RANDOMS 

[9]  MATlhl+SUM 

[10]  RANDOMFITS 

[11]  MATL2;1*-SUM 

[12]  (••) 

[13]  SSRANDOMFIT 

[14]  MATL3  O+SUM 

[15]  LSRANDOMFIT 

[16]  MATH;1*-SUM 

[17]  SPFITS 

[18]  MATl5il+SUM 

[19]  LPFITS 

[20]  MAT16;  ~\*-SUM 

[21]  (”) 

[22]  SPSCFITS 

[23]  MATH  ;  ~\*-SUM 

[24]  SPLCFITS 

[25]  MATlBil+SUM 

[26]  LPSCFITS 

[27]  MATlB\~\*rSUM 

[28]  LPLCFITS 

[29]  MATt  1 0  ;  "]*-SUM 

[30]  (»») 

[31]  SSSPFIT 

[32]  MATlll\]+SUM 

[33]  SSLPFIT 

[34]  MAT112;1+SUM 

[35]  LSSPFIT 

[36]  MAT113 il+SUM 

[37]  LSLPFIT 

[38]  MATHml+SUM 

[39]  STATS 

V 
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Documentation  (for  MAIN  on  Page  6) 

*  This  is  the  main  scheduling  program 

*  Accepts  definition  of  the  roster  and  shop  from  the  keyboard 
and  incorporates  the  scheduling  routines.  The  following  is 
automatic  when  the  program  MAIN  is  called. 

*  Displays  a  0  and  unlocks  the  keyboard  four  times.  The 
following  are  defined  by  the  operator  in  sequence: 

INPDUR  -  roster  durations  vector 
INPCAP  -  roster  capacities  vector 
INPNRJOB  -  number  of  jobs  in  roster 
SHOP  -  shop  capacity  vector 

*  Executes  each  of  the  fourteen  scheduling  rules  in  turn  with 
these  global  variables  as  input 

*  Stores  the  summary  characteristics  calculated  when  each 
rule  is  implemented  in  the  matrix  MAT 

*  Incorporates  a  program  STATS  which  ranks  the  columns  of 
MAT  and  prints  the  ranked  summary  matrix 


, 

' 

' 

i 
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A-  l.b  Scheduling  Programs 


1RAUD0MSI D]V 

V  RANDOMS 
Cl]  P+-INPDUR 

[2]  X+INPCAP 

[3]  NRJOB+INPNRJOB 

[4]  C+-SHOP 

[5]  ASGNXT 

[6]  SUMARY 

V 


VRANDOMFITS CD]V 

V  RAN DOME  ITS 

[1]  P+INPDUR 

[2]  X+INPCAP 

[3]  NRJOB+INPNRJOB 

[4]  C+SHOP 

[5]  ASGNAHD 

[6]  SUMARY 

V 


Documentation 


*  Specifies  the  roster  and  shop  as  the  global  variables  defined  in  MAIN 

*  Incorporates  the  assignment  routines  ASGNXT  for  RANDOMS  and 
ASGNAHD  for  RANDOMFITS 


*  Incorporates  the  summary  statistics  generating  routine  SUMARY 


VSSSPFITIU1V 

V  SSSPFIT 

[1]  CAPLIN PC AP 

[2]  DURBIN PD  UR 

[3]  NRJOB+INPNRJOB 

[4]  ASSORT 

[5]  DUR+AA 

[6]  CAP+-BB 

[7]  ASSORTS 

[8]  X+-BB 
[  9  ]  P+AA 

[10]  C+-SH0P 

[11]  ASGNAHD 

[12]  SUM ARY 


V 


i  . 
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VSSRANDOMFITLUlV 

V  SSRANDOMFIT 

Cl]  CAP+IN PC AP 

[2]  DUR+-INPDUR 

[3]  NRJOB+INPNRJOB 
C4]  ASSORTS 

[  5  ]  X+-BB 

[6]  P+-AA 

[7]  C+SHOP 

[8]  ASGNAHD 

[9]  SUMARY 

V 


USSLPFIT [D]V 

V  SSLPFIT 

[1]  CAPLIN PC AP 

[2]  DURBIN PD  UR 

[3]  IIRJ  OB+-INPNRJOB 

[4]  DEC  SORT 

[5]  DUR+AA 
t  6 ]  CAP+BB 

[7]  ASSORTS 

[8]  X+BB 

[9]  P+AA 

[10]  C+SHOP 

[11]  ASGNAHD 

[12]  SUMARY 

V 


Documentation 


*  Scheduling  rules  based  on  priority  to  smallest  size  job 

*  Specifies  the  roster  and  shop  as  the  global  variables  defined  by 
MAIN 

*  Performs  the  tie-breaking  sort  using  ASSORT,  nothing,  and 
DECSORT  to  get  smallest,  random,  or  largest  duration  first 

*  Uses  the  routine  ASSORT  3  to  arrange  the  roster  in  ascending  size 
order 

*  Specifies  the  input  arguments  for  the  scheduling  routine  ASGNAHD 

*  Incorporates  ASGNAHD  and  SUMARY  to  compute  the  summary 
statistic  s 


VLSSPFITlUlV 

V  LSSPFIT 

[1]  CAP+IN PC AP 

[2]  DURBIN PD UR 

[3]  NRJOB+INPNRJOB 

[4]  ASSORT 

[5]  DVR+AA 

[6]  CAP+-BB 

[7]  DECS0RT3 

[8]  X+BB 

[9]  P+AA 

[10]  C+-SHOP 

[11]  ASGNAHD 

[12]  SUMARY 

V 


VLSRANDOMFITlUlV 

V  LSRANDOMFIT 

[1]  CAP+INPCAP 

[2]  DUR+INPDUR 

[  3  ]  NRJ OB+-INPNRJ OB 

[4]  DEC  SORT  3 

[  5  ]  X+BB 

[6]  P+AA 

[7]  C+SHOP 

[8]  ASGNAHD 

[9]  SUMARY 

V 


VLSLPFITl D]V 
V  LSLPFIT 

[1]  CAP+-INPCAP 

[2]  DUR+INPDUR 

[3]  NRJOB+INPNRJOB 

[4]  DECSORT 

[5]  DUR+AA 

[6]  CAP+-BB 

[7]  DECSORT 3 
[  8  ]  X+-BB 

[ 9 ]  P+AA 

[10]  C+SHOP 

[11]  ASGNAHD 

[12]  SUMARY 


V 
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Documentation 


*  Scheduling  rules  based  on  priority  to  largest  size  job 

*  Specifies  the  roster  and  shop  as  the  global  veMables  defined  by 

MAIN 

*  Performs  the  tie -breaking  sort  using  ASSORT,  nothing,  and 

DECSORT  to  get  smallest,  random,  or  largest  duration  first 

*  Uses  the  routine  DECSORT  3  to  arrange  the  roster  in  descending 

size  order 

*  Specifies  the  input  arguments  for  the  scheduling  routine  ASGNAHD 

*  Incorporates  ASGNAHD  and  SUMARY  to  compute  the  summary 

statistic  s 


USPSCFITS CD]V 

V  SPSCFITS 

[1]  CAP+INPDUR 

[2]  DUR+INPCAP 

[3]  NRJOB+INPNRJOB 

[4]  ASSORT 

[5]  DUR+BB 

[6]  CAP+AA 

[7]  ASSORT 

[8]  X+BB 

[9]  P+AA 

[10]  C+SHOP 

[11]  ASGNAHD 

[12]  SUM ARY 

V 


VSPFITSlW]  V 
V  SPFITS 

[1]  DUR+-INPDUR 

[2]  CAP+INPCAP 

[3]  NRJOB+INPNRJOB 

[4]  ASSORT 

[5]  P+AA 

[6]  X+-BB 

[7]  C+SHOP 

[8]  ASGNAHD 

[9]  SUMARY 


V 


, 
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VSPLCFITSIW]  V 

V  SPLCFITS 

[1]  CAP+-INPDUR 

[2]  DUR+INPCAP 

[3]  NRJ  OB+-INPNRJ  OB 

[4]  DECSORT 
C  5 ]  DUR+BB 

[6]  CAP* A  A 

[7]  ASSORT 

[8]  X+BB 

[9]  P+AA 

[10]  C+-SROP 

[11]  ASGNAHD 

[12]  SUMARY 

V 

Documentation 


*  Scheduling  rules  based  on  priority  to  the  smallest  duration  job 


*  Specifies  the  roster  and  shop  as  the  global  variables  defined  by  MAIN 


*  Performs  the  tie-breaking  sort  using  ASSORT,  nothing  and  DECSORT 
to  get  the  smallest,  random  or  largest  capacity  first 


*  Uses  the  routine  ASSORT  to  arrange  the  roster  in  ascending 
duration  order 


*  Specifies  the  input  arguments  for  the  scheduling  routine  ASGNAHD 


*  Incorporates  ASGNAHD  and  SUMARY  to  compute  the  summary 
statistic  s 


VLPSCFITSIUIV 
V  LPSCFITS 

[1]  CAP+-INPDUR 

[2]  DUR+INPCAP 

[3]  NRJ  0B+-INPNRJ  OB 

[4]  ASSORT 

[5]  DUR+-BB 

[6]  CAP* A  A 

[7]  DECSORT 

[8]  X+BB 

[9]  P+AA 

[10]  C+-SH0P 

[11]  ASGNAHD 

[12]  SUM ARY 


V 
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VLPFITSIU1V 

V  LPFITS 

[1]  DUR+INPDUR 

[2]  CAP+INPCAP 

[ 3 ]  NRJOB+INPNRJOB 

[4]  DECSORT 

[5]  P+AA 

[6]  X+BB 

[7]  C+SHOP 

[8]  ASGNAHD 

[9]  SUMARY 

V 


VLPLCFITSlUlV 

V  LPLCFITS 

Cl]  CAP+-INPDUR 

[2]  DUR+-INPCAP 

[3]  NRJOB+INPNRJOB 

[4]  DECSORT 

[5]  DUR+BB 

[6]  CAP+AA 

[7]  DECSORT 

[8]  X+BB 

[9]  P+AA 

[10]  C+-SHOP 

[11]  ASGNAHD 

[12]  SUMARY 

V 


Documentation 


*  Scheduling  rules  based  on  priority  to  the  largest  duration  job 

*  Specifies  the  roster  and  shop  as  the  global  variables  defined  by 
MAIN 


*  Perform  the  tie- breaking  sort  using  ASSORT,  nothing  and 
DECSORT  to  get  the  smallest,  random  or  largest  size  first 

*  Uses  the  routine  DECSORT  to  arrange  the  roster  in  descending 
duration  order 

*  Specifies  the  input  arguments  for  the  scheduling  routine  ASGNAHD 

*  Incorporates  ASGNAHD  and  SUMARY  to  compute  the  summary 
statistic  s 


N 
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A  -  2c)  Dependent  Sorter,  Scheduler  and  Summary  Routines 
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A  -  2c  -  I)  Sorter  Routines 

VDECSORTIU l]V 

V  DECSORT 

[1]  AA* *-BB+0p  0 

[2]  A+DUR 

[3]  B+CAP 

[4]  DECS0RT1 : C+\ / A 

[5]  •>(  C  =  0  )  /DECS0RT2 

[ 6 ]  D+A i C 

[7]  AA+-AA.C 

[8]  BB+BB.BID] 

[9]  4[£>0 

[10]  +DECS0RT 1 

[11]  DECSORT 2:-*0 

V 


VASS0RTIU1V 
7  ASSORT 

[1]  AA+BB+ OpO 

[2]  4«-Z?£/i? 

[3] 

[4]  ASS0RT1 \C+l/A 

[5]  +(C>U)0QQ) /ASS0RT2 
[  6  ]  D+A  i  C 

[7]  AA+AA, C 

[8]  BB+BB,BID1 

[9]  4[Z?]^1000’t)0000 

[10]  ^ASSORT  1 

[11]  4SS0PT2  :-*0 
V 


Documentation 


*  Input  is  vectors  DUR  and  CAP 


*  Program  specifies  a  vector  AA  corresponding  to  DUR  and  a 
vector  BB  corresponding  to  CAP 

*  AA  is  the  vector  DUR  rearranged  in  descending  order  by 
DECSORT  and  ascending  order  by  ASSORT 

*  If  the  element  DUR  (ZlH  becomes  AA  Q J I]  ,  then  the  element 
CAP  Cll  will  become  BB  (UJ  for  all  elements 


- 


' 
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VDECS0RT3\_W]V 

V  DEC SORT 3 

[1]  AA"* *~BB +-S+- 0 p  0 

[2]  A+DUR 

[3]  B+CAP 

[4]  SIZE+AxB 

[5]  DECSORT32 :X+[ /SIZE 

[6]  ->(  X=0)  /DECS0RT31 

[7]  IND+SIZExX 

[8]  S+-S  tSIZElLND] 

[9]  AA+AA,AlIND] 

[10]  BB+BB ,B[ INDJ 

[11]  SIZE[IND]+0 

[12]  -+DECSORT  3  2 

[13]  DECSORT 31:+0 

V 


V4S,S0i?T3[[]]V 
v  4£so;?T3 

[1]  0  p  0 

[2]  A+DUR 

[3]  B+CAP 

[4]  SIZE+AxB 

[5]  ASSORT32 : X+L /SIZE 

[6]  +( Z=10000 )/ASSORT31 

[7]  IND+SIZE \X 

[8]  S+S ,SIZEIIND] 

[9]  AA+AA^IIND] 

[10]  BB+BB  ,BiIND~\ 

[11]  SIZElINDl+lOOOO 

[12]  ^ASSORT  3  2 

[13]  ASSORTS l:+0 
V 


Documentation 

*  Input  arguments  are  vectors  DUR  and  CAP 

*  Generates  a  vector  SIZE  DD  =  dur  Cil  x  cap  Ci3 
DECSORT  3  sorts  SIZE  in  decending  order  and  ASSORT  3 
sorts  SIZE  in  ascending  order 

*  Output  arguments  are  vectors  S,  AA,  and  BB  corresponding  to 
SIZE,  DUR  and  CAP  respectively 

*  If  SIZE  Cil  becomes  S  CjD  ,  DUR  ClU  and  CAP  ClU 
become  AA  QjJ  and  BB  CJ]  respectively  for  all  elements 
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A-2C-II  Summary  Routines 

VSUMARYIU1V 

V  SUMARY 

[1]  T++/(PxX) 

[2]  U+CC* *\ /F 

[3]  EFF+-TtU 

[4]  MS<r(+/NUM)*pNUM 

[5]  SUM* 9p0 

[6]  SUMlll+MS 

[7]  Si/M[2>[/F 

[8]  S£/M[3]«-( +/JOB) v(  \ /F) 

[9]  SUMlUl+EFF 

[10]  BEGIN<-(  (T/F)-(  4  |  T/F))*4 

[11]  SLACK+  +  /ClBEGIN  +  \  2*BEGIN'} 

[12]  SUMl  S]+-l-SLACKiCCx2*BEGIN 

[13]  S  UMl6  ]«-+  /  Px  Jx  (  RUM-1  ) 

[14]  SUM171++/NUM 

[15]  VE  CTMS •«-(  p  NUM )  p  MS’ 

[16]  SUMl  8>(  (+/(NUM-VECTMS)*2)ipNUM)*0.  5 

[17]  S£/M[9X(  +/NUM*2)ipNUM)*Q.  5 

[18]  S£/M 

V 


Documentation 


Input  arguments  required  are  as  follows: 

NUM  =  vector  of  start  times 
F  =  vector  of  finish  times 
P  =  vector  of  job  durations 

X  =  vector  of  job  capacity  requirements 

CC  =  capacity  of  a  uniform  shop  at  start 
JOB  =  number  of  jobs  running  each  time  period 
C  =  unused  shop  capacity  each  time  period 


*  Computes  a  vector  SUM  consisting  of  the  nine  summary 
characteristics  described  in  Section  4.  2.  2.  CII  * 


*  Displays  the  vector  SUM 


'  X 
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A  -  2C.  Ill  Scheduling  Routines 

VASGNXTlUlV 

V  ASGNXT 
Cl]  HEADING 

[2]  AVCAP+(+/X)ipX 

[3]  JOB*'  5  0  0  p  0 

[4]  N+J+ 1 

[5]  CC+C [1] 

[6]  NUM+F+NRJOBpO 

[7]  ASGNXT1 : NN+N 

[  8  ]  ASGNXT 2  :  -►(  Cl  N]  <Xl  </] ) /ASGNXT 5 

[9]  i45{7ra3:C'[^]^C'[^]-Z[^] 

[10]  e70B[tf]«-«T0B[ff]  +  l 

[11]  M/  +  1 

[12]  -►(  N<NN+PlJ'}  )  /ASGNXT3 

[13]  N+NN 

[14]  FlJl+N+PlJl-1 

[15]  NUMIJ1+N 

[16]  DISPLAY 

[17]  J+J  +  l 

[18]  -►(  J >NRJ OB  ) /ASGNXT1* 

[19]  -+ASGNXT2 

[20]  ASGNXT5 : N+N+l 

[21]  +ASGNXT 1 

[22]  ASGNXT^  i  (  •  » 

[23]  ASGNDETL 

V 

Documentation 

*  Input :  X  = 

P 
C 

NR JOB = 

*  Output:  AVCAP= 

NUM  = 

F 

JOB 
C 


) 


vector  of  job  capacities 
vector  of  job  durations 
shop  capacity  in  each  time  period 
number  of  jobs  in  the  roster 

average  job  capacity 
vector  of  start  times 
vector  of  finish  times 

number  of  jobs  running  in  each  time  period 
unused  shop  capacity  in  each  time  period 


*  Jobs  are  assigned  in  order  of  occurance  in  X  and  P  as  described 
in  Section  4.  2.  2.  CIII 


*  A  complete  display  of  the  schedule  is  made.  To  avoid  this 

modify  statements  [1]  and  [Ll6D  to— ^2  and— >17  respectively. 
The  summary  statistics  only  will  then  be  printed 


. 

. 
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VASGNAHDIWM 

V  ASGNAHD 

[1]  HEADING 

[2]  AVCAP+(+/X)ipX 

[3]  CC* *-C  [1] 

[4]  JOB*-  5  0  0  p  0 

[5]  N*-J*-l 

C  6 ]  NUM+M+F+NRJOB p 0 

[7]  AS GNAHD12  :  NN+-N 

C8]  +(CIN]<XIJ])/ASGNAHD 14 

[9]  ASGNAHD13:ClN]*-ClN]-XlJ] 

[ 10]  JOBlNl+JOBlNl+1 

[11]  N+N+ 1 

[12]  +(N<NN+PlJ]) /ASGNAHD  13 

[13] 

[14]  MIJ]*- 1 

[15]  F[J]^  +  P[e7]-l 

[16]  NUM  [  e7  ] 

[17]  DISPLAY 

[18]  4SG/IMtfZ?14  :J+J  +  1 

[19]  -*(  J>NRJOB  )  / ASGNAHD15 

[20]  +  (M[e7G=l  )/4SG/lM#P14 

[21]  -MSGMPZm 

[  22]  ASGNAHD  15  :+(  [  /M  =  l  ) /4£GW4tfZ?l  6 

[23]  <7«-l 

[24]  N+N+l 

[25]  +ASGNAHD1* 

[ 26 ]  ASGNAHD16 : (  »  »  ) 

£27]  ASGNDETL 

V 


Documentation 


*  Input  and  output  arguments  have  identical  labels  to  those  of 
ASGNXT  and,  therefore,  this  program  is  interchangeable  with  it  in 
the  scheduling  routines 

*  Assignment  is  attempted  in  the  order  the  jobs  appear  in  X  and  P, 
but  a  search  ahead  for  a  job  which  does  not  violate  capacity- 
restrictions  is  made  as  described  in  Section  4.  2.  2.  CIII 

*  A  complete  display  of  the  schedule  is  made.  To  avoide  this,  modify 
statements  to  and  Cl  73  to  — >2  and  -^18  respectively. 

*  The  schedule  shows  job  duration,  capacity,  start  time,  and  finish 
time  in  that  order.  If  the  option  is  taken,  only  the  summary 
characteristics  of  the  schedule  produced  by  SUMMARY  are 
displayed 


. 
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1ASGNDETLIW}  V 

V  ASGNDETL 
Cl]  OCC+\/F 

[2]  ('JOBS  IN  PROCESS  ' ;JOBl \0CC1 ; ) 

[3]  (» UNUSED  CAPACITY  ' ;Cl\OCC];) 

V 


Documentation 

*  Input  vectors  are:  F  -  job  finish  times,  JOB  -  number  of 
jobs  in  process  and  C  -  unused  capacity 

*  Displays  the  npmber  of  jobs  in  process  and  unused  capacity 
for  the  time  the  roster  is  in  the  shop 

A- 2a.  Summary  Ranking  Routines 


V  STATS' CDJV 

V  STATS 

[  1  ]  N<- 1 

[2]  STATS2  :MAX+\ /MAT [  ; N ] 

[3]  MIN+l/MATl  ;Nl 

[4]  DIFF+MAX-MIN 

[5]  J4- 1 

[6]  STATS1  :MATII  ;N1  +  (MATII  ;N']-MIN)* *DIFF 

[7]  I«-I  +  l 

[8]  ->(  J<15  j/SJVlTSl 

[9]  N<-N+l 

[10]  -^ ( iY <  1  0  )/STATS2 

[11]  M4T+10000x(4  /M?) 

[12]  MAT 

V 


Documentation 

*  Input  is  matrix  MAT  of  summary  statistics  as  generated  by 
MAIN 


*  Performs  the  column  by  column  ranking  of  MAT  described 
in  Section  4.  2.  2.  e 


*  Output  is  a  matrix  MAT  consisting  of  the  rankings  in  place 
of  the  summary  statistics 


A  -  3  SUMMARY  COMPUTATION  PROGRAMS 


A  -  3a  Summary  Generators 


VASYN CD]V 

V  ASYN 

[1]  ATOTALH  14  6  0 ) p  0 

[2]  MSDATOTAL+{  14  3)p0 

[3]  N+l 

[4]  AS  YN1  lATOTALlNi  1+AF3  [/!/;]  ,4F4[tf ;  ]  tAF5lN;  ]  tAP3lN\  ] 
i4P4[tf;  ]  ,;4P5[W;  ] 

[5]  ]«-MSD  ATOTAL  [  N  ;  ] 

[6]  (  "  ) 

C7]  (  ”  ) 

[8]  (  •  »  ) 

C  9  ]  (  '  »  ) 

[10]  ->(//  =  14)/0 

[11]  N+N+ 1 

[12]  -MSYtfl 

V 


V3S7/'/[D]  V 
V  BSYN 

[1]  BT0TAL<-{  14  6  0  )  p  0 

[2]  MSDBTOTALH  14  3)p0 

[  3  ]  N+- 1 

[4]  BSTtfl  :BT02\4L[B;  >BF3[W;  ]  ,BF4  [A7;  ]  ,BF5  [/!/;]  ,BP3[B;  ] 
BP4[B; ] ,BP5[ff ; ] 

[5]  MSDBTOTALIN;  ~\+MSD  BTOTALlN ;] 

[6]  (  »  •  ) 

[7]  (”) 

[8]  (  ”  ) 

[9]  (  •  »  ) 

[10]  ->(  77  =  1 4  )/0 

[11]  N+N+l 

[12]  +BSYN 1 


V 


•  c 
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VC5Y/y[  □  ]  V 
V  CSYN 


Cl] 

CTOT AL+-(  1 4  6  0  )  p  0 

[2] 

MSDCTOTAL*-{ 14  3)p0 

[3] 

N+l 

[4] 

CSYN1:CT0TALLN  i]+CF3tN ;]  ,CF4[tf;  ],  CF5  [  /!/;],  CP3  [  N  ;]  , 

CT4 [F ; ]  ,  CP5  [  /!/ ;  ] 

[5] 

MSDCTOTALIN ; ]+MSD  CTOTALlN ; 

] 

[6] 

(  '  1  ) 

[7] 

(  1  *  ) 

[8] 

(  '  '  ) 

[9] 

(  ”) 

[10] 

N 

[11] 

->(/\7=14)/0 

[12] 

N+N  +  l 

[13] 

V 

■+CSYN1 

VDSYN [□] V 

V 

DSYN 

[1] 

DTOT  AL+-(  14  6  0  )  p  0 

[2] 

MSDDTOTAL+( 14  3)p0 

[3] 

N+ 1 

[4] 

DSYN1  iDTOTALlN;  ]«-DF3[/17;  ]  ,FF4[/I/;  ]  ,  FF5  [  F  ;  ]  ,  DP3  [  N  ;  ]  , 
FP4[F;  ]  ,FP5[F;  ] 

[5] 

N 

[6] 

MS DDTOTAL  [  F ;  ]  +-MSD  DTOTAL  [  N  ; 

] 

[7] 

(  ’  '  ) 

[8] 

(  ’  ’  ) 

[9] 

(  '  ’  ) 

[10] 

(  '  ‘  ) 

[11] 

O 

J" 

tH 

II 

N— ' 

\ 

[12] 

N+N  +  l 

[13] 

V 

•+DSYN1 

Program  Name  Input 

Output 

ASYN 

AF3,  AF4,  AF 5,  AP3, 

AP4,  AP5  MSDATOTAL 

BSYN 

BF3,  BF4,  BF 5,  BP3, 

BP4,  BP5  MSDBTOTAL 

CSYN 

CF3,  CF4,  CF5,  CP3, 

CP4,  CP5  MSDC TOTAL 

DSYN 

DF3,  DF4,  DF5,  DP3, 

DP4,  DP5  MSDD TOTAL 

*  Input  is  rankings  for  each  roster  according  to  each  measure  of 
effec  tivene  ss 


*  Output  is  mean  and  standard  deviation  for  each  measure  of 
effec  tivene  s  s 


. 
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VMSDIUIV 

V  MSDA+MSD  AiMEAN \STDEV \AA\I \S 

[1]  MEAE+{+/A)* *pA 

[2]  STDEV+( (+/ (A~(pA)pMEAN)*2 )tpA)* 0.5 

[3]  MSDA+( pA ) %MEAN ,STDEV 

[4]  U+MSDA 

[5]  ( »• ) 

C  6  ]  AA*-(  pA  )  p  0 

[7]  1+1 

[8]  MSD2:S+l/A 

[9]  -*•(  5  =  10  0  00  00  )  /MSD3 

[10]  AA[_I~\ 

[11]  ALA \S]+1 000000 

[12]  I«-I+l 

[13]  +MSD 2 

[14]  MSD3:U+AA 

V 


Documentation 


*  Required  input  is  a  vector  A. 

*  Computes  the  order,  mean  and  standard  deviation  of  A  and  places 
the  result  in  a  vector  MSDA 


*  Displays  MSDA 

*  Arranges  the  vector  A  in  ascending  order  in  AA 

*  Displays  A  A 

*  This  program  is  incorporated  in  the  programs  (A,  B,  C,  or  D)  SYN 
which  specify  the  vector  A  as  each  row  of  the  matrix  formed  by 
joining  (A,  B,  C,  or  D)  F3,  F4,  F5,  P3,  P4,  P5.  MSDA  then 
replaces,  column  by  column,  the  zeros  in  the  zero  matrix  MSD 
(A,  B,  C,  or  D)  TOTAL 
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VCOMPARElUl! 

V  COMPARE 

[1]  N+- 1 

[2]  AN-1 

[3]  COMPARE  2 :(’ COMPARING  ’ ; N;'  AND  '\M\) 

[4]  MSDATOTAL  IN  ;  hDIFFMEANS  MSDATOT  AL\_M ;  ] 

[5]  (  »  «  ) 

[6]  M+M+l 

[7]  -►(  Af=l  5  )  / COMPAREl 

[8]  ^COMPARE  2 

[9]  COMPAREl :N+N+1 

[10]  ->(/!/  = 15)/0 

[11]  (  "  ) 

[12]  (  "  ) 

[13]  A/^-yV 

[14]  ^COMPARE  2 

V 


V5IFFAf5A/l/S[[]]V 

V  A  DIFFMEANS  B  ;  N 1 ;  N 2  ;  S 1 ;  S 2  ;  X 1 ;  X 2  ;  S H AT ; SD ;T  ;Z 

[1]  N1+-A11] 

[2]  Zl+A[2] 

[3]  51«-4[3] 

[4]  N2+-B  [  1  ] 

[5] *  *2«-B[2] 

[6]  5  2  <-B  [  3  ] 

[7]  Z+N1+N2-2 

[8]  SHAT+(  (  ( 71/1x51x51 )  +  (  N2*S2*S2  ) ) *Z  )*0  .  5 
[  9]  SD+SHATx( ( ( Nl+N2)*(NlxN2 )  )*0. 5 ) 

[10]  T+{ X1-X2 ) iSD 

[11]  ->(  TZ0)/D IFFMEA NS1 

[12]  T+-T 

[13]  D  IFF  ME  AN  SI :  (  'T  IS  '  ;  T ;  »  J/ITF  *  ;  Z  ;  *  FF’) 

V 


Documentation 


*  DIFFMEANS  performs  the  t  test  of  Section  4.  2.  3  on  input  vectors 
A  and  B  and  displays  the  result 

*  COMPARE  specifies  A  and  B  for  DIFFMEANS  to  be  each  possible 
pairs  of  rows  in  MSD  (A,  B,  C  or  D)  TOTAL 

*  The  output  is  a  statement  of  the  row  numbers  of  MSD  (A,  B,  C  or  D) 
TOTAL  being  compared  and  the  corresponding  t  value 


. 
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Vi?4M[]]V 

V  RANK 

[1]  N* * 1 

[2]  INDEX* OpO 

[3]  P4PZ3  :  7tfCTOi?«-j4F3  [tf ;  ]  ,  4F4  [/!/;]  ,  4F5  [//;],  4P3  [  N  ;  ]  ,  A  P4  [  N  ;  ]  , 
4P5[tf; ] 

[4]  INDEX* OpO 

[ 5 ]  RANK1 : A*l J VECTOR 

[6]  ->(4  =  1000000)/  RANK2 

[7]  IND*9+( VECTORxA ) 

[ 8 ]  INDEX*INDEX , ( ( IND - ( 1 0 | IND ) ) i 1 0 ) 

[9]  VECTORlIND- 9]«-i  000000 

[10]  +RANK1 

[11]  RANK2 :U*INDEX 

[12]  (”) 

[13]  ->(  /1/  =  1 4  )  /  0 

[14]  N*N+ 1 

[15]  ->i?4P^3 

V 

Documentation 


*  Input  vector  is  row  by  row  from  matrices  (A,  B,  C,  or  D)  F3,  F4, 
F5,  P3,  P4,  P5 

*  Each  row  is  rearranged  in  ascending  numerical  order  and  the  index 
each  element  had  in  the  original  vector  is  recorded.  The  indices  are 
computed  so  a  first  number  1  to  6  inclusive  represents  an  element 
from  F3,  F4,  F5,  P3,  P4  or  P5  respectively.  The  second  number 

if  printed  represents  the  position  of  the  element  in  the  matrix 

*  Output  is  a  display  of  these  indices  in  the  order  the  elements  appear 
in  the  rearranged  vector 
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